Search

The Markov chain Monte Carlo (MCMC) method is widely used in various fields as a powerful numerical integration technique for systems with many degrees of freedom. In MCMC methods, probabilistic state transitions can be considered as a random walk in state space, and random walks allow for sampling from complex distributions. However, paradoxically, it is necessary to carefully suppress the randomness of the random walk to improve computational efficiency. By breaking detailed balance, we can create a probability flow in the state space and perform more efficient sampling along this flow. Motivated by this idea, practical and efficient nonreversible MCMC methods have been developed over the past ten years. In particular, the lifting technique, which introduces probability flows in an extended state space, has been applied to various systems and has proven more efficient than conventional reversible updates. We review and discuss several practical approaches to implementing nonreversible MCMC methods, including the shift method in the cumulative distribution and the directed-worm algorithm. [ABSTRACT FROM AUTHOR]

Molecules under strong or ultra-strong light–matter coupling present an intriguing route to modify chemical structure, properties, and reactivity. A rigorous theoretical treatment of such systems requires handling matter and photon degrees of freedom on an equal quantum mechanical footing. In the regime of molecular electronic strong or ultra-strong coupling to one or a few molecules, it is desirable to treat the molecular electronic degrees of freedom using the tools of ab initio quantum chemistry, yielding an approach referred to as ab initio cavity quantum electrodynamics (ai-QED), where the photon degrees of freedom are treated at the level of cavity QED. We analyze two complementary approaches to ai-QED: (1) a parameterized ai-QED, a two-step approach where the matter degrees of freedom are computed using existing electronic structure theories, enabling the construction of rigorous ai-QED Hamiltonians in a basis of many-electron eigenstates, and (2) self-consistent ai-QED, a one-step approach where electronic structure methods are generalized to include coupling between electronic and photon degrees of freedom. Although these approaches are equivalent in their exact limits, we identify a disparity between the projection of the two-body dipole self-energy operator that appears in the parameterized approach and its exact counterpart in the self-consistent approach. We provide a theoretical argument that this disparity resolves only under the limit of a complete orbital basis and a complete many-electron basis for the projection. We present numerical results highlighting this disparity and its resolution in a particularly simple molecular system of helium hydride cation, where it is possible to approach these two complete basis limits simultaneously. In this same helium hydride system, we examine and compare the practical issue of the computational cost required to converge each approach toward the complete orbital and many-electron bases limit. Finally, we assess the aspect of photonic convergence for polar and charged species, finding comparable behavior between parameterized and self-consistent approaches. [ABSTRACT FROM AUTHOR]

Modeling the dynamics of a quantum system coupled to a dissipative environment becomes particularly challenging when the system's dimensionality is too high to permit the computation of its eigenstates. This problem is addressed by introducing an eigenstate-free formalism, where the open quantum system is represented as a mixture of high-dimensional, time-dependent wave packets governed by coupled Schrödinger equations, while the environment is described by a multi-component quantum master equation. An efficient computational implementation of this formalism is presented, employing a variational mixed Gaussian/multiconfigurational time-dependent Hartree (G-MCTDH) ansatz for the wave packets and propagating the environment dynamics via hierarchical equations, truncated at the first or second level of the hierarchy. The effectiveness of the proposed methodology is demonstrated on a 61-dimensional model of phonon-driven vibrational relaxation of an adsorbate. G-MCTDH calculations on 4- and 10-dimensional reduced models, combined with truncated hierarchical equations for the mean fields, nearly quantitatively replicate the full-dimensional quantum dynamical results on vibrational relaxation while significantly reducing the computational time. This approach thus offers a promising quantum dynamical method for modeling complex system–bath interactions, where a large number of degrees of freedom must be explicitly considered. [ABSTRACT FROM AUTHOR]

Electron attachment to pyridine results in electronic resonances, metastable states that can decay through electronic or nuclear degrees of freedom. This study uses orbital stabilization techniques combined with bound electronic structure methods, based on equation of motion coupled cluster or multi-reference methods, to calculate positions and widths of electronic resonances in pyridine that exist below 10 eV. We report four 2B1 and four 2A2 resonances, including one 2B1 not previously reported experimentally and two 2A2 resonances not reported at all in the literature. The two lower energy resonances are one-particle shape resonances, while the remaining are mixed or primarily core-excited resonances. Multi-reference perturbation theory provides the best description of these resonances, especially when their character is mixed. We describe the character of these resonances qualitatively and calculate Dyson orbitals, which provide information about their decay channels. [ABSTRACT FROM AUTHOR]

We present an ab initio method for computing vibro-polariton and phonon-polariton spectra of molecules and solids coupled to the photon modes of optical cavities. We demonstrate that if interactions of cavity photon modes with both nuclear and electronic degrees of freedom are treated on the level of the cavity Born–Oppenheimer approximation, spectra can be expressed in terms of the matter response to electric fields and nuclear displacements, which are readily available in standard density functional perturbation theory implementations. In this framework, results over a range of cavity parameters can be obtained without the need for additional electronic structure calculations, enabling efficient calculations on a wide range of parameters. Furthermore, this approach enables results to be more readily interpreted in terms of the more familiar cavity-independent molecular electric field response properties, such as polarizability and Born effective charges, which enter into the vibro-polariton calculation. Using corresponding electric field response properties of bulk insulating systems, we are also able to obtain the Γ point phonon-polariton spectra of two dimensional (2D) insulators. Results for a selection of cavity-coupled molecular and 2D crystal systems are presented to demonstrate the method. [ABSTRACT FROM AUTHOR]

Transition state theory (TST) provides a framework to estimate the rate of chemical reactions. Despite its great success with many reaction systems, the underlying assumptions such as local equilibrium and nonrecrossing do not necessarily hold in all cases. Although dynamical systems theory can provide the mathematical foundation of reaction tubes existing in phase space that enables us to predict the fate of reactions free from the assumptions of TST, numerical demonstrations for large systems have been yet one of the challenges. Here, we propose a dimensionality reduction algorithm to demonstrate structures in phase space (called reactive islands) that predict reactivity in systems with many degrees of freedom. The core of this method is the application of supervised principal component analysis, where a coordinate transformation is performed to preserve the dynamical information on reactivity (i.e., to which potential basin the system moves from a region of interest) as much as possible. The reactive island structures are expected to be reflected in the transformed, low-dimensional phase space. As an illustrative example, the algorithm is scrutinized using a modified Hénon–Heiles Hamiltonian system extended to many degrees of freedom, which has three channels leading to three different products from one stable potential basin. It is shown that our algorithm can predict the reactivity in the transformed, low-dimensional coordinate system better than a naïve coordinate system and that the reactivity distribution in the transformed low-dimensional space is considered to reflect the underlying reactive islands. [ABSTRACT FROM AUTHOR]

