Physics

Date: Thu, 9 May 2024 | Total: 78

#1 How a table modulates the risk of airborne transmission between facing individuals [PDF] [Copy] [Kimi]

Authors: Oğuzhan Kaplan ; Manouk Abkarian ; Simon Mendez

Airborne transmission has been recognized as an important route of transmission for SARS-CoV-2, the virus responsible for the COVID-19 pandemic. While coughing and sneezing are spectacular sources of production of infected aerosols with far-reaching airflows, the prevalence of asymptomatic transmissions highlighted the importance of social activities. Gathering around a table, a common scenario of human interactions, not only fixes a typical distance between static interlocutors, but influences airborne transmission, by serving both as a flow diverter and a surface for droplet deposition. Here, we use high-fidelity large-eddy simulations to characterize short-range airborne transmission when two people face each other at a typical table. We show that compared to the natural distance travelled by free buoyant puffs and jets, the distance between the table and the emission constitutes a new length scale that modifies downward exhaled flows, common during nose breathing, speech, and laughter. When the table is close to the emitter, its main effect is to restrict the forward spread of emitted particles. However, if the distance between individuals is too short, particles reaching the recipient become more concentrated, rising transmission risks. Additionally, simulations of forceful exhalations, like laughter, demonstrate that the table acts as a filter that collects medium-sized particles that would have remained in the free jet otherwise, but can in that case be involved in the fomite transmission route. The table introduces a cut-off size for particles that depends on the inertia of the exhaled material, thereby modifying the size distribution of particles suspended in the air.

#2 Multimode amplitude squeezing through cascaded nonlinear optical processes [PDF] [Copy] [Kimi]

Authors: Sahil Pontula ; Yannick Salamin ; Charles Roques-Carmes ; Marin Soljacic

Multimode squeezed light is enticing for several applications, from squeezed frequency combs for spectroscopy to signal multiplexing in optical computing. To generate squeezing in multiple frequency modes, optical parametric oscillators have been vital in realizing multimode squeezed vacuum states through second-order nonlinear processes. However, most work has focused on generating multimode squeezed vacua and squeezing in mode superpositions (supermodes). Bright squeezing in multiple discrete frequency modes, if realized, could unlock novel applications in quantum-enhanced spectroscopy and optical quantum computing. Here, we show how $Q$ factor engineering of a multimode nonlinear cavity with cascaded three wave mixing processes creates strong, spectrally tunable single mode output amplitude noise squeezing over 10 dB below the shot noise limit. In addition, we demonstrate squeezing for multiple discrete frequency modes above threshold. This bright squeezing arises from enhancement of the (noiseless) nonlinear rate relative to decay rates in the system due to the cascaded generation of photons in a single idler "bath" mode. A natural consequence of the strong nonlinear coupling in our system is the creation of an effective cavity in the synthetic frequency dimension that sustains Bloch oscillations in the modal energy distribution. Bloch mode engineering could provide an opportunity to better control nonlinear energy flow in the synthetic frequency dimension, with exciting applications in quantum random walks and topological photonics. Lastly, we show evidence of long-range correlations in amplitude noise between discrete frequency modes, pointing towards the potential of long-range entanglement in a synthetic frequency dimension.

#3 Unclocklike biological oscillators with frequency memory [PDF] [Copy] [Kimi]

Authors: Christian Mauffette Denis ; Paul François

Entrainment experiments on the vertebrate segmentation clock have revealed that embryonic oscillators actively change their internal frequency to adapt to the driving signal. This is neither consistent with a one-dimensional clock model nor with a limit-cycle model, but rather suggests a new "unclocklike" behavior. In this work, we propose simple biologically realistic descriptions of such internal frequency adaptation, where a phase oscillator activates a memory variable controlling the oscillator's frequency. We study two opposite limits for the control of the memory variable, one with a smooth phase-averaging memory field, and the other with a pulsatile, phase-dependent activation. Both models recapitulate intriguing properties of the entrained segmentation clock, such as very broad Arnold tongues and an entrainment phase plateauing with detuning. We compute analytically multiple properties of such systems, such as the entrainment phases and cycle shapes. We further describe new phenomena, including hysteresis in entrainment, bistability in the frequency of the entrained oscillator, and probabilistic entrainment. Our work shows oscillators with frequency memory can exhibit new classes of unclocklike properties, that can be tested experimentally.

