Date: Fri, 9 Aug 2024 | Total: 58

Molecular representation is a foundational element in our understanding of the physical world. Its importance ranges from the fundamentals of chemical reactions to the design of new therapies and materials. Previous molecular machine learning models have employed strings, fingerprints, global features, and simple molecular graphs that are inherently information-sparse representations. However, as the complexity of prediction tasks increases, the molecular representation needs to encode higher fidelity information. This work introduces a novel approach to infusing quantum-chemical-rich information into molecular graphs via stereoelectronic effects. We show that the explicit addition of stereoelectronic interactions significantly improves the performance of molecular machine learning models. Furthermore, stereoelectronics-infused representations can be learned and deployed with a tailored double graph neural network workflow, enabling its application to any downstream molecular machine learning task. Finally, we show that the learned representations allow for facile stereoelectronic evaluation of previously intractable systems, such as entire proteins, opening new avenues of molecular design.

While intuitive for humans, the concept of visual complexity is hard to define and quantify formally. We suggest adopting the multi-scale structural complexity (MSSC) measure, an approach that defines structural complexity of an object as the amount of dissimilarities between distinct scales in its hierarchical organization. In this work, we apply MSSC to the case of visual stimuli, using an open dataset of images with subjective complexity scores obtained from human participants (SAVOIAS). We demonstrate that MSSC correlates with subjective complexity on par with other computational complexity measures, while being more intuitive by definition, consistent across categories of images, and easier to compute. We discuss objective and subjective elements inherently present in human perception of complexity and the domains where the two are more likely to diverge. We show how the multi-scale nature of MSSC allows further investigation of complexity as it is perceived by humans.