Owing to charge free property, magnon is highly promising to achieve dissipationless transport without Joule heating and, thus, potentially applicable to energy efficient devices. In this paper, using the non-equilibrium Green's function, we present the bulk-boundary correspondence for magnonic Kagome lattices by studying the edge magnons transport. With staggered exchange interaction and Dzyaloshinskii–Moriya interaction in the Kagome lattices, one can observe valley contrasting magnon Hall effect, which endows magnon transport with the valley degree of freedom and adds a new dimension to regulate magnon excitation. In particular, we demonstrate that the valley splitting in the Kagome lattice enables a tunable single edge chiral transport. Thermal rectification is a direction-dependent asymmetric heat transfer phenomenon; here, we report the tunable thermal rectification by asymmetric nonlinear effect, and it is, indeed, regulated by the Dzyaloshinskii–Moriya interaction direction (D → − D) and the exchange of J 1 and J 2 (J 1 ↔ J 2). Moreover, we show that the topological edge state mainly localizes around edges and leaks into the bulk with oscillatory decay. These give full play to spin and valley degrees of freedom and provide various avenues for information encoding and manipulation based on valley related magnonic flux. [ABSTRACT FROM AUTHOR]

In this work, we present a statistical adiabatic channel model for termolecular reactions, A + B + C → Products. Our approach relies on hyperspherical coordinates, where the adiabatic channels are readily defined in the hyper-radius after averaging the hyperangular degrees of freedom. In this way, we find a general expression for termolecular rate constants. We focus on ion-neutral association reactions to test our approach's accuracy and predictive power, finding a good agreement between theory and experiment, especially in those reactions' temperature dependence. [ABSTRACT FROM AUTHOR]

Chiral-induced-spin-selectivity of electron transport and its interplay with DNA's mechanical motion are explored in a double stranded DNA helix with spin–orbit-coupling. The mechanical degree of freedom is treated as a stochastic classical variable experiencing fluctuations and dissipation induced by the environment as well as force exerted by nonequilibrium, current-carrying electrons. Electronic degrees of freedom are described quantum mechanically using nonequilibrium Green's functions. Nonequilibrium Green's functions are computed along the trajectory for the classical variable taking into account dynamical, velocity dependent corrections. This mixed quantum-classical approach enables calculations of time-dependent spin-resolved currents. We showed that the electronic force may significantly modify the classical potential, which, at sufficient voltage, creates a bistable potential with a considerable effect on electronic transport. The DNA's mechanical motion has a profound effect on spin transport; it results in chiral-induced spin selectivity, increasing spin polarization of the current by 9% and also resulting in temperature-dependent current voltage characteristics. We demonstrate that the current noise measurement provides an accessible experimental means to monitor the emergence of mechanical instability in DNA motion. The spin resolved current noise also provides important dynamical information about the interplay between vibrational and spin degrees of freedom in DNA. [ABSTRACT FROM AUTHOR]

We describe how to simulate charge diffusion in organic semiconductors using a recently introduced mixed quantum–classical method, the mapping approach to surface hopping. In contrast to standard fewest-switches surface hopping, this method propagates the classical degrees of freedom deterministically on the most populated adiabatic electronic state. This correctly preserves the equilibrium distribution of a quantum charge coupled to classical phonons, allowing one to time-average along trajectories to improve the statistical convergence of the calculation. We illustrate the method with an application to a standard model for the charge transport in the direction of maximum mobility in crystalline rubrene. Because of its consistency with the equilibrium distribution, the present method gives a time-dependent diffusion coefficient that plateaus correctly to a long-time limiting value. The resulting mobility is somewhat higher than that of the relaxation time approximation, which uses a phenomenological relaxation parameter to obtain a non-zero diffusion coefficient from a calculation with static phonon disorder. However, it is very similar to the mobility obtained from Ehrenfest dynamics, at least in the parameter regimes we have investigated here. This is somewhat surprising because Ehrenfest dynamics overheats the electronic subsystem and is, therefore, inconsistent with the equilibrium distribution. [ABSTRACT FROM AUTHOR]

Understanding the dynamics of open quantum systems in strong coupling and non-Markovian regimes remains a formidable theoretical challenge. One popular and well-established method of approximation in these circumstances is provided by the polaron master equation (PME). In this work, we re-evaluate and extend the validity of the PME to capture the impact of non-Markovian polaron dressing, induced by non-equilibrium open system dynamics. By comparing with numerically exact techniques, we confirm that while the standard PME successfully predicts the dynamics of system observables that commute with the polaron transformation (e.g., populations in the Pauli z-basis), it can struggle to fully capture those that do not (e.g., coherences). This limitation stems from the mixing of system and environment degrees of freedom inherent to the polaron transformation, which affects the accuracy of calculated expectation values within the polaron frame. Employing the Nakajima–Zwanzig projection operator formalism, we introduce correction terms that provide an accurate description of observables that do not commute with the transformation. We demonstrate the significance of the correction terms in two cases, the canonical spin-boson model and a dissipative time-dependent Landau–Zener protocol, where they are shown to impact the system dynamics on both short and long timescales. [ABSTRACT FROM AUTHOR]

Although progress has been made in high-performance computing, there are still limitations on temporal and spatial scales of molecular dynamic calculations. A major issue in molecular dynamic (MD) simulations is the computational cost, and coarse-grained methods can save computational costs and accelerate calculations by reducing the degrees of freedom in the system. This method takes a selected group of representative atoms in the atomic microstructure as a bead and uses the proportional relationship of atomic scale atomic potentials in MD simulation to define the interaction between beads and solve the motion equation of beads. This article proposed a coarse-grained potential for graphite based on the modification of Airebo potential, with an n3: 1 mapping, and established a corresponding n3: 1 coarse-grained model, where one bead represents the number of n3 atoms. The results indicated that the coarse-grained model well reflected the basic thermal and mechanical properties of graphite. In addition, this article also proposed a coarsening method for the Lennard–Jones (L–J) potential function parameters and established a coarse-grained wetting model with n = 2. The results indicate that the coarse-grained wetting model can effectively predict the wetting performance of copper droplets on graphite in only 60% of the time using the all-atom model. The coarsening potential function and model proposed in this article are also applicable to graphite with rough surfaces. It can be anticipated that coarse-grained molecular dynamics methods will have more applications in the future, as they can handle large-scale calculations more quickly. [ABSTRACT FROM AUTHOR]

The large dimensionality of environments is the limiting factor in applying optimal control to open quantum systems beyond the Markovian approximation. Various methods exist to simulate non-Markovian systems, which effectively reduce the environment to a number of active degrees of freedom. Here, we show that several of these methods can be expressed in terms of a process tensor in the form of a matrix-product-operator, which serves as a unifying framework to show how they can be used in optimal control and to compare their performance. The matrix-product-operator form provides a general scheme for computing gradients using back propagation and allows the efficiency of the different methods to be compared via the bond dimensions of their respective process tensors. [ABSTRACT FROM AUTHOR]