#4 Network mutual information measures for graph similarity [PDF] [Copy] [Kimi]

Authors: Helcio Felippe ; Federico Battiston ; Alec Kirkley

A wide range of tasks in exploratory network analysis and machine learning, such as clustering network populations or identifying anomalies in temporal graph streams, require a measure of the similarity between two graphs. To provide a meaningful data summary for downstream scientific analyses, the graph similarity measures used in these unsupervised settings must be principled, interpretable, and capable of distinguishing meaningful overlapping network structure from statistical noise at different scales of interest. Here we derive a family of graph mutual information measures that satisfy these criteria and are constructed using only fundamental information theoretic principles. Our measures capture the information shared among networks according to different encodings of their structural information, with our mesoscale mutual information measure allowing for network comparison under any specified network coarse-graining. We test our measures in a range of applications on real and synthetic network data, finding that they effectively highlight intuitive aspects of network similarity across scales in a variety of systems.

#5 Stable electro-optic modulators using thin-film lithium tantalate [PDF] [Copy] [Kimi]

Authors: Keith Powell ; Xudong Li ; Daniel Assumpcao ; Letícia Magalhães ; Neil Sinclair ; Marko Lončar

We demonstrate electro-optic modulators realized in low-loss thin-film lithium tantalate with superior DC-stability (<1 dB power fluctuation from quadrature with 12.1 dBm input) compared to equivalent thin-film lithium niobate modulators (5 dB fluctuation) over 46 hours.

#6 Data-Error Scaling in Machine Learning on Natural Discrete Combinatorial Mutation-prone Sets: Case Studies on Peptides and Small Molecules [PDF] [Copy] [Kimi]

Authors: Vanni Doffini ; O. Anatole von Lilienfeld ; Michael A. Nash

We investigate trends in the data-error scaling behavior of machine learning (ML) models trained on discrete combinatorial spaces that are prone-to-mutation, such as proteins or organic small molecules. We trained and evaluated kernel ridge regression machines using variable amounts of computationally generated training data. Our synthetic datasets comprise i) two na\"ive functions based on many-body theory; ii) binding energy estimates between a protein and a mutagenised peptide; and iii) solvation energies of two 6-heavy atom structural graphs. In contrast to typical data-error scaling, our results showed discontinuous monotonic phase transitions during learning, observed as rapid drops in the test error at particular thresholds of training data. We observed two learning regimes, which we call saturated and asymptotic decay, and found that they are conditioned by the level of complexity (i.e. number of mutations) enclosed in the training set. We show that during training on this class of problems, the predictions were clustered by the ML models employed in the calibration plots. Furthermore, we present an alternative strategy to normalize learning curves (LCs) and the concept of mutant based shuffling. This work has implications for machine learning on mutagenisable discrete spaces such as chemical properties or protein phenotype prediction, and improves basic understanding of concepts in statistical learning theory.

#7 Principal Component Analysis for Spatial Phase Reconstruction in Atom Interferometry [PDF] [Copy] [Kimi]

Authors: Stefan Seckmeyer ; Holger Ahlers ; Jan-Niclas Kirsten-Siemß ; Matthias Gersemann ; Ernst M. Rasel ; Sven Abend ; Naceur Gaaloul