Calcium carbonate plays a crucial role in the global carbon cycle, and its phase diagram has always been of significant scientific interest. In this study, we used molecular dynamics (MD) to investigate several structural phase transitions of calcium carbonate. Using the Raiteri potential model, we explored the structural transitions occurring at a constant pressure of 1 bar with temperatures ranging from 300 K to 2500 K, and at a constant temperature of 1600 K with pressures ranging from 0 to 13 GPa. At increasing temperatures, the transitions calcite, CaCO$_3$-IV, and CaCO$_3$-V are observed and characterized. Within the calcite structure, CO$_3^{2-}$ ions are ordered between layers. As temperature increases, the calcite to CaCO$_3$-IV transition occurs, determining the partial disordering of CO$_3^{2-}$ ions. At a higher temperature, CaCO$_3$-IV transforms into CaCO$_3$-V. By applying free energy analysis, we have classified the last transition as a continuous order-disorder phase transition. At a temperature of 2000 K, it appears a `disordered CaCO$_3$' structure, characterized by low order within the calcium and carbonate sublattices and the free rotation of CO$_3^{2-}$ ions. At increasing pressures, two calcium carbonate transformations were observed. At $P=$ 2 GPa, the CaCO$_3$-IV phase undergoes a phase transition into CaCO$_3$-V, demonstrating that the model can describe the transition between these two phases as pressure and temperature-driven. Another phase transition was detected at $P=$ 4.25 GPa -- CaCO$_3$-V transits into the recently discovered CaCO$_3$-Vb phase. This transition is classified as a first-order phase transition by structural analysis and free energy-based arguments.

In this work, topological resonance behaviors of surface acoustic wave (SAW) under surface liquid-layer loadings are investigated. The results show that through the coupling of waveguide and resonant cavity, a resonant peak with high Q-factor can emerge, and the frequency of which is significantly sensitive to the liquid parameters. Based on that, a topological-resonance SAW liquid-phase sensor is proposed. The results show that the formation of the topological phase transition is closely related to the liquid-layer loading. Furthermore, by changing the thickness of the liquid-layer loading, the operating-frequency range of the device can be significantly expanded, which is vital for biomedical detecting applications. Based on which, the sensing performances of topological-resonance SAW liquid-phase sensor are simulated, which are used to sensing the concentration of albumin and hemoglobin concentration in blood, and high sensitivities and Q-factors can be obtained for this device. The results presented in this paper can provide an important basis for the realization of highly sensitive and stable SAW biomedical sensors in the future.

The widely used thermal Hartree-Fock (HF) theory is generalized to include the effect of electron correlation while maintaining its quasi-independent-particle framework. An electron-correlated internal energy (or grand potential) is defined by the second-order finite-temperature many-body perturbation theory (MBPT), which then dictates the corresponding thermal orbital (quasi-particle) energies in such a way that all thermodynamic relations are obeyed. The associated density matrix is of the one-electron type, whose diagonal elements take the form of the Fermi-Dirac distribution functions, when the grand potential is minimized. The formulas for the entropy and chemical potential are unchanged from those of Fermi-Dirac or thermal HF theory. The theory thus postulates a finite-temperature extension of the second-order Dyson self-energy of one-particle many-body Green's function theory and can be viewed as a second-order, diagonal, frequency-independent, thermal inverse Dyson equation. At low temperature, the theory approaches finite-temperature MBPT of the same order, but it outperforms the latter at intermediate temperature by including additional electron-correlation effects through orbital energies. A physical meaning of these thermal orbital energies (including that of thermal HF orbital energies, which has been elusive) is proposed.

In the past decade, it has been demonstrated that monolayers of metal dichalcogenides are well-suited for thermoelectric applications. ZrX2N4 (X = Si, Ge) is a reasonable choice for thermoelectric applications when considering a favorable value of the figure of merit in two-dimensional (2D) layered materials. In this study, we examined the thermoelectric characteristics of the two-dimensional monolayer of ZrX2N4 (where X can be either Si or Ge) using a combination of Density Functional Theory (DFT) and the Boltzmann Transport Equation (BTE). A thermoelectric figure of merit (ZT) of 0.90 was achieved at a temperature of 900 K for p-type ZrGe2N4, while a ZT of 0.83 was reported for n-type ZrGe2N4 at the same temperature. In addition, the ZrGe2N4 material exhibited a thermoelectric figure of merit (ZT) of around 0.7 at room temperature for the p-type. Conversely, the ZrSi2N4 exhibited a relatively lower thermoelectric figure of merit (ZT) at ambient temperature. At higher temperatures, the ZT value experiences a substantial increase, reaching 0.89 and 0.82 for p-type and n-type materials, respectively, at 900 K. Through our analysis of the electronic band structure, we have determined that ZrSi2N4 and ZrGe2N4 exhibit indirect bandgaps (BG) of 2.74 eV and 2.66 eV, respectively, as per the Heyd-Scuseria-Ernzerhof (HSE) approximation.

Crush curves are of fundamental importance to numerical modeling of small and porous astrophysical bodies. The empirical literature often measures them for silica grains, and different studies have used various methods, sizes, textures, and pressure conditions. Here we review past studies and supplement further experiments in order to develop a full and overarching understanding of the silica crush curve behavior. We suggest a new power-law function that can be used in impact simulations of analog materials similar to micro-granular silica. We perform a benchmarking study to compare this new crush curve to the parametric quadratic crush curve often used in other studies, based on the study case of the DART impact onto the asteroid Dimorphos. We find that the typical quadratic crush curve parameters do not closely follow the silica crushing experiments, and as a consequence they under (over) estimate compression close (far) from the impact site. The new crush curve presented here, applicable to pressures between a few hundred Pa and up to 1.1 GPa, might therefore be more precise. Additionally, it is not calibrated by case-specific parameters, and can be used universally for comet- or asteroid-like bodies, given an assumed composition similar to micro-granular silica.

Following the recent interest in van der Waals (vdW) ferroelectrics, topologically nontrivial polar structures have been predicted to form in twisted bilayers. Due to the unconventional nature of vdW ferroelectricity, these topological polar structures have proven difficult to observe experimentally. Here we propose that these textures may be probed optically by showing that topological polarization textures result in exotic nonlinear optical responses. We derive this relationship analytically using non-Abelian Berry connections and a quantum-geometric framework, supported by tight-binding and first-principles calculations. For the case of moir\'e materials without centrosymmetry, which form networks of polar merons and antimerons, the shift photoconductivity forms a vortex-like structure in real space. For a range of frequencies where transitions occur at the Brillouin zone edge, the shift photocurrents are anti-parallel to the in-plane electronic polarization field. Our findings highlight the interplay between complex polarization textures and nonlinear optical responses in vdW materials and provide a sought-after strategy for the experimental detection of topological polarization structures.

The R-value is a measure of the strength of photospheric magnetic Polarity Inversion Lines (PILs) in Active Regions (ARs). This work investigates the possibility of a relation between R-value variations and the occurrence of X-class flares in ARs, not in the solar photosphere, as usual, but above it in regions, closer to where flares occur. The modus operandi is to extrapolate the Solar Dynamic Observatory's (SDO) Helioseismic and Magnetic Imager (HMI) magnetogram data up to a height of 3.24 Mm above the photosphere and then compute the R-value based on the extrapolated magnetic field. Recent studies have shown that certain flare-predictive parameters such as the horizontal gradient of the vertical magnetic field and magnetic helicity may improve flare prediction lead times significantly if studied at a specific height range above the photosphere, called the Optimal Height Range (OHR). Here we define the OHR as a collection of heights where a sudden but sustained increase in R-value is found. For the eight case studies discussed in this paper, our results indicate that it is possible for OHRs to exist in the low solar atmosphere (between 0.36 - 3.24 Mm), where R-value spikes occur 48-68 hrs before the first X-class flare of an emerging AR. The temporal evolution of R-value before the first X-class flare for an emerging AR is also found to be distinct from that of non-flaring ARs. For X-class flares associated with non-emerging ARs, an OHR could not be found.

In this work, we investigate the implications of the differential Hebbian learning rule known as Input-Correlations (ICO) learning in the classification of synchronization in coupled nonlinear oscillator systems. We are investigating the parity-time symmetric coupled Duffing oscillator system with nonlinear dissipation/amplification. In our investigation of the temporal dynamics of this system, it is observed that the system exhibits chaotic as well as quasiperiodic dynamics. On further investigation, it is found that the chaotic dynamics is distorted anti-phase synchronized, whereas the quasiperiodic dynamics is desynchronized. So, on the application of the ICO learning in these two parametric regimes, we observe that the weight associated with the stimulus remains constant when the oscillators are anti-phase synchronized, in spite of there being distortion in the synchronization. But when the oscillators exhibit quasiperiodic dynamics, there is erratic evolution of the weight with time. So, from this, it could be ascertained that the ICO learning could be made use of in the classification of synchronization dynamics in nonlinear systems.

The total energy of a fireball is commonly obtained from optical measurements with an assumed value for luminous efficiency. Acoustic energy measurements offer an independent means of energy estimation. Here we combine optical and acoustic methods to validate the luminous efficiency model of Borovi\v{c}ka et al. (2020). Our goal is to compare these models with acoustic measurements of meteoroid energy deposition. Employing theoretical blast scaling laws following the approach of McFadden et al. (2021), we determine explosive yields for both fireball fragmentation events and cylindrical shocks for four different bright fireballs. We model fireballs using the MetSim software (Vida et al., 2023) and find that the Borovi\v{c}ka et al. (2020) model produces agreement better than a factor of two for our three chondritic fireball case studies. The major exception is an iron meteorite-producing fireball where the luminous efficiency is an order of magnitude higher than model predictions calibrated with stony fireballs. We suggest that large disparities between optical and acoustic energies could be a signature of iron fireballs and hence useful as a discriminant of that population.

We propose an approach for characterizing source-impact pathways, the interactions of a set of variables in space-time due to an external forcing, in climate models using in-situ analyses that circumvent computationally expensive read/write operations. This approach makes use of a lightweight open-source software library we developed known as CLDERA-Tools. We describe how CLDERA-Tools is linked with the U.S. Department of Energy's Energy Exascale Earth System Model (E3SM) in a minimally invasive way for in-situ extraction of quantities of interested and associated statistics. Subsequently, these quantities are used to represent source-impact pathways with time-dependent directed acyclic graphs (DAGs). The utility of CLDERA-Tools is demonstrated by using the data it extracts in-situ to compute a spatially resolved DAG from an idealized configuration of the atmosphere with a parameterized representation of a volcanic eruption known as HSW-V.

Halide double perovskites are a chemically-diverse and growing class of compound semiconductors that are promising for optoelectronic applications. However, the prediction of their fundamental gaps and optical properties with density functional theory (DFT) and {\it ab initio} many-body perturbation theory has been a significant challenge. Recently, a nonempirical Wannier-localized optimally-tuned screened range-separated hybrid (WOT-SRSH) functional has been shown to accurately produce the fundamental band gaps of a wide set of semiconductors and insulators, including lead halide perovskites. Here we apply the WOT-SRSH functional to five halide double perovskites, and compare the results with those obtained from other known functionals and previous $GW$ calculations. We also use the approach as a starting point for $GW$ calculations and we compute the band structures and optical absorption spectrum for Cs\textsubscript{2}Ag{Bi}Br\textsubscript{6}, using both time-dependent DFT and the $GW$-Bethe-Salpeter equation approach. We show that the WOT-SRSH functional leads to accurate fundamental and optical band gaps, as well as optical absorption spectra, consistent with spectroscopic measurements, thereby establishing WOT-SRSH as a viable method for the accurate prediction of optoelectronic properties of halide double perovskites.

Photodetectors based on two-dimensional (2D) atomically thin semiconductors suffer from low light absorption, limiting their potential for practical applications. In this work, we demonstrate a high-performance MoS2 phototransistors by integrating few-layer MoS2 on a PN junction formed in a silicon (Si) substrate. The photovoltage created in the PN junction under light illumination electrically gates the MoS2 channel, creating a strong photoresponse in MoS2. We present an analytical model for the photoresponse of our device and show that it is in good agreement with measured experimental photocurrent in MoS2 and photovoltage in the Si PN junction. This device structure separates light absorption and electrical response functions, which provides us an opportunity to design new types of photodetectors. For example, incorporating ferroelectric materials into the gate structure can produce a negative capacitance that boosts gate voltage, enabling low power, high sensitivity phototransistor; this, combined with separating light absorption and electrical functions, enables advanced high-performance photodetectors.

Magnetic fields are an energetically important component of star-formation galaxies, but it is often difficult to measure their properties from observations. One of the complexities stems from the fact that the magnetic fields, especially in spiral galaxies, have a two-scale nature: a large-scale field, coherent over ${\rm kpc}$ scales and a small-scale, random field with a scale of $\lesssim$ $100~{\rm pc}$. Moreover, it is known that the strength of small- and large-scale fields are comparable and this makes it even harder to find their imprints in radio polarisation observations such as the Faraday rotation measure, ${\rm RM}$, which is the integral over the path length of the product of the thermal electron density and the parallel component of the magnetic field to the line of sight. Here, we propose and demonstrate the use of second-order structure functions of ${\rm RM}$ computed with multiple higher-order stencils as a powerful analysis to separate the small- and large-scale magnetic field components. In particular, we provide new methods and calibrations to compute the scale and the strength of the large-scale magnetic field in the presence of small-scale magnetic fluctuations. We then apply the method to find the scale of large-scale magnetic fields in the nearby galaxies M51 and NGC 6946, using archival data and further discuss the need for computing the ${\rm RM}$ structure functions with higher-order stencils. With multiple modern radio polarisation observatories and eventually the Square Kilometre Array, ${\rm RM}$ observations will significantly improve in quantity and quality, and the higher-order stencil structure function techniques developed here can be used to extract information about multiscale magnetic fields in galaxies.

We theoretically propose a method for implementing the Hamiltonian incorporating Heisenberg and Dzyaloshinskii-Moriya (DM) interactions within Rydberg atoms arranged in a two-dimensional square lattice, utilizing Floquet engineering. In our scheme, we use both global and local operations of the spins. The global operations can be realized by applying the microwave and the local operations can be realized by the locally addressing lasers, which yields the ac-Stark shift. Since our engineered Hamiltonian contains bond-dependent DM interactions, we expect the emergence of quantum skyrmions in the ground state.

Recent cutting-edge experiments have provided $in\,situ$ structure characterization and measurements of the pressure ($P$), density ($\bar{\rho}$) and temperature ($T$) of shock compressed silicon in the 100 GPa range of pressures and up to $\sim$10,000K. We present first-principles calculations in this $P,T,\bar{\rho}$ regime to reveal a plethora of novel liquid-liquid phase transitions (LPTs), making the interpretation of these experiments very challenging. The short-ranged ionic structure of the fluid is preserved under collective adjustments of many distant atoms when subject to compression and heating, with surprisingly little change in electrical and thermal conductivities $\sigma$ and $\kappa$. We match the experimental X-Ray Thompson scattering and X-ray diffraction data theoretically, and provide pressure isotherms, ionization data and compressibilities that support the above picture of liquid silicon as a highly complex LPT-driven ``glassy'' metallic liquid. These novel results are relevant to materials research, studies of planetary interiors, high-energy-density physics, and in laser-fusion studies.

Operating on a surface with high noise rates and requiring susceptible, yet small-scale detectors, the Neutrinos-Angra detector is an excellent platform for technological development and expertise in new detection methods. This report details the primary features of the detector, the electronics involved, and preliminary physics results from the operational phase, particularly the ON-OFF analysis comparing the signals with the reactor in operation (ON) and during the maintenance shut-down (OFF), demonstrating the detector's capability to monitor reactor activity. Additionally, we will briefly discuss the prospects of using a cryogenic calorimeter to detect neutrinos via Coherent Elastic Neutrino-Nucleus Scattering (CEvNS), highlighting potential advancements in neutrino detection technology. Looking ahead, the project promises to play a crucial role in the integration of Latin American scientists and engineers into global scientific collaborations, significantly contributing to the LASF4RI and the broader HECAP strategic framework.

Beyond the Sun-Earth line, spacecraft equipped with various solar telescopes are intended to be deployed at several different vantage points in the heliosphere to carry out coordinated, multi-view observations of the Sun and its dynamic activities. In this context, we investigate solar visibility by imaging instruments onboard the spacecraft orbiting the Sun-Earth Lagrange points L1, L4 and L5, respectively. An optimal arrival time for vertical periodic orbits stationed at L4 and L5 is determined based on geometric considerations that ensure maximum visibility of solar poles or higher latitudes per year. For a different set of orbits around the three Lagrange points (L1, L4 and L5), we calculate the visibility of the solar surface (i.e., observation days per year) as a function of the solar latitude. We also analyze where the solar limb viewed from one of the three Sun-Earth Lagrange points under consideration is projected onto the solar surface visible to the other two. This analysis particularly aims at determining the feasibility of studying solar eruptions, such as flares and coronal mass ejections, with coordinated observations of off-limb erupting coronal structures and their on-disk magnetic footpoints. In addition, visibility analysis of a feature (such as sunspots) on the solar surface is made for multiple spacecraft in various types of orbits with different inclinations to quantify the improvement in continuous tracking of the target feature for studying its long-term evolution from emergence, growth and to decay. A comprehensive comparison of observations from single (L1), double (L1 and L4) and multi-space missions (L1, L4 and L5) is carried out through our solar visibility analysis, and this may help us to design future space missions of constructing multiple solar observatories at the Sun-Earth Lagrange points.

We report the first experimental determination of the interspecies $s$-wave scattering length\,($a_{34}$) between the $2\,^3S_1\,(F=3/2,m_F=3/2)$ state of $^3$He$^*$ and the $2\,^3S_1\,(m_J=1)$ state of $^4$He$^*$. We determine $a_{34}$ by inducing oscillations in a trapped Bose-Einstein condensate of $^4$He$^*$ and measuring the damping rate of these oscillations due to the presence of $^3$He$^*$ atoms. The deduced value of $a_{34}=29\pm3$\,nm is in good agreement with theoretical predictions. The knowledge of this scattering length is important for many fundamental experiments between these helium isotopes.

We prove finite-time vorticity blowup in the compressible Euler equations in $\mathbb{R}^d$ for any $d \geq 3$, starting from smooth, localized, and non-vacuous initial data. This is achieved by lifting the vorticity blowup result from [CCSV24] in $\mathbb{R}^2$ to $\mathbb{R}^d$ and utilizing the axisymmetry in $\mathbb{R}^d$. At the time of the first singularity, both vorticity blowup and implosion occur on a sphere $S^{d-2}$. Additionally, the solution exhibits a non-radial implosion, accompanied by a stable swirl velocity that is sufficiently strong to initially dominate the non-radial components and to generate the vorticity blowup.

When we place conducting bodies in electrolyte solutions, their surface potential $\Phi_s$ appears to be much smaller in magnitude than the intrinsic one $\Phi_0$ and normally does not obey the classical electrostatic boundary condition of a constant surface potential expected for conductors. In this paper, we demonstrate that an explanation of these observations can be obtained by postulating that diffuse ions condense at the "wall" due to a reduced permittivity of a solvent. For small values of $\Phi_0$ the surface potential responds linearly. On increasing $\Phi_0$ further $\Phi_s$ augments nonlinearly and then saturates to a constant value. Analytical approximations for $\Phi_s$ derived for these three distinct modes show that it always adjusts to salt concentration, which is equivalent to a violation of the constant potential condition. The latter would be appropriate for highly dilute solutions, but only if $\Phi_0$ is small. Surprisingly, when the plateau with high $\Phi_s$ is reached, the conductor surface switches to a constant charge density condition normally expected for insulators. Our results are directly relevant for conducting electrodes, mercury drops, colloidal metallic particles and more.

Quasi two-dimensional (2D) colloidal synthesis made quantum confinement readily accessible in perovskites, generating additional momentum in perovskite LED research and lasing. Ultrathin perovskite layers exhibit high exciton binding energies and beneficial charge transport properties interesting for solar cells. In 2D perovskites, the combination of layers with different thickness helps to direct charge carriers in a targeted manner toward thicker layers with a smaller bandgap. However, detailed knowledge about the mechanisms by which excitons and charge carriers funnel and recombine in these structures is lacking. Here, we characterize colloidal 2D methylammonium lead bromide (MAPbBr3) Ruddlesden-Popper perovskites with a broad combination of layers (n = 3 to 10, and bulk fractions with n > 10) in one stack by femtosecond transient absorption spectroscopy and time-resolved photoluminescence, which gives comprehensive insights into the complexity of funneling and recombination processes. We find that after photoexcitation second- and third-order processes dominate in MAPbBr3 nanosheets, which indicates exciton-exciton annihilation (EEA) and Auger recombination. Long-lived excitons in thin layers (e.g., n = 5, Eb = 136 meV) funnel into high n with t = 10-50 ps, which decreases their exciton binding energy below kB T = 26 meV ( T = 300K) and leads to radiative recombination. Parallel and consecutive funneling compete with exciton trapping processes, making funneling an excellent tool to overcome exciton self-trapping when high-quality n-n interfaces are present. Free charge carriers in high n regions on the other hand facilitate radiative recombination and EEA is bypassed, which is desirable for LED and lasing applications.

We propose a scheme to perform optical pulses that suppress the effect of photon recoil by three orders of magnitude compared to ordinary pulses in the Lamb-Dicke regime. We derive analytical insight about the fundamental limits to the fidelity of optical qubits for trapped atoms and ions. This paves the way towards applications in quantum computing for realizing $>1000$ of gates with an overall fidelity above 99\%.

Using nonlinear simulations in two settings, we demonstrate that QG$^\mathrm{+1}$, a potential-vorticity (PV) based next-order-in-Rossby balanced model, captures several aspects of ocean submesoscale physics. In forced-dissipative 3D simulations under baroclinically unstable Eady-type background states, the statistical equilibrium turbulence exhibits long cyclonic tails and a plethora of rapidly-intensifying ageostrophic fronts. Despite that the model requires setting an explicit, small value for the fixed scaling Rossby number, the emergent flows are nevertheless characterized by $O(f)$ vorticity and convergence, as observed in upper-ocean submesoscale flows. Simulations of QG$^\mathrm{+1}$ under the classic strain-induced frontogenesis set-up show realistic frontal asymmetry and a finite time blow-up, quantitatively comparable to simulations of the semigeostrophic equations. The inversions in the QG$^\mathrm{+1}$ model are straightforward elliptic problems, allowing for the reconstruction of all flow fields from the PV and surface buoyancy, while avoiding the semigeostrophic coordinate transformation. Together, these results suggest QG$^\mathrm{+1}$ as a useful tool for studying upper-ocean submesoscale dynamics.