Quantum states in complex aggregates are unavoidably affected by environmental effects, which typically cannot be accurately modeled by simple Markovian processes. As system sizes scale up, nonperturbative simulation becomes thus unavoidable, but they are extremely challenging due to the intimate interplay of intrinsic many-body interaction and time-retarded feedback from environmental degrees of freedom. In this work, we utilize the recently developed quantum dissipation with minimally extended state space approach to address reservoir induced long-ranged temporal correlations in finite size Ising-type spin chains. For thermal reservoirs with ohmic and subohmic spectral density, we simulate the quantum time evolution from finite to zero temperature. The competition between thermal fluctuations, quantum fluctuations, and anti-/ferromagnetic interactions reveals a rich pattern of dynamical phases, including dissipative induced phase transitions and spatiotemporal correlations. [ABSTRACT FROM AUTHOR]

We present a multiscale method coupling the theory of open quantum systems with real-time ab initio treatment of electronic structure to study hot-carrier dynamics in photoexcited plasmonic systems. We combine the Markovian Stochastic Schrödinger equation with an ab initio GW coupled to the Bethe–Salpeter (BSE) equation description of the electronic degrees of freedom, interacting with a metallic nanoparticle modeled classically according to the polarizable continuum model. We apply this methodology to study the effect of relaxation (T1) and pure dephasing (T2) times on the hot-carrier dynamics in a system composed of a quantum portion described at GW/BSE level, i.e., a CHO fragment adsorbed on a vertex of a rhodium nanocube, and of the rest of the nanocube, treated classically, when irradiated with a 2.7 eV light pulse, inspired by the experimental results on plasmon-driven CO2 photoreduction. A net hole injection from rhodium to CHO is observed, with and without the classical portion of the nanocube. The nanocube effect is to enhance the generated charge population by two orders of magnitude. The nonradiative decay, via a relaxation time T1 based on the energy-gap law, produces a rapid decrease of the charge population. Results with T2 only show that a charge injection retarded with respect to the pulse, which is present in the coherent dynamics, disappears when coherence is erased. [ABSTRACT FROM AUTHOR]

Generating and manipulating magnetic skyrmions at ultrafast time scales is essential for future skyrmion-based racetrack memory and logic gate applications. Using the atomistic spin dynamics simulations, we demonstrate the nucleation of ultra-small field-free magnetic skyrmions in amorphous GdFeCo at picosecond time scales by femtosecond laser heating. The ultrafast nature of laser heating and subsequent cooling from a high-temperature state is crucial for forming magnetic skyrmion. The magnon localization and magnon coalescence are the key driving mechanisms responsible for stabilizing the magnetic skyrmions at zero-field conditions. The polarization and, hence, the topological charge can be switched by exploiting the all-optical switching observed in GdFeCo. The skyrmion sizes and numbers can be controlled by varying pulse width and fluence of incident laser pulses. Applying an external magnetic field provides an additional degree of freedom to tune the skyrmion radius during the ultrafast optical creation of magnetic skyrmions. Our results provide a detailed understanding of the ultrafast creation of magnetic skyrmions using femtosecond laser pulses, a vital step in advancing next-generation skyrmion-based memory technologies. [ABSTRACT FROM AUTHOR]

The recent discovery of high-temperature superfluorescence in hybrid perovskite thin films has opened new possibilities for harnessing macroscopic quantum phenomena in nanotechnology. This study aimed to elucidate the mechanism that enables high-temperature superfluorescence in these systems. The proposed model describes a quasi-2D Wannier exciton in a thin film that interacts with phonons via the longitudinal optical phonon–exciton Fröhlich interaction. We show that the super-radiant properties of the coherent state in hybrid perovskites are stable against perturbations caused by the longitudinal optical phonon–exciton Fröhlich interaction. Using the multiconfiguration Hartree approach, we derive semiclassical equations of motion for a single-exciton wavefunction, where the vibrational degrees of freedom interact with the Wannier exciton through a mean-field Hartree term. Super-radiance is effectively described by a non-Hermitian term in the Hamiltonian. The analysis was then extended to multiple excited states using the semiclassical Hamiltonian as the basic model. We demonstrate that the ground state of the model exciton Hamiltonian with long-range interactions is a symmetric Dicke super-radiant state, where the Fröhlich interaction is nullified. The additional density matrix-based consideration draws an analogy between this system and stable systems, where the conservation laws determine the nullification of the constant (momentum-independent) decay rate part. In the exciton–phonon system, nullification is associated with the absence of a momentum-independent component in the Wannier exciton–phonon interaction coupling function. [ABSTRACT FROM AUTHOR]

In single vibronic level (SVL) fluorescence experiments, the electronically excited initial state is also excited in one or several vibrational modes. Because computing such spectra by evaluating all contributing Franck–Condon factors becomes impractical (and unnecessary) in large systems, here we propose a time-dependent approach based on Hagedorn wavepacket dynamics. We use Hagedorn functions—products of a Gaussian and carefully generated polynomials—to represent SVL initial states because in systems whose potential is at most quadratic, Hagedorn functions are exact solutions to the time-dependent Schrödinger equation and can be propagated with the same equations of motion as a simple Gaussian wavepacket. Having developed an efficient recursive algorithm to compute the overlaps between two Hagedorn wavepackets, we can now evaluate emission spectra from arbitrary vibronic levels using a single trajectory. We validate the method in two-dimensional global harmonic models by comparing it with quantum split-operator calculations. In addition, we study the effects of displacement, distortion (squeezing), and Duschinsky rotation on SVL fluorescence spectra. Finally, we demonstrate the applicability of the Hagedorn approach to high-dimensional systems on a displaced, distorted, and Duschinsky-rotated harmonic model with 100 degrees of freedom. [ABSTRACT FROM AUTHOR]

Stimuli-responsive macromolecules display large conformational changes during their dynamics, sometimes switching between states. Such a multi-stability is useful for the development of soft functional materials. Here, we introduce a mean-field dynamical density functional theory for a model of responsive colloids to study the nonequilibrium dynamics of a colloidal dispersion in time-dependent external fields, with a focus on the coupling of translational and conformational dynamics during their relaxation. Specifically, we consider soft Gaussian particles with a bimodal size distribution between two confining walls with time-dependent (switching-on and off) external gravitational and osmotic fields. We find a rich relaxation behavior of the systems in excellent agreement with particle-based Brownian dynamics computer simulations. In particular, we find time-asymmetric relaxations of integrated observables (wall pressures, mean size, and liquid center-of-mass) for activation/deactivation of external potentials, respectively, which are tunable by the ratio of translational and conformational diffusion time scales. Our work thus paves the way for studying the nonequilibrium relaxation dynamics of complex soft matter with multiple degrees of freedom and hierarchical relaxations. [ABSTRACT FROM AUTHOR]