Atom interferometers are sensitive to a wide range of forces by encoding their signals in interference patterns of matter waves. To estimate the magnitude of these forces, the underlying phase shifts they imprint on the atoms must be extracted. Up until now, extraction algorithms typically rely on a fixed model of the patterns' spatial structure, which if inaccurate can lead to systematic errors caused by, for example, wavefront aberrations of the used lasers. In this paper we employ an algorithm based on Principal Component Analysis, which is capable of characterizing the spatial phase structure and per image phase offsets of an atom interferometer from a set of images. The algorithm does so without any prior knowledge about the specific spatial pattern as long as this pattern is the same for all images in the set. On simulated images with atom projection noise we show the algorithm's reconstruction performance follows distinct scaling laws, i.e., it is inversely-proportional to the square-root of the number atoms or the number of images respectively, which allows a projection of its performance for experiments. We also successfully extract the spatial phase patterns of two experimental data sets from an atom gravimeter. This algorithm is a first step towards a better understanding and complex spatial phase patterns, e.g., caused by inhomogeneous laser fields in atom interferometry.

#8 Knowledge Gaps and Research Needs for Modeling CO2 Mineralization in the Basalt-CO2-Water System: A Review of Laboratory Experiments [PDF] [Copy] [Kimi]

Authors: Peng Lu ; John Apps ; Guanru Zhang ; Alexander Gysi Chen Zhu

Carbon capture and storage in basalt is being actively investigated as a scalable climate change mitigation option. Accurate geochemical modeling prediction of the extent and rate of CO2 mineralization is a critical component in assessing the local and global feasibility and efficacy of this strategy. In this study, we review basalt-CO2-water interaction experimental studies conducted during the last two decades to determine whether they provide useable information for geochemical modeling. Most of the cited experiments generate data on the temporal evolution of water composition, and a few provide identification of secondary precipitates and their compositions, offering empirical and semi-quantitative information about the reactivity of basalts and the likelihood of secondary carbonate mineralization at various temperatures, pHs, and pCO2 conditions. However, most experiments provide insufficient information on the properties and quantity of secondary minerals formed, prohibiting accurate mass balance calculations and hence more quantitative geochemical modeling studies. Primary Ca, Mg, and Fe-bearing minerals in basalt control the availability of major ions released into aqueous solution for carbonate precipitation, and many secondary minerals, i.e., smectites, Ca-Mg-Fe carbonates, and zeolites, provide sinks for the same major ions, some of which are difficult to quantify experimentally. Thus, we have a multi-source and multi-sink inverse mass balance problem with insufficient constraints on the bulk system in which the temporal evolution of major ions does not provide sufficient information on which mineral(s) dissolve or the sequence of dissolution and precipitation reactions. Going forward, we propose that future experimental work should focus on trace elements and multiple isotopic tracers and better characterize the solid reaction products with modern analytical instruments.

#9 Universality of opinions disappearing in sociophysical models of opinion dynamics: From initial multitude of opinions to ultimate consensus [PDF] [Copy] [Kimi]

Authors: Maciej Wołoszyn ; Tomasz Masłyk ; Szymon Pająk ; Krzysztof Malarz

Possibility of reaching a consensus in social systems with strong initial fragmentation is one of the most interesting issues in sociopysics. It is also intriguing what the dynamics of such processes is. To address those problems, we performed computer simulations using well-established models of social opinion formation, namely the voter, Sznajd, and Latan\'e models. We investigated opinion dynamics in cases where the initial number of opinions is very large, equal to the number of actors (the voter and Latan\'e models) or when every second actor has their own opinion (Sznajd model), with some variations on the update schemes, lattice topologies, effective ranges of interaction, and information noise levels. For all considered models, the number of opinions assumed by the actors is finally almost always reduced to only one. However, while the voter and Latan\'e models exhibit a power-law time decrease in the number of opinions, the Sznajd model follows a complex three-stage behavior. We also demonstrated that the mean/median time of reaching the consensus scales with system size according to a power law for voter and Sznajd models, while for the Latan\'e model this increase is even faster. Our results show that in the studied models the consensus is possible, provided that a long enough and model-dependent time to reach this state is available.

#10 Understanding solid nitrogen through machine learning simulation [PDF] [Copy] [Kimi]