We present an approximate semiclassical (SC) framework for mixed quantized dynamics in Wigner phase space in a two-part series. In the first article, we introduced the Adiabatic Hybrid Wigner Dynamics (AHWD) method that allows for a few important "system" degrees of freedom to be quantized using high-level double Herman–Kluk SC theory while describing the rest (the "bath") using classical-limit linearized SC theory. In this second article, we extend our hybrid Wigner dynamics to nonadiabatic processes. The resulting Nonadiabatic Hybrid Wigner Dynamics (NHWD) has two variants that differ in the choice of degrees of freedom to be quantized. Specifically, we introduce NHWD(E) where only the electronic state variables are quantized and the NHWD(V) where both electronic state variables and a handful of strongly coupled nuclear modes are quantized. We show that while NHWD(E) proves accurate for a wide range of scattering models and spin-boson models, systems where a few nuclear modes are strongly coupled to electronic states require NHWD(V) to accurately capture the long-time dynamics. Taken together, we show that AHWD and NHWD represent a new framework for SC simulations of high-dimensional systems with significant quantum effects. [ABSTRACT FROM AUTHOR]

The Wigner phase space formulation of quantum mechanics is a complete framework for quantum dynamic calculations that elegantly highlights connections with classical dynamics. In this series of two articles, building upon previous efforts, we derive the full hierarchy of approximate semiclassical (SC) dynamic methods for adiabatic and non-adiabatic problems in Wigner phase space. In Paper I, focusing on adiabatic single surface processes, we derive the well-known double Herman–Kluk (DHK) approximation for real-time correlation functions in Wigner phase space and connect it to the linearized SC (LSC) approximation through a stationary phase approximation. We exploit this relationship to introduce a new hybrid SC method, termed Adiabatic Hybrid Wigner Dynamics (AHWD) that allows for a few important "system" degrees of freedom (dofs) to be treated at the DHK level, while treating the rest of the dofs (the "bath") at the LSC level. AHWD is shown to accurately capture quantum interference effects in models of coupled oscillators and the decoherence of vibrational probability density of a model I2 Morse oscillator coupled to an Ohmic thermal bath. We show that AHWD significantly mitigates the sign problem and employs reduced dimensional prefactors bringing calculations of complex system–bath problems within the reach of SC methods. Paper II focuses on extending this hybrid SC dynamics to nonadiabatic processes. [ABSTRACT FROM AUTHOR]

A precise dynamical characterization of quantum impurity models with multiple interacting orbitals is challenging. In quantum Monte Carlo methods, this is embodied by sign problems. A dynamical sign problem makes it exponentially difficult to simulate long times. A multi-orbital sign problem generally results in a prohibitive computational cost for systems with multiple impurity degrees of freedom even in static equilibrium calculations. Here, we present a numerically exact inchworm method that simultaneously alleviates both sign problems, enabling simulation of multi-orbital systems directly in the equilibrium or nonequilibrium steady-state. The method combines ideas from the recently developed steady-state inchworm Monte Carlo framework [Erpenbeck et al., Phys. Rev. Lett. 130, 186301 (2023)] with other ideas from the equilibrium multi-orbital inchworm algorithm [Eidelstein et al., Phys. Rev. Lett. 124, 206405 (2020)]. We verify our method by comparison with analytical limits and numerical results from previous methods. [ABSTRACT FROM AUTHOR]

A comparison is made between three simple approximate formulas for the configurational adiabat (i.e., constant excess entropy, sex) lines in a Lennard-Jones (LJ) fluid, one of which is an analytic formula based on a harmonic approximation, which was derived by Heyes et al. [J. Chem. Phys. 159, 224504 (2023)] (analytic isomorph line, AIL). Another is where the density is normalized by the freezing density at that temperature (freezing isomorph line, FIL). It is found that the AIL formula and the average of the freezing density and the melting density ("FMIL") are configurational adiabats at all densities essentially down to the liquid–vapor binodal. The FIL approximation departs from a configurational adiabat in the vicinity of the liquid–vapor binodal close to the freezing line. The self-diffusion coefficient, D, shear viscosity, ηs, and thermal conductivity, λ, in macroscopic reduced units are essentially constant along the AIL and FMIL at all fluid densities and temperatures, but departures from this trend are found along the FIL at high liquid state densities near the liquid–vapor binodal. This supports growing evidence that for simple model systems with no or few internal degrees of freedom, isodynes are lines of constant excess entropy. It is shown that for the LJ fluid, ηs and D can be predicted accurately by an essentially analytic procedure from the high temperature limiting inverse power fluid values (apart from at very low densities), and this is demonstrated quite well also for the experimental argon viscosity. [ABSTRACT FROM AUTHOR]

We employ the exact-factorization formalism to study the coupled dynamics of photons, electrons, and nuclei at the quantum mechanical level, proposing illustrative examples of model situations of nonadiabatic dynamics and spontaneous emission of electron–nuclear systems in the regime of strong light–matter coupling. We make a particular choice of factorization for such a multi-component system, where the full wavefunction is factored as a conditional electronic amplitude and a marginal photon–nuclear amplitude. Then, we apply the coupled-trajectory mixed quantum–classical (CTMQC) algorithm to perform trajectory-based simulations, by treating photonic and nuclear degrees of freedom on equal footing in terms of classical-like trajectories. The analysis of the time-dependent potentials of the theory along with the assessment of the performance of CTMQC allows us to point out some limitations of the current approximations used in CTMQC. Meanwhile, comparing CTMQC with other trajectory-based algorithms, namely multi-trajectory Ehrenfest and Tully surface hopping, demonstrates the better quality of CTMQC predictions. [ABSTRACT FROM AUTHOR]

The pyrite phase of sodium superoxide, NaO 2 , is studied using equilibrium molecular dynamics simulations and lattice dynamics calculations to understand the impacts of static disorder and dynamic disorder on its thermal conductivity. Three structural regimes are observed based on the rotational dynamics and orientations of O 2 − ions. At low temperatures, where the O 2 − ions librate and the system is fully ordered, thermal conductivity exhibits a crystal-like temperature dependence, decreasing with increasing temperature. As temperature increases, the static disorder regime emerges, where the O 2 − ions transition between different orientations on a time scale larger than the librational period. In this regime, the thermal conductivity continues to decrease and then becomes temperature independent. At higher temperatures, where the O 2 − ions freely rotate, the system is dynamically disordered and the thermal conductivity is temperature independent, as in an amorphous solid. Using instantaneous normal mode analysis and Allen–Feldman theory, 80% of the thermal conductivity in the dynamic disorder regime is attributed to diffusons, vibrational modes that are non-propagating and non-localized. When increasing the lattice constant at a constant temperature, transitions from librations to static disorder to dynamic disorder are also observed, with the thermal conductivity decreasing monotonically. The presented methodology can be applied to other crystals with rotational degrees of freedom, offering strategies for the design of thermal conductivity switches that are responsive to external stimuli. [ABSTRACT FROM AUTHOR]