Authors: Marcin Kirsz ; Ciprian G. Pruteanu ; Graeme J. Ackland

We construct a fast, transferable, general purpose, machine-learning interatomic potential suitable for large-scale simulations of $N_2$. The potential is trained only on high quality quantum chemical molecule-molecule interactions, no condensed phase information is used. The potential reproduces the experimental phase diagram including the melt curve and the molecular solid phases of nitrogen up to 10 GPa. This demonstrates that many-molecule interactions are unnecessary to explain the condensed phases of $N_2$. With increased pressure, transitions are observed from cubic ($\alpha-N_2$), which optimises quadrupole-quadrupole interactions, through tetragonal ($\gamma-N_2$) which allows more efficient packing, through to monoclinic ($\lambda-N_2$) which packs still more efficiently. On heating, we obtain the hcp 3D rotor phase ($\beta-N_2$) and, at pressure, the cubic $\delta-N_2$ phase which contains both 3D and 2D rotors, tetragonal $\delta^\star-N_2$ phase with 2D rotors and the rhombohedral $\epsilon-N_2$. Molecular dynamics demonstrates where these phases are indeed rotors, rather than frustrated order. The model does not support the existence of the wide range of bondlengths reported for the complex $\iota-N_2$ phase. The thermodynamic transitions involve both shifts of molecular centres and rotations of molecules. We simulate these phase transitions between finding that the onset of rotation is rapid whereas motion of molecular centres is inhibited and the cause of the observed sluggishness of transitions. Routine density functional theory calculations give a similar picture to the potential.

#11 Transition frequencies, isotope shifts, and hyperfine structure in $4s \rightarrow 4p$ transitions of Ti$^+$ ions [PDF] [Copy] [Kimi]

Authors: Tim Ratajczyk ; Kristian König ; Philipp Bollinger ; Tim Lellinger ; Victor Varentsov ; Wilfried Nörtershäuser ; Julien Spahn

We have measured transition frequencies, isotope shifts and hyperfine structure splittings in the $3 d^{2}\left({ }^{3\!}F\right) 4 s\,{ }^{4} F_J\rightarrow 3 d^{2}\left({ }^{3\!} F\right) 4 p \,^{4} G_{J+1}$ transitions in Ti$^+$ ions for $J=\frac{3}{2},\, \frac{5}{2},\, \frac{7}{2}$ using collinear laser spectroscopy. Ions were generated by laser ablation in a buffer gas atmosphere and extracted into vacuum through a nozzle and a pair of radiofrequency (RF) funnels. The obtained results are of interest as reference values for on-line measurements of short-lived titanium isotopes and for astrophysical searches for temporal or spatial variations of the fine structure constant $\alpha$ using quasar absorption spectra.

#12 Beam-induced backgrounds measured in the ATLAS detector during local gas injection into the LHC beam vacuum [PDF] [Copy] [Kimi]

Author: ATLAS Collaboration

Inelastic beam-gas collisions at the Large Hadron Collider (LHC), within a few hundred metres of the ATLAS experiment, are known to give the dominant contribution to beam backgrounds. These are monitored by ATLAS with a dedicated Beam Conditions Monitor (BCM) and with the rate of fake jets in the calorimeters. These two methods are complementary since the BCM probes backgrounds just around the beam pipe while fake jets are observed at radii of up to several metres. In order to quantify the correlation between the residual gas density in the LHC beam vacuum and the experimental backgrounds recorded by ATLAS, several dedicated tests were performed during LHC Run 2. Local pressure bumps, with a gas density several orders of magnitude higher than during normal operation, were introduced at different locations. The changes of beam-related backgrounds, seen in ATLAS, are correlated with the local pressure variation. In addition the rates of beam-gas events are estimated from the pressure measurements and pressure bump profiles obtained from calculations. Using these rates, the efficiency of the ATLAS beam background monitors to detect beam-gas events is derived as a function of distance from the interaction point. These efficiencies and characteristic distributions of fake jets from the beam backgrounds are found to be in good agreement with results of beam-gas simulations performed with the FLUKA Monte Carlo programme.