The pyrite phase of sodium superoxide, NaO 2 , is studied using equilibrium molecular dynamics simulations and lattice dynamics calculations to understand the impacts of static disorder and dynamic disorder on its thermal conductivity. Three structural regimes are observed based on the rotational dynamics and orientations of O 2 − ions. At low temperatures, where the O 2 − ions librate and the system is fully ordered, thermal conductivity exhibits a crystal-like temperature dependence, decreasing with increasing temperature. As temperature increases, the static disorder regime emerges, where the O 2 − ions transition between different orientations on a time scale larger than the librational period. In this regime, the thermal conductivity continues to decrease and then becomes temperature independent. At higher temperatures, where the O 2 − ions freely rotate, the system is dynamically disordered and the thermal conductivity is temperature independent, as in an amorphous solid. Using instantaneous normal mode analysis and Allen–Feldman theory, 80% of the thermal conductivity in the dynamic disorder regime is attributed to diffusons, vibrational modes that are non-propagating and non-localized. When increasing the lattice constant at a constant temperature, transitions from librations to static disorder to dynamic disorder are also observed, with the thermal conductivity decreasing monotonically. The presented methodology can be applied to other crystals with rotational degrees of freedom, offering strategies for the design of thermal conductivity switches that are responsive to external stimuli. [ABSTRACT FROM AUTHOR]

The ultrafast quantum dynamics of photophysical processes in complex molecules is an extremely challenging computational problem with a broad variety of fascinating applications in quantum chemistry and biology. Inspired by recent developments in open quantum systems, we introduce a pure-state unraveled hybrid-bath method that describes a continuous environment via a set of discrete, effective bosonic degrees of freedom using a Markovian embedding. Our method is capable of describing both, a continuous spectral density and sharp peaks embedded into it. Thereby, we overcome the limitations of previous methods, which either capture long-time memory effects using the unitary dynamics of a set of discrete vibrational modes or use memoryless Markovian environments employing a Lindblad or Redfield master equation. We benchmark our method against two paradigmatic problems from quantum chemistry and biology. We demonstrate that compared to unitary descriptions, a significantly smaller number of bosonic modes suffices to describe the excitonic dynamics accurately, yielding a computational speed-up of nearly an order of magnitude. Furthermore, we take into account explicitly the effect of a δ-peak in the spectral density of a light-harvesting complex, demonstrating the strong impact of the long-time memory of the environment on the dynamics. [ABSTRACT FROM AUTHOR]

A complex interplay of structural, electronic, and vibrational degrees of freedom underpins the fate of molecular excited states. Organic assemblies exhibit a myriad of excited-state decay processes, such as symmetry-breaking charge separation (SB-CS), excimer (EX) formation, singlet fission, and energy transfer. Recent studies of cofacial and slip-stacked perylene-3,4:9,10-bis(dicarboximide) (PDI) multimers demonstrate that slight variations in core substituents and H- or J-type aggregation can determine whether the system follows an SB-CS pathway or an EX one. However, questions regarding the relative importance of structural properties and molecular vibrations in driving the excited-state dynamics remain. Here, we use a combination of two-dimensional electronic spectroscopy, femtosecond stimulated Raman spectroscopy, and quantum chemistry computations to compare the photophysics of two PDI dimers. The dimer with 1,7-bis(pyrrolidin-1′-yl) substituents (5PDI2) undergoes ultrafast SB-CS from a photoexcited mixed state, while the dimer with bis-1,7-(3′,5′-di-t-butylphenoxy) substituents (PPDI2) rapidly forms an EX state. Examination of their quantum beating features reveals that SB-CS in 5PDI2 is driven by the collective vibronic coupling of two or more excited-state vibrations. In contrast, we observe signatures of low-frequency vibrational coherence transfer during EX formation by PPDI2, which aligns with several previous studies. We conclude that key electronic and structural differences between 5PDI2 and PPDI2 determine their markedly different photophysics. [ABSTRACT FROM AUTHOR]

Boltzmann's H-function H(t), often regarded as an analog of time-dependent entropy, holds a venerable place in the history of science. However, accurate numerical evaluation of H(t) for particles other than atoms is rare. To remove this lacuna, we generalize Boltzmann's H-function to a gas of molecules with orientational degrees of freedom and evaluate H(t) from the time-dependent joint probability distribution function f(p, L, t) for linear (p) and angular (L) momenta, evolving from an initial nonequilibrium state, by molecular dynamics simulations. We consider both prolate- and oblate-shaped particles, interacting via the well-known Gay–Berne potential and obtain the relaxation of the generalized molecular H(t) from initial (t = 0) nonequilibrium states. In the long-time limit, the H function saturates to its exact equilibrium value, which is the sum of translational and rotational contributions to the respective entropies. Both the translational and rotational components of H(t) decay nearly exponentially with time; the rotational component is more sensitive to the molecular shape that enters through the aspect ratio. A remarkable rapid decrease in the rotational relaxation time is observed as the spherical limit is approached, in a way tantalizingly reminiscent of Hu–Zwanzig hydrodynamic prediction with the slip boundary condition. In addition, we obtain H(t) analytically by solving the appropriate translational and rotational Fokker–Planck equation and obtain a modest agreement with simulations. We observe a remarkable signature of translation–rotation coupling as a function of molecular shape, captured through a physically meaningful differential term that quantifies the magnitude of translation–rotation coupling. [ABSTRACT FROM AUTHOR]

Policy-guided Monte Carlo is an adaptive method to simulate classical interacting systems. It adjusts the proposal distribution of the Metropolis–Hastings algorithm to maximize the sampling efficiency, using a formalism inspired by reinforcement learning. In this work, we first extend the policy-guided method to deal with a general state space, comprising, for instance, both discrete and continuous degrees of freedom, and then apply it to a few paradigmatic models of glass-forming mixtures. We assess the efficiency of a set of physically inspired moves whose proposal distributions are optimized through on-policy learning. Compared to conventional Monte Carlo methods, the optimized proposals are two orders of magnitude faster for an additive soft sphere mixture but yield a much more limited speed-up for the well-studied Kob–Andersen model. We discuss the current limitations of the method and suggest possible ways to improve it. [ABSTRACT FROM AUTHOR]