#13 Photon Squeezing in Photonic Time Crystals [PDF] [Copy] [Kimi]

Authors: Jaime Echave-Sustaeta ; Francisco J. García-Vidal ; P. A. Huidobro

Time-varying media offer a platform to realize novel and exotic wave effects, including photonic time crystals characterized by momentum band gaps with exponential wave amplification. Here we focus on the quantum electrodynamical properties of time-varying media, in particular vacuum amplification and squeezing. For that purpose, we present a theory of photon pair generation in photonic time crystals that unveils the link between the classical and quantum electrodynamical properties of these systems, that is, a direct relation between reflectivity and pair generation through the squeezing parameter. By working within an Hermitian framework, we are able to characterize quantum pair generation processes in photonic time crystals, showing how momentum bandgaps result in a non-resonant exponential enhancement of dynamical Casimir processes.

#14 Range separation of the interaction potential in intermolecular and intramolecular symmetry-adapted perturbation theory [PDF] [Copy] [Kimi]

Authors: Du Luu ; Clemence Corminboeuf ; Konrad Patkowski

Symmetry-adapted perturbation theory (SAPT) is a popular and versatile tool to compute and decompose noncovalent interaction energies between molecules. The intramolecular SAPT (ISAPT) variant provides a similar energy decomposition between two nonbonded fragments of the same molecule, covalently connected by a third fragment. In this work, we explore an alternative approach where the noncovalent interaction is singled out by a range separation of the Coulomb potential. We investigate two common splittings of the $1/r$ potential into long-range and short-range parts based on the Gaussian and error functions, and approximate either the entire intermolecular/interfragment interaction or only its attractive terms by the long-range contribution. These range separation schemes are tested for a number of intermolecular and intramolecular complexes. We find that the energy corrections from range-separated SAPT or ISAPT are in reasonable agreement with complete SAPT/ISAPT data. This result should be contrasted with the inability of the long-range multipole expansion to describe crucial short-range charge penetration and exchange effects; it shows that the long-range interaction potential does not just recover the asymptotic interaction energy but also provides a useful account of short-range terms. The best consistency is attained for the error-function separation applied to all interaction terms, both attractive and repulsive. This study is the first step towards a fragmentation-free decomposition of intramolecular nonbonded energy.

#15 Monitoring of neoadjuvant chemotherapy through time domain diffuse optics: Breast tissue composition changes and collagen discriminative potential [PDF] [Copy] [Kimi]

Authors: Nikhitha Mule ; Giulia Maffeis ; Rinaldo Cubeddu ; Carolina Santangelo ; Giampaolo Bianchini ; Pietro Panizza ; Paola Taroni

The purpose of this clinical study is to test broad spectral range (635-1060 nm) time domain diffuse optical spectroscopy in monitoring the response of breast cancer patients to neoadjuvant chemotherapy (NAC). The broadband operation allows us to fully analyze tissue composition in terms of hemoglobin, water, lipids and collagen concentration, which has never been systematically studied until now during the course of therapy. Patients are subjected to multiple breast optical imaging sessions, each one performed at different stages of NAC, both on tumor-bearing and contralateral healthy breast. We correlate the optical results with conventional imaging techniques and pathological response. Preliminary outcomes on data of 10 patients show an average significant reduction in the concentrations of oxy-hemoglobin (-53%, p = 0.0020), collagen (-36%, p = 0.0039) and water (-15%, p = 0.0195), and increase in lipids (+39%, p = 0.0137) from baseline to the end of therapy in the tumor-bearing breast of patients who responded to therapy at least partially. With respect to scattering, the scattering amplitude, a, increases slightly by the end of the therapy compared to the baseline, while the scattering slope, b, shows no significant change. Some change in concentrations of constituents was noticed also in the contralateral healthy breast, even though it was significant only for oxy-hemoglobin concentration. Collagen seems to be the only component distinguishing between complete (p = 0.0423) and partial responders (p = 0.1012) by the end of 2-3 weeks after the start of therapy, suggesting that it could serve as a potential biomarker to measure NAC effectiveness early during treatment. These initial results are consistent with those of standard medical modalities and highlight the sensitivity of the technique to changes that occur in breast composition during NAC.