Quantum dynamics simulation and computational spectroscopy serve as indispensable tools for the theoretical understanding of various fundamental physical and chemical processes, ranging from charge transfer to photochemical reactions. When simulating realistic systems, the primary challenge stems from the overwhelming number of degrees of freedom and the pronounced many-body correlations. Here, we present Kylin-V, an innovative quantum dynamics package designed for accurate and efficient simulations of dynamics and spectroscopic properties of vibronic Hamiltonians for molecular systems and their aggregates. Kylin-V supports various quantum dynamics and computational spectroscopy methods, such as time-dependent density matrix renormalization group and our recently proposed single-site and hierarchical mapping approaches, as well as vibrational heat-bath configuration interaction. In this paper, we introduce the methodologies implemented in Kylin-V and illustrate their performances through a diverse collection of numerical examples. [ABSTRACT FROM AUTHOR]

In this work, we present a mixed quantum–classical open quantum system dynamics method for studying rate modifications of ground-state chemical reactions in an optical cavity under vibrational strong-coupling conditions. In this approach, the cavity radiation mode is treated classically with a mean-field nuclear force averaging over the remaining degrees of freedom, both within the system and the environment, which are handled quantum mechanically within the hierarchical equations of motion framework. Using this approach, we conduct a comparative analysis by juxtaposing the mixed quantum–classical results with fully quantum-mechanical simulations. After eliminating spurious peaks that can occur when not using the rigorous definition of the rate constant, we confirm the crucial role of the quantum nature of the cavity radiation mode in reproducing the resonant peak observed in the cavity frequency-dependent rate profile. In other words, it appears necessary to explicitly consider the quantized photonic states in studying reactivity modification in vibrational polariton chemistry (at least for the model systems studied in this work), as these phenomena stem from cavity-induced reaction pathways involving resonant energy exchanges between photons and molecular vibrational transitions. [ABSTRACT FROM AUTHOR]

Saturated packing is a random packing state of particles widely applied in investigating the physicochemical properties of granular materials. Optimizing particle shape to maximize packing density is a crucial challenge in saturated packing research. The known optimal three-dimensional shape is an ellipsoid with a saturated packing density of 0.437 72(51). In this work, we generate saturated packings of three-dimensional asymmetric shapes, including spherocylinders, cones, and tetrahedra, via the random sequential adsorption algorithm and investigate their packing properties. Results show that the optimal shape of asymmetric spherocylinders gives the maximum density of 0.4338(9), while cones achieve a higher value of 0.4398(10). Interestingly, tetrahedra exhibit two distinct optimal shapes with significantly high densities of 0.4789(19) and 0.4769(18), which surpass all previous results in saturated packing. The study of adsorption kinetics reveals that the two optimal shapes of tetrahedra demonstrate notably higher degrees of freedom and faster growth rates of the particle number. The analysis of packing structures via the density pair-correlation function shows that the two optimal shapes of tetrahedra possess faster transitions from local to global packing densities. [ABSTRACT FROM AUTHOR]

The emergence of many efficient optical field modulation methods and planar optical devices is attributed to the continuous research of geometric phase. Nevertheless, conjugate symmetry of the geometric phase limits the multiplexing of metasurfaces. To overcome this limitation, integrating the propagation phase and the geometric phase to achieve self-decouped metasurfaces can effectively double channel capacity. Herein, a more in-depth derivation of Jones matrices is conducted, leveraging two degrees of freedom offered by the propagation phase and the geometric phase. This approach enables complete modulation in circular polarization channels. By designing the phase difference between a fast axis and a slow axis, energy allocation between different channels can be controlled. It means independent complex amplitude modulation is achieved. On this basis, the geometric phase is introduced to realize tri-channel multiplexing metasurfaces. To verify the feasibility of this method, two metasurfaces are designed, including a bifocal metasurface with adjustable energy allocation and a tri-channel multiplexing metasurface. The proposed multifunctional metasurface offers new insights into wavefront multiplexing for communication systems. [ABSTRACT FROM AUTHOR]

Electronic spins provide an additional degree of freedom that can be used in modern spin-based electronic devices. Some benefits of spintronic devices include nonvolatility, energy efficiency, high endurance, and CMOS compatibility, which can be leveraged for data processing and storage applications in today's digital era. To implement such functionalities, controlling and manipulating electron spins is of prime interest. One of the efficient ways of achieving this in spintronics is to use the electric field to control electron spin or magnetism through the voltage-controlled magnetic anisotropy (VCMA) effect. VCMA avoids the movement of charges and significantly reduces the Ohmic loss. This article reviews VCMA-based spintronic devices for magnetic memory applications. First, we briefly discuss the VCMA effect and various mechanisms explaining its physical origin. We then mention various challenges in VCMA that impede it for practical VCMA-based magnetic memory. We review various techniques to address them, such as field-free switching operation, write error rate improvement, widening the operation window, enhancing the VCMA coefficient, and ensuring fast-read operation with low read disturbance. Finally, we draw conclusions outlining the future perspectives. [ABSTRACT FROM AUTHOR]

The use of cavities to impact molecular structure and dynamics has become popular. As cavities, in particular plasmonic nanocavities, are lossy and the lifetime of their modes can be very short, their lossy nature must be incorporated into the calculations. The Lindblad master equation is commonly considered an appropriate tool to describe this lossy nature. This approach requires the dynamics of the density operator and is thus substantially more costly than approaches employing the Schrödinger equation for the quantum wave function when several or many nuclear degrees of freedom are involved. In this work, we compare numerically the Lindblad and Schrödinger descriptions discussed in the literature for a molecular example where the cavity is pumped by a laser. The laser and cavity properties are varied over a range of parameters. It is found that the Schrödinger description adequately describes the dynamics of the polaritons and emission signal as long as the laser intensity is moderate and the pump time is not much longer than the lifetime of the cavity mode. Otherwise, it is demonstrated that the Schrödinger description gradually fails. We also show that the failure of the Schrödinger description can often be remedied by renormalizing the wave function at every step of time propagation. The results are discussed and analyzed. [ABSTRACT FROM AUTHOR]

We optimize the internuclear geometry and electronic structure of a model chiral system to achieve a maximal photoelectron circular dichroism (PECD) in its one-photon ionization by circularly polarized light. The electronic structure calculations are performed by the single center method, while the optimization is done using quantum alchemy employing a Taylor series expansion. Thereby, the effect of bond lengths and uncompensated charge distributions on the chiral response of the model is investigated theoretically in some detail. It is demonstrated that manipulating a chiral asymmetry of the ionic potential may enhance the dichroic parameter (i.e., the PECD) of the randomly oriented model system well beyond β1 = 25%. Furthermore, we demonstrate that quantum alchemy is applicable to PECD despite the unusually strong coupling of spatial and electronic degrees of freedom and discuss the relative impact of the individual degrees of freedom in this model system. We define the necessary conditions for the computational design of PECD for real (non-model) chiral molecules using our approach. [ABSTRACT FROM AUTHOR]

We present a general and practical theoretical framework to investigate how energy is dissipated in open quantum system dynamics. This is performed by quantifying the contributions of individual bath components to the overall dissipation of the system. The framework is based on the Nakajima–Zwanzig projection operator technique, which allows us to express the rate of energy dissipation into a specific bath degree of freedom by using traces of operator products. The approach captures system-bath interactions to all orders, but is based on second-order perturbation theory on the off-diagonal subsystem's couplings and a Markovian description of the bath. The usefulness of our theory is demonstrated by applying it to various models of open quantum systems involving harmonic oscillators or spin baths and connecting the outcomes to existing results such as our previously reported formula derived for locally coupled harmonic baths [Kim and Franco, J. Chem. Phys. 154, 084109 (2021)]. We also prove that the dissipation calculated by our theory rigorously satisfies thermodynamic principles such as energy conservation and detailed balance. Overall, the strategy can be used to develop the theory and simulation of dissipation pathways to interpret and engineer the dynamics of open quantum systems. [ABSTRACT FROM AUTHOR]

Microgels exhibit the ability to undergo reversible swelling in response to shifts in environmental factors that include variations in temperature, concentration, and pH. While several models have been put forward to elucidate specific aspects of microgel swelling and its impact on bulk behavior, a consistent theoretical description that chains throughout the microscopic degrees of freedom with suspension properties and deepens into the full implications of swelling remains a challenge yet to be met. In this work, we extend the mean-field swelling model of microgels from Denton and Tang [J. Chem. Phys. 145, 164901 (2016)] to include the finite extensibility of the polymer chains. The elastic contribution to swelling in the original work is formulated for Gaussian chains. By using the Langevin chain model, we modify this elastic contribution in order to account for finite extensibility effects, which become prominent for microgels containing highly charged polyelectrolytes and short polymer chains. We assess the performance of both elastic models, namely for Gaussian and Langevin chains, comparing against coarse-grained bead–spring simulations of ionic microgels with explicit electrostatic interactions. We examine the applicability scope of the models under a variation of parameters, such as ionization degree, microgel concentration, and salt concentration. The models are also tested against experimental results. This work broadens the applicability of the microgel swelling model toward a more realistic description, which brings advantages when describing the suspensions of nanogels and weak-polyelectrolyte micro-/nanogels. [ABSTRACT FROM AUTHOR]

In molecular dynamics simulations, dynamically consistent coarse-grained (CG) models commonly use stochastic thermostats to model friction and fluctuations that are lost in a CG description. While Markovian, i.e., time-local, formulations of such thermostats allow for an accurate representation of diffusivities/long-time dynamics, a correct description of the dynamics on all time scales generally requires non-Markovian, i.e., non-time-local, thermostats. These thermostats typically take the form of a Generalized Langevin Equation (GLE) determined by a memory kernel. In this work, we use a Markovian embedded formulation of a position-independent GLE thermostat acting independently on each CG degree of freedom. Extracting the memory kernel of this CG model from atomistic reference data requires several approximations. Therefore, this task is best understood as an inverse problem. While our recently proposed approximate Newton scheme allows for the iterative optimization of memory kernels (IOMK), Markovian embedding remained potentially error-prone and computationally expensive. In this work, we present an IOMK-Gauss–Newton scheme (IOMK-GN) based on IOMK that allows for the direct parameterization of a Markovian embedded model. [ABSTRACT FROM AUTHOR]

In this paper, we present a combined experimental and theoretical study that explored the initial sticking of water on cooled surfaces. Specifically, these ultra-high vacuum gas–surface scattering experiments utilized supersonic molecular beam techniques in conjunction with a cryogenically cooled highly oriented pyrolytic graphite crystal, giving control over incident kinematic conditions. The D2O translational energy spanning 300–750 meV, the relative D2O flux, and the incident angle could all be varied independently. Three different experimental measurements were made. One involved measuring the total amount of D2O scattering as a function of surface temperature to determine the onset of sticking under non-equilibrium gas–surface collision conditions. Another measurement used He specular scattering to assess structural and coverage information for the interface during D2O adsorption. Finally, we used time-of-flight (TOF) measurements of the scattered D2O to determine how energy is exchanged with the graphite surface at surface temperatures above and near the conditions needed for gaseous condensation. For comparison and elaboration of the roles that internal degrees of freedom play in this process, we also did similar TOF measurements using another mass 20 incident particle, atomic neon. Enriching this study are precise molecular dynamics simulations that elaborate on gas–surface energy transfer and the roles of molecular degrees of freedom in gas–surface collisional energy exchange processes. This study furthers our fundamental understanding of energy exchange and the onset of sticking and ultimately gaseous condensation for gas–surface encounters occurring under high-velocity flows. [ABSTRACT FROM AUTHOR]

Porous, stacked two-dimensional covalent organic frameworks (2D COFs) bearing semiconducting linkers can support directional charge transfer across adjacent layers of the COF. To better inform the current and possible future design rules for enhancing electron and hole transport in such materials, an understanding of how linker selection and functionalization affects interlayer electronic couplings is essential. We report electronic structure simulation and analysis of electronic couplings across adjacent linker units and to encapsulated species in functionalized electroactive 2D COFs. The detailed dependence of these electronic couplings on interlayer interactions is examined through scans along key interlayer degrees of freedom and through configurational sampling from equilibrium molecular dynamics on semiempirical potential energy surfaces. Beyond affirming the sensitivity of the electronic coupling to interlayer distance and orientation, these studies offer guidance toward linker functionalization strategies for enhancing charge carrier transport in electroactive 2D COFs. [ABSTRACT FROM AUTHOR]

A novel mixed quantum–classical approach to simulating nonadiabatic dynamics of molecules at metal surfaces is presented. The method combines the numerically exact hierarchical equations of motion approach for the quantum electronic degrees of freedom with Langevin dynamics for the classical degrees of freedom, namely, low-frequency vibrational modes within the molecule. The approach extends previous mixed quantum–classical methods based on Langevin equations to models containing strong electron–electron or quantum electronic–vibrational interactions, while maintaining a nonperturbative and non-Markovian treatment of the molecule–metal coupling. To demonstrate the approach, nonequilibrium transport observables are calculated for a molecular nanojunction containing strong interactions. [ABSTRACT FROM AUTHOR]