#16 An anti-noise seismic inversion method based on diffusion model [PDF] [Copy] [Kimi]

Authors: Yingtian Liu ; Yong Li ; Xingan Hao ; Huating Li ; Zhangquan Liao ; Junheng Peng

Seismic impedance inversion is one of the most important part of geophysical exploration. However, due to random noise, the traditional semi-supervised learning (SSL) methods lack generalization and stability. To solve this problem, some authors have proposed SSL methods with anti-noise function to improve noise robustness and inversion accuracy. However, such methods are often not ideal when faced with strong noise. In addition, Low-frequency impedance models can mitigate this problem, but creating accurate low-frequency models is difficult and error-prone when well-log data is sparse and subsurface structures are complex. To address those issues, we propose a novel deep learning inversion method called DSIM-USSL (Unsupervised and Semi-supervised joint Learning for Seismic Inversion based on diffusion model). Specifically, we are the first to introduce a diffusion model with strong noise tendency and construct a diffusion seismic inversion model (DSIM). In the reverse diffusion of DSIM, we design the encoder-decoder which combines CNN for capturing local features and GRU for global sequence modeling; and we choose U-net to learn the distribution of random noise, enhancing the generalization and stability of proposed method. Furthermore, to further improve generalization of the proposed method, a two-step training approach (USSL) is utilized. First, an unsupervised trained encoder-decoder is used as the initial network model in place of the traditional low-frequency wave impedance model that is difficult to accurately acquire. Then, the SSL is employed to further optimize the encoder-decoder model. Experimental results on the Marmousi2 model and field data demonstrate that the DSIM-USSL method achieves higher accuracy in the presence of seismic data with random noise, and maintains high stability even under strong noise conditions.

#17 Kerr nonlinearity effect on the stability of Wannier-Stark states in active optical systems [PDF] [Copy] [Kimi]

Authors: Alexey Verbitskiy ; Alexey Yulin

The paper provides an analytical and numerical investigation of the dynamics of a one-dimensional chain of coupled optical resonators with conservative cubic nonlinearity and the gain saturated by nonlinear losses. The linear dependency of the resonator eigenfrequencies on their indexes makes it possible to use Wannier-Stark states as lasing modes. Numerical simulations have shown that the dependency of the resonant frequencies on the light intensity strongly affects the stability of Wannier-Stark states. To explain the observed destabilization of monochromatic lasing based on Wannier-Stark states a simple perturbation theory has been developed and compared with the data obtained in the numerical simulations.

#18 Tunable Localisation in Parity-Time-Symmetric Resonator Arrays with Imaginary Gauge Potentials [PDF] [Copy] [Kimi]

Authors: Habib Ammari ; Silvio Barandun ; Ping Liu ; Alexander Uhlmann

The aim of this paper is to illustrate both analytically and numerically the interplay of two fundamentally distinct non-Hermitian mechanisms in a deep subwavelength regime. Considering a parity-time symmetric system of one-dimensional subwavelength resonators equipped with two kinds of non-Hermiticity - an imaginary gauge potential and on-site gain and loss - we prove that all but two eigenmodes of the system decouple when going through an exceptional point. By tuning the gain-to-loss ratio, the system changes from a phase with unbroken parity-time symmetry to a phase with broken parity-time symmetry. At the macroscopic level, this is observed as a transition from symmetrical eigenmodes to condensated eigenmodes at one edge of the structure. Mathematically, it arises from a topological state change. The results of this paper open the door to the justification of a variety of phenomena arising from the interplay between non-Hermitian reciprocal and non-reciprocal mechanisms not only in subwavelength wave physics but also in quantum mechanics where the tight binding model coupled with the nearest neighbour approximation can be analysed with the same tools as those developed here.