Two computational approaches for computing the rates of internal conversions in molecular systems where a large set of nuclear degrees of freedom plays a role are discussed and compared. One approach is based on the numerical solution of the time-dependent Schrödinger equation and allows us to include almost the whole set of vibrational coordinates, thanks to the employment of effective procedures for selecting those elements of the Hilbert space which play a significant role in dynamics. The other approach, based on the time-dependent perturbation theory and limited to the use of the harmonic approximation, allows us to include the whole Hilbert space spanned by the vibrational states of the system. The two approaches are applied to the photophysics of azulene, whose anti-Kasha behavior caused by anomalous internal conversion rates is well assessed. The calculated rates for the decays of the first two excited singlet states are in very good agreement with experimental data, indicating the reliability of both methodologies. [ABSTRACT FROM AUTHOR]

Coarse-grained (CG) molecular models greatly reduce the computational cost of simulations allowing for longer and larger simulations, but come with an artificially increased acceleration of the dynamics when compared to the parent atomistic (AA) simulation. This impedes their use for the quantitative study of dynamical properties. During coarse-graining, grouping several atoms into one CG bead not only reduces the number of degrees of freedom but also reduces the roughness on the molecular surfaces, leading to the acceleration of dynamics. The RoughMob approach [M. K. Meinel and F. Müller-Plathe, J. Phys. Chem. B 126(20), 3737–3747 (2022)] quantifies this change in geometry and correlates it to the acceleration by making use of four so-called roughness volumes. This method was developed using simple one-bead CG models of a set of hydrocarbon liquids. Potentials for pure components are derived by the structure-based iterative Boltzmann inversion. In this paper, we find that, for binary mixtures of simple hydrocarbons, it is sufficient to use simple averaging rules to calculate the roughness volumes in mixtures from the roughness volumes of pure components and add a correction term quadratic in the concentration without the need to perform any calculation on AA or CG trajectories of the mixtures themselves. The acceleration factors of binary diffusion coefficients and both self-diffusion coefficients show a large dependence on the overall acceleration of the system and can be predicted a priori without the need for any AA simulations within a percentage error margin, which is comparable to routine measurement accuracies. Only if a qualitatively accurate description of the concentration dependence of the binary diffusion coefficient is desired, very few additional simulations of the pure components and the equimolar mixture are required. [ABSTRACT FROM AUTHOR]

Deposition of Mo S 2 on Au(111) alters the electronic properties of Mo S 2 . In this study, we investigate the free-standing MoS 2 monolayer and the MoS 2 /Au(111) heterostructure, with and without strain, as well as defects of interest in memristive and neuromorphic applications. We focus on the so-called atomristor devices based on monolayer materials that achieve resistive switching characteristics with the adsorption and desorption of metal adatoms. Our study confirms that the formation of midgap states is the primary mechanism behind the resistive switching. Our results show that strain lowers the adsorption/desorption energies of Au+defect structures of interest, leading to more favorable switching energies, but simultaneously reduces the switching ratio between states of differing conductivities. The presence of the Au(111) substrate additionally introduces non-uniform amounts of strain and charge transfer to the MoS 2 monolayer. We propose that the induced strain contributes to the experimentally observed n - to p -type transition and Ohmic to Schottky transition in the MoS 2 monolayer. The charge transfer leads to a permanent polarization at the interface, which can be tuned by strain. Our study has important implications on the role of the electrode as being a source of the observed variability in memristive devices and as an additional degree of freedom for tuning the switching characteristics of the memristor device. [ABSTRACT FROM AUTHOR]

Enhanced sampling algorithms are indispensable when working with highly disconnected multimodal distributions. An important application of these is the conformational exploration of particular internal degrees of freedom of molecular systems. However, despite the existence of many commonly used enhanced sampling algorithms to explore these internal motions, they often rely on system-dependent parameters, which negatively impact efficiency and reproducibility. Here, we present fully adaptive simulated tempering (FAST), a variation of the irreversible simulated tempering algorithm, which continuously optimizes the number, parameters, and weights of intermediate distributions to achieve maximally fast traversal over a space defined by the change in a predefined thermodynamic control variable such as temperature or an alchemical smoothing parameter. This work builds on a number of previously published methods, such as sequential Monte Carlo, and introduces a novel parameter optimization procedure that can, in principle, be used in any expanded ensemble algorithms. This method is validated by being applied on a number of different molecular systems with high torsional kinetic barriers. We also consider two different soft-core potentials during the interpolation procedure and compare their performance. We conclude that FAST is a highly efficient algorithm, which improves simulation reproducibility and can be successfully used in a variety of settings with the same initial hyperparameters. [ABSTRACT FROM AUTHOR]

Among lattice configurations of densely packed hard ellipses, Monte Carlo simulations are used to identify the so-called parallel and diagonal lattices as the two favorable states. The free energies of these two states are computed for several system sizes employing the Einstein crystal method. An accurate calculation of the free energy difference between the two states reveals the parallel lattice as the state with the lowest free energy. The origin of the entropic difference between the two states is further elucidated by assessing the roles of the translational and rotational degrees of freedom. [ABSTRACT FROM AUTHOR]

Entangled two-photon absorption (ETPA) could form the basis of nonlinear quantum spectroscopy at very low photon fluxes, since, at sufficiently low photon fluxes, ETPA scales linearly with the photon flux. When different pairs start to overlap temporally, accidental coincidences are thought to give rise to a "classical" quadratic scaling that dominates the signal at large photon fluxes and, thus, recovers a supposedly classical regime, where any quantum advantage is thought to be lost. Here, we scrutinize this assumption and demonstrate that quantum-enhanced absorption cross sections can persist even for very large photon numbers. To this end, we use a minimal model for quantum light, which can interpolate continuously between the entangled pair and a high-photon-flux limit, to analytically derive ETPA cross sections and the intensity crossover regime. We investigate the interplay between spectral and spatial degrees of freedom and how linewidth broadening of the sample impacts the experimentally achievable enhancement. [ABSTRACT FROM AUTHOR]

Two-dimensional (2D) transition metal dichalcogenides have emerged as promising quantum functional blocks benefitting from their unique combination of spin, valley, and layer degrees of freedom, particularly for the tremendous flexibility of moiré superlattices formed by van der Waals stacking. These degrees of freedom coupled with the enhanced Coulomb interaction in 2D structures allow excitons to serve as on-chip information carriers. However, excitons are spatially circumscribed due to their low mobility and limited lifetime. One way to overcome these limitations is through the coupling of excitons with surface plasmon polaritons (SPPs), which facilitates an interaction between remote quantum states. Here, we showcase the successful coupling of SPPs with interlayer excitons in molybdenum diselenide/tungsten diselenide heterobilayers. Our results indicate that the valley polarization can be efficiently transferred to SPPs, enabling preservation of polarization information even after propagating tens of micrometers. [ABSTRACT FROM AUTHOR]