#19 Optical characterization of size- and substrate-dependent performance of ultraviolet hybrid plasmonic nanowire lasers [PDF] [Copy] [Kimi]

Authors: Francesco Vitale ; Stephen A. Church ; Daniel Repp ; Karthika S. Sunil ; Mario Ziegler ; Marco Diegel ; Andrea Dellith ; Thi-Hien Do ; Sheng-Di Lin ; Jer-Shing Huang ; Thomas Pertsch ; Patrick Parkinson ; Carsten Ronning

Nanowire-based plasmonic lasers are now established as nano-sources of coherent radiation, appearing as suitable candidates for integration into next-generation nanophotonic circuitry. However, compared to their photonic counterparts, their relatively high losses and large lasing thresholds still pose a burdening constraint on their scalability. In this study, the lasing characteristics of ZnO nanowires on Ag and Al substrates, operating as optically-pumped short-wavelength plasmonic nanolasers, are systematically investigated in combination with the size-dependent performance of the hybrid cavity. A hybrid nanomanipulation-assisted single nanowire optical characterization combined with high-throughput PL spectroscopy enables the correlation of the lasing characteristics to the metal substrate and the nanowire diameter. The results evidence that the coupling between excitons and surface plasmons is closely tied to the relationship between substrate dispersive behavior and nanowire diameter. Such coupling dictates the degree to which the lasing character, be it more plasmonic- or photonic-like, can define the stimulated emission features and, as a result, the device performance.

#20 Improved calculation of the second-order ocean Doppler spectrum for sea state inversion [PDF] [Copy] [Kimi]

Author: Charles-Antoine Guérin

We propose a simple method, based on an original change of variables, for the fast and accurate calculation of the second-order ocean Doppler spectrum that describes the sea echo of High-Frequency radars. A byproduct of the technique is the derivation of an improved weighting function which can be used for the inversion of the main sea state parameters. For this we revisit Barrick's method for the estimation of the significant wave height and the mean period from the ocean Doppler spectrum. On the basis of numerical simulations we show that a better estimation of these parameters can be reached but necessitates a preliminary bias correction that depends only on the radar frequency. A second consequence of our formulation is the derivation of a simple yet analytical nonlinear approximation of the second-order ocean Doppler spectrum when the Doppler frequency is larger than the Bragg frequency. This opens new perspective for the inversion of directional wave spectra from High-Frequency radar measurements.

#21 Application of SPOD to the analysis of sound field of aeroacoustic sources [PDF] [Copy] [Kimi]

Authors: Oleg Bychkov ; Georgy Faranosov

Spectral proper orthogonal decomposition (SPOD) is proposed for identification of the multipole structure of aeroacoustic noise sources from far-field measurements. The method is verified via tests with point multipoles and validated using experimental data on flow-induced cylinder noise and turbulent jet noise.

#22 Accurate estimation of the normalized mutual information of multidimensional data [PDF] [Copy] [Kimi]

Authors: Daniel Nagel ; Georg Diez ; Gerhard Stock

While the linear Pearson correlation coefficient represents a well-established normalized measure to quantify the interrelation of two stochastic variables $X$ and $Y$, it fails for multidimensional variables such as Cartesian coordinates. Avoiding any assumption about the underlying data, the mutual information $I(X, Y)$ does account for multidimensional correlations. However, unlike the normalized Pearson correlation, it has no upper bound ($I \in [0, \infty)$), i.e., it is not clear if say, $I = 0.4$ corresponds to a low or a high correlation. Moreover, the mutual information (MI) involves the estimation of high-dimensional probability densities (e.g., six-dimensional for Cartesian coordinates), which requires a k-nearest neighbor algorithm, such as the estimator by Kraskov et al. [Phys. Rev. E 69, 066138 (2004)]. As existing methods to normalize the MI cannot be used in connection with this estimator, a new approach is presented, which uses an entropy estimation method that is invariant under variable transformations. The algorithm is numerically efficient and does not require more effort than the calculation of the (un-normalized) MI. After validating the method by applying it to various toy models, the normalized MI between the $C_{\alpha}$ -coordinates of T4 lysozyme is considered and compared to a correlation analysis of inter-residue contacts.

#23 Trapping and extreme clustering of finitely-dense inertial particles near a rotating vortex pair [PDF] [Copy] [Kimi]

Authors: Saumav Kapoor ; Divya Jaganathan ; Rama Govindarajan

Small heavy particles cannot get attracted into a region of closed streamlines in a non-accelerating frame (Sapsis & Haller 2010). In a rotating system, however, particles can get trapped (Angilella 2010) near vortices. We perform numerical simulations examining trapping of inertial particles in a prototypical rotating flow: an identical pair of rotating Lamb-Oseen vortices, without gravity. Our parameter space includes the particle Stokes number $St$, measuring the particle's inertia, and a density parameter $R$, measuring the particle-to-fluid relative density. We focus on inertial particles that are finitely denser than the fluid. Particles can get indefinitely trapped near the vortices and display extreme clustering into smaller dimensional objects: attracting fixed-points, limit cycles and chaotic attractors. As $St$ increases for a given $R$, we may have an incomplete or complete period-doubling route to chaos, as well as an unusual period-halving route back to a fixed-point attractor. The fraction of trapped particles can vary non-monotonically with $St$. We may even have windows in $St$ for which no particle trapping occurs. At $St$ larger than a critical value, beyond no trapping occurs, significant fractions of particles can spend long but finite times in the vortex vicinity. The inclusion of the Basset-Boussinesq history (BBH) force is imperative in our study due to particle's finite density. BBH force significantly increases the basin of attraction as well as the range of $St$ where trapping can occur. Extreme clustering can be physically significant in planetesimal formation by dust aggregation in protoplanetary disks, phytoplankton aggregation in oceans, etc.

#24 Predicting the binding of small molecules to proteins through invariant representation of the molecular structure [PDF2] [Copy] [Kimi1]

Authors: R. Beccaria ; A. Lazzeri ; G. Tiana

We present a computational scheme for predicting the ligands that bind to a pocket of known structure. It is based on the generation of a general abstract representation of the molecules, which is invariant to rotations, translations and permutations of atoms, and has some degree of isometry with the space of conformations. We use these representations to train a non-deep machine learning algorithm to classify the binding between pockets and molecule pairs, and show that this approach has a better generalization capability than existing methods.

#25 Urban Boundary Delineation from Commuting Data with Bayesian Stochastic Blockmodeling: Scale, Contiguity, and Hierarchy [PDF] [Copy] [Kimi]

Authors: Sebastian Morel-Balbi ; Alec Kirkley

A common method for delineating urban and suburban boundaries is to identify clusters of spatial units that are highly interconnected in a network of commuting flows, each cluster signaling a cohesive economic submarket. It is critical that the clustering methods employed for this task are principled and free of unnecessary tunable parameters to avoid unwanted inductive biases while remaining scalable for high resolution mobility networks. Here we systematically assess the benefits and limitations of a wide array of Stochastic Block Models (SBMs)$\unicode{x2014}$a family of principled, nonparametric models for identifying clusters in networks$\unicode{x2014}$for delineating urban spatial boundaries with commuting data. We find that the data compression capability and relative performance of different SBM variants heavily depends on the spatial extent of the commuting network, its aggregation scale, and the method used for weighting network edges. We also construct a new measure to assess the degree to which community detection algorithms find spatially contiguous partitions, finding that traditional SBMs may produce substantial spatial discontiguities that make them challenging to use in general for urban boundary delineation. We propose a fast nonparametric regionalization algorithm that can alleviate this issue, achieving data compression close to that of unconstrained SBM models while ensuring spatial contiguity, benefiting from a deterministic optimization procedure, and being generalizable to a wide range of community detection objective functions.