Search

Understanding the time evolution of strongly correlated open quantum systems (OQSs) in response to perturbations (quenches) is of fundamental importance to the precise control of quantum devices. It is, however, rather challenging in multi-impurity quantum systems because such evolution often involves multiple intricate dynamical processes. In this work, we apply the numerically exact hierarchical equations of motion approach to explore the influence of two different types of perturbations, i.e., sudden swapping of the energy levels of impurity systems and activating the inter-impurity spin-exchange interaction, on the dissipation dynamics of the Kondo-correlated two-impurity Anderson model over a wide range of energetic parameters. By evaluating the time-dependent impurity spectral function and other system properties, we analyze the time evolution of the Kondo state in detail and conclude a phenomenologically scaling relation for Kondo dynamics driven by these perturbations. The evolutionary scaling relationship is not only related to the Kondo characteristic energy TK but also significantly affected by the simultaneous non-Kondo dynamic characteristic energy. We expect these results will inspire subsequent theoretical studies on the dynamics of strongly correlated OQSs. [ABSTRACT FROM AUTHOR]

We describe the new method that can be useful for calculation of the excitation dynamics in large molecular arrays that can be split into compartments with weak exciton coupling between them. In this method, the dynamics within each compartment is evaluated nonperturbatively using hierarchical equations of motion (HEOM), whereas transfers between the exciton states belonging to different compartments are treated by the generalized Förster (gF) theory. In a combined HEOM-gF approach, the number of equations increases linearly when adding new compartments as opposed to pure HEOM, where a depth of hierarchy exhibits strong non-linear grows when scaling the total number of molecules. Comparing the combined HEOM-gF method with an exact HEOM solution enabled us to estimate the parameters corresponding to a validity range of the proposed theory. The possibility of using the method for modeling of energy transfers in photosynthetic antenna supercomplexes is discussed. [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]

Transport and dispersion of active particles in structured environments, such as corrugated channels and porous media, are important for the understanding of both natural and engineered active systems. Owing to their continuous self-propulsion, active particles exhibit rectified transport under spatially asymmetric confinement. While progress has been made in experiments and particle-based simulations, a theoretical understanding of the effective long-time transport dynamics in spatially periodic geometries remains less developed. In this paper, we apply generalized Taylor dispersion theory to analyze the long-time effective transport dynamics of active Brownian particles (ABPs) in periodic channels and fields. We show that the long-time transport behavior is governed by an effective advection–diffusion equation. The derived macrotransport equations allow us to characterize the average drift and effective dispersion coefficient. For the case of ABPs subject to a no-flux boundary condition at the channel wall, we show that regardless of activity, the average drift is given by the net diffusive flux along the channel. For ABPs, their activity is the driving mechanism that sustains a density gradient, which ultimately leads to rectified motion along the channel. Our continuum theory is validated against direct Brownian dynamics simulations of the Langevin equations governing the motion of each ABP. [ABSTRACT FROM AUTHOR]

This study benchmarks the newly re-implemented single-reference excited-state methods, IP-EOM-CCSD, EA-EOM-CCSD, and STEOM-CCSD, in ORCA6.0, with a focus on open-shell systems. We compare STEOM against EOM-CCSD, CC3, and CCSDT across a range of systems, including small organic radicals, hydrated transition metal (TM) ions, and TM diatomic systems with both closed and open-shell configurations. For organic radicals, STEOM and EOM-CCSD show comparable performance, aligning closely with CC3 and CCSDT results. In the case of hydrated TM ions, IP-EOM closely matches DLPNO-CCSD results, while deviations from DLPNO-CCSD(T) are consistent. For open-shell TM systems, IP-EOM exhibits a blueshift relative to both the DLPNO-CCSD methods, while EA-EOM-CCSD shows better agreement. When comparing STEOM and CC3 to CCSDT, STEOM shows slightly larger deviations in closed-shell systems but shows excellent agreement in open-shell systems. Computational efficiency is also assessed, revealing a significant speedup in ORCA 6.0 compared to ORCA 5.0, with optimizations improving computation times. This study provides valuable insights into the performance and efficiency of STEOM in various chemical environments, highlighting its strengths and limitations. [ABSTRACT FROM AUTHOR]

We investigate the quantum dynamics of a spin coupling to a bath of independent spins via the dissipaton equation of motion (DEOM) approach. The bath, characterized by a continuous spectral density function, is composed of spins that are independent level systems described by the s u (2) Lie algebra, representing an environment with a large magnitude of anharmonicity. Based on the previous work by Suarez and Silbey [J. Chem. Phys. 95, 9115 (1991)] and by Makri [J. Chem. Phys. 111, 6164 (1999)] that the spin bath can be mapped to a Gaussian environment under its linear response limit, we use the time-domain Prony fitting decomposition scheme to the bare–bath time correlation function (TCF) given by the bosonic fluctuation–dissipation theorem to generate the exponential decay basis (or pseudo modes) for DEOM construction. The accuracy and efficiency of this strategy have been explored by a variety of numerical results. We envision that this work provides new insights into extending the hierarchical equations of motion and DEOM approach to certain types of anharmonic environments with arbitrary TCF or spectral density. [ABSTRACT FROM AUTHOR]

An important class of methods for modeling dynamics in open quantum systems is based on the well-known influence functional (IF) approach to solving path-integral equations of motion. Within this paradigm, path-filtering schemes based on the removal of IF elements that fall below a certain threshold aim to reduce the effort needed to calculate and store the influence functional, making very challenging simulations possible. A filtering protocol of this type is considered acceptable as long as the simulation remains mathematically stable. This, however, does not guarantee that the approximated dynamics preserve the physics of the simulated process. In this paper, we explore the possibility of training Generative Flow Networks (GFlowNets) to produce filtering protocols while optimizing for mathematical stability and for physical accuracy. Trained using the trajectory balance objective, the model produces sets of paths to be added to a truncated initial set; it is rewarded if the combined set of paths gives rise to solutions in which the trace of the density matrix is conserved, the populations remain real, and the dynamics approach the exact reference. Using a simple two-level system coupled to a dissipative reservoir, we perform proof-of-concept simulations and demonstrate the elegant and surprising filtering solutions proposed by the GFlowNet. [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]

Calculating molecular properties using quantum devices can be performed through the quantum linear response (qLR) or, equivalently, the quantum equation of motion (qEOM) formulations. Different parameterizations of qLR and qEOM are available, namely naïve, projected, self-consistent, and state-transfer. In the naïve and projected parameterizations, the metric is not the identity, and we show that it depends on redundant orbital rotations. This dependency may lead to divergences in the excitation energies for certain choices of the redundant orbital rotation parameters in an idealized noiseless setting. Furthermore, this leads to a significant variance when calculations include statistical noise from finite quantum sampling. [ABSTRACT FROM AUTHOR]

We introduce a Python package based on matrix product states (MPS) to simulate both the time-dependent Schrödinger equation (TDSE) and the hierarchical equations of motion (HEOM). The wave function in the TDSE or the reduced density operator/auxiliary density operators in the HEOM are represented using MPS. A matrix product operator (MPO) is then constructed to represent the Hamiltonian in the TDSE or the generalized Liouvillian in the HEOM. The fourth-order Runge–Kutta method and the time-dependent variational principle are used to propagate the MPS. Several examples, including the nonadiabatic interconversion dynamics of the pyrazine molecule, excitation energy transfer dynamics in molecular aggregates and photosynthetic light-harvesting complexes, the spin-boson model, a laser driven two-state model, the Holstein model, and charge transport in the Anderson impurity model, are presented to demonstrate the capability of the package. [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]

Five time-dependent orbital optimized coupled-cluster methods, of which four can converge to the time-dependent complete active space self-consistent-field method, are presented for fermion-mixtures with arbitrary fermion kinds and numbers. Truncation schemes maintaining the intragroup orbital rotation invariance, as well as equations of motion of coupled-cluster (CC) amplitudes and orbitals, are derived. Present methods are compact CC-parameterization alternatives to the time-dependent multiconfiguration self-consistent-field method for systems consisting of arbitrarily different kinds and numbers of interacting fermions. Theoretical analysis of applications of present methods to various chemical systems is reported. [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]

In quantum chemistry, Lowest Unoccupied Molecular Orbital (LUMO) is important for studying various chemical processes, including photochemical reactions, electron attached states, and electron excites states. Recently, an effective method has been introduced that involves the use of the Parametric Equation of Motion (PEM) in conjunction with the nuclear charge stabilization method for precise identification of true LUMO. However, the inclusion of extra diffuse functions in the basis set, which is necessary for describing electron-attached and electron-excited states, can cause issues due to the presence of the same symmetry states, leading to avoided crossing. Identifying the true LUMO among these avoided crossings is challenging due to the mixing of states and the exchange of their orbital character. This article introduces a modification of the PEM to identify the true LUMO by preventing the stabilization of specific states involved in avoided crossings. The present method is highly effective and requires minimal computational cost. [ABSTRACT FROM AUTHOR]

Evaluation of the charge transport property of organic semiconductors requires exact quantum dynamics simulation of large systems. We present a numerically nearly exact approach to investigate carrier transport dynamics in organic semiconductors by extending the non-Markovian stochastic Schrödinger equation with complex frequency modes to a forward–backward scheme and by solving it using the matrix product state (MPS) approach. By utilizing the forward–backward formalism for noise generation, the bath correlation function can be effectively treated as a temperature-independent imaginary part, enabling a more accurate decomposition with fewer complex frequency modes. Using this approach, we study the carrier transport and mobility in the one-dimensional Peierls model, where the nonlocal electron–phonon interaction is taken into account. The reliability of this approach was validated by comparing carrier diffusion motion with those obtained from the hierarchical equations of motion method across various parameter regimes of the phonon bath. The efficiency was demonstrated by the modest virtual bond dimensions of MPS and the low scaling of the computational time with the system size. [ABSTRACT FROM AUTHOR]

We here present the normal dynamics technique, which recasts the Newton's equations of motion in terms of phonon normal modes by exploiting a proper sampling of the reciprocal space. After introducing the theoretical background, we discuss how the reciprocal space sampling enables us to (i) obtain a computational speedup by selecting which and how many wave vectors of the Brillouin zone will be considered and (ii) account for distortions realized across large atomic distances without the use of large simulation cells. We implemented the approach into an open-source code, which we used to present three case studies: in the first one, we elucidate the general strategy for the sampling of the reciprocal space; in the second one, we illustrate the potential of the approach by studying the stabilization effect of temperature in α-uranium; and in the last one, we investigate the characterization of Raman spectra at different temperatures in MoS2/MX2 transition metal dichalcogenide heterostructures. Finally, we discuss how the procedure is general and can be used to simulate periodic, semiperiodic, and finite systems such as crystals, slabs, nanoclusters, or molecules. [ABSTRACT FROM AUTHOR]

We theoretically study the conformational and dynamical properties of semiflexible active polar ring polymers under linear shear flow. A ring is described as a continuous semiflexible Gaussian polymer with a tangential active force of a constant density along its contour. The linear but non-Hermitian equation of motion is solved using an eigenfunction expansion, which yields activity-independent, but shear-rate-dependent, relaxation times and activity-dependent frequencies. As a consequence, the ring's stationary-state properties are independent of activity, and its conformations and rheological properties are equal to those of a passive ring under shear. The presence of characteristic time scales by relaxation and the activity-dependent frequencies give rise to a particular dynamical behavior. A tank-treading-like motion emerges for long relaxation times and high activities, specifically for stiff rings. In the case of very flexible polymers, the relaxation behavior dominates over activity contributions suppressing tank-treading. Shear strongly affects the crossover from a tank-treading to a relaxation-dominated dynamics, and the ring polymer exhibits tumbling motion at high shear rates. This is reflected in the tumbling frequency, which displays two shear-rate dependent regimes, with an activity-dependent plateau at low shear rates followed by a power-law regime with increasing tumbling frequency for high shear rates. [ABSTRACT FROM AUTHOR]

The motional narrowing effect has been extensively studied for cavity exciton–polariton systems in recent decades both experimentally and theoretically, which is featured by (1) the subaverage behavior and (2) the asymmetric linewidths for the upper polariton and the lower polariton. However, a minimal theoretical model that is clear and adequate to address all these effects as well as the linewidth scaling relations remains missing. In this work, based on the single mode 1D Holstein–Tavis–Cummings (HTC) model, we studied the motional narrowing effect of the polariton linear absorption spectra via both semi-analytic derivations and numerically exact quantum dynamics simulations using the hierarchical equations of motion approach. The results reveal that under collective light–matter coupling between a cavity mode and N molecules, the polariton linewidth scales as 1 / N under the slow limit, while scales as 1/N under the fast limit, due to the polaron decoupling effect. Furthermore, by varying the detunings, the polariton linewidths exhibit significant motional narrowing, covering both characters mentioned above. Our analytic linewidth expressions [Eqs. (34) and (35)] agree well with the numerical exact simulations in all the parameter regimes we explored. These results indicate that the physics of motional narrowing is adequately accounted for by the single-mode 1D HTC model. We envision that both the numerical results and the analytic polariton linewidths expression presented in this work will offer great theoretical value for providing a better understanding of the exciton–polariton motional narrowing based on the HTC model. [ABSTRACT FROM AUTHOR]

Strong coupling between the molecular system and photon inside the cavity generates polaritons, which can alter reaction rates by orders of magnitude. In this work, we benchmark the surface hopping method to simulate non-adiabatic dynamics in a cavity. The comparison is made against a numerically exact method (the hierarchical equations of motion) for a model system investigating excitonic energy transfer for a broad range of parameters. Surface hopping captures the effects of the radiation mode well, both at resonance and off-resonance. We have further investigated parameters that can increase or decrease the rate of population transfer, and we find that surface hopping in general can capture both effects well. Finally, we show that the dipole self-energy term within our parameter regime does not significantly affect the system's dynamics. [ABSTRACT FROM AUTHOR]

Conical intersections (CIs) play a crucial role in photochemical reactions, offering an efficient channel for ultrafast non-adiabatic relaxation of excited states. This significantly influences the reaction pathways and the resulting products. In this work, we utilize the non-Markovian stochastic Schrödinger equation with complex modes method to explore the dynamics of electronic transitions through conical intersections (CIs) in pyrazine. The linear vibronic coupling model serves as the foundational framework, incorporating both intra-state and inter-state electron–vibrational interactions. The dynamics of the excited electronic transitions are analyzed across varying strengths of system–bath coupling and different bath relaxation times. The accuracy of this method is demonstrated by comparing its predictions with those from the hierarchical equations of motion method. [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]

Voltage distribution in sub-cellular micro-domains such as neuronal synapses, small protrusions, or dendritic spines regulates the opening and closing of ionic channels, energy production, and thus, cellular homeostasis and excitability. Yet how voltage changes at such a small scale in vivo remains challenging due to the experimental diffraction limit, large signal fluctuations, and the still limited resolution of fast voltage indicators. Here, we study the voltage distribution in nano-compartments using a computational approach based on the Poisson–Nernst–Planck equations for the electro-diffusion motion of ions, where inward and outward fluxes are generated between channels. We report a current–voltage (I–V) logarithmic relationship generalizing Nernst law that reveals how the local membrane curvature modulates the voltage. We further find that an influx current penetrating a cellular electrolyte can lead to perturbations from tens to hundreds of nanometers deep, depending on the local channel organization. Finally, we show that the neck resistance of dendritic spines can be completely shunted by the transporters located on the head boundary, facilitating ionic flow. To conclude, we propose that voltage is regulated at a subcellular level by channel organization, membrane curvature, and narrow passages. [ABSTRACT FROM AUTHOR]

The projected complex absorbing potential (CAP) technique is one of the methods that allow one to extend the bound state methods for computing resonances' energies and widths. Here, we explore the accuracy of the potential energy curves generated with different electronic structure theory methods in combination with the projected CAP technique by considering resonant vibrational excitation (RVE) of N2 by electron impact as a model process. We report RVE cross sections computed using the boomerang model with potential energy curves obtained with CAP-based extended multistate complete active space perturbation theory (XMS-CASPT2) and equation of motion coupled-cluster method for electron attachment with single and double substitution (EOM-EA-CCSD) methods. We also compare potential energy curves computed with several electronic structure methods, including XMS-CASPT2, EOM-EA-CCSD, multireference configuration interaction with singles (MR-CIS) and singles and doubles (MR-CISD). A good agreement is observed between the experiment and simulated RVE cross sections obtained with the potential energy curves generated with XMS-CASPT2 and EOM-EA-CCSD methods, thus highlighting the potential of the projected CAP technique combined with accurate electronic structure methods for dynamical simulations of the processes that proceed through metastable electronic states. [ABSTRACT FROM AUTHOR]

The projected complex absorbing potential (CAP) technique is one of the methods that allow one to extend the bound state methods for computing resonances' energies and widths. Here, we explore the accuracy of the potential energy curves generated with different electronic structure theory methods in combination with the projected CAP technique by considering resonant vibrational excitation (RVE) of N2 by electron impact as a model process. We report RVE cross sections computed using the boomerang model with potential energy curves obtained with CAP-based extended multistate complete active space perturbation theory (XMS-CASPT2) and equation of motion coupled-cluster method for electron attachment with single and double substitution (EOM-EA-CCSD) methods. We also compare potential energy curves computed with several electronic structure methods, including XMS-CASPT2, EOM-EA-CCSD, multireference configuration interaction with singles (MR-CIS) and singles and doubles (MR-CISD). A good agreement is observed between the experiment and simulated RVE cross sections obtained with the potential energy curves generated with XMS-CASPT2 and EOM-EA-CCSD methods, thus highlighting the potential of the projected CAP technique combined with accurate electronic structure methods for dynamical simulations of the processes that proceed through metastable electronic states. [ABSTRACT FROM AUTHOR]

In non-perturbative non-Markovian open quantum systems, reaching either low temperatures with the hierarchical equations of motion (HEOM) or high temperatures with the Thermalized Time Evolving Density Operator with Orthogonal Polynomials Algorithm (T-TEDOPA) formalism in Hilbert space remains challenging. We compare different ways of modeling the environment. Sampling the Fourier transform of the bath correlation function, also called temperature dependent spectral density, proves to be very effective. T-TEDOPA [Tamascelli et al., Phys. Rev. Lett. 123, 090402 (2019)] uses a linear chain of oscillators with positive and negative frequencies, while HEOM is based on the complex poles of an optimized rational decomposition of the temperature dependent spectral density [Xu et al., Phys. Rev. Lett. 129, 230601 (2022)]. Resorting to the poles of the temperature independent spectral density and of the Bose function separately is an alternative when the problem due to the huge number of Bose poles at low temperatures is circumvented. Two examples illustrate the effectiveness of the HEOM and T-TEDOPA approaches: a benchmark pure dephasing case and a two-bath model simulating the dynamics of excited electronic states coupled through a conical intersection. We show the efficiency of T-TEDOPA to simulate dynamics at a finite temperature by using either continuous spectral densities or only all the intramolecular oscillators of a linear vibronic model calibrated from ab initio data of a phenylene ethynylene dimer. [ABSTRACT FROM AUTHOR]

In the previous paper [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111-1–214111-13 (2024)], we developed a theory called MQME-D, which allows us to decompose the overall energy dissipation process in open quantum system dynamics into contributions by individual components of the bath when the subsystem dynamics is governed by a Markovian quantum master equation (MQME). Here, we contrast the predictions of MQME-D against the numerically exact results obtained by combining hierarchical equations of motion (HEOM) with a recently reported protocol for monitoring the statistics of the bath. Overall, MQME-D accurately captures the contributions of specific bath components to the overall dissipation while greatly reducing the computational cost compared to exact computations using HEOM. The computations show that MQME-D exhibits errors originating from its inherent Markov approximation. We demonstrate that its accuracy can be significantly increased by incorporating non-Markovianity by exploiting time scale separations (TSS) in different components of the bath. Our work demonstrates that MQME-D combined with TSS can be reliably used to understand how energy is dissipated in realistic open quantum system dynamics. [ABSTRACT FROM AUTHOR]

Static and time-dependent quantum-mechanical approaches have been employed in the literature to characterize the physics of light-emitting molecules and nanostructures. However, the electromagnetic emission induced by an input current has remained beyond the realm of molecular simulations. This is the challenge addressed here with the help of an equation of motion for the density matrix coupled to a photon bath based on a Redfield formulation. This equation is evolved within the framework of the driven-Liouville von Neumann approach, which incorporates open boundaries by introducing an applied bias and a circulating current. The dissipated electromagnetic power can be computed in this context from the time derivative of the energy. This scheme is applied in combination with a self-consistent tight-binding Hamiltonian to investigate the effects of bias and molecular size on the electroluminescence of metallic and semiconducting chains. For the latter, a complex interplay between bias and molecular length is observed: there is an optimal number of atoms that maximizes the emitted power at high voltages but not at low ones. This unanticipated behavior can be understood in terms of the band bending produced along the semiconducting chain, a phenomenon that is captured by the self-consistency of the method. A simple analytical model is proposed that explains the main features revealed by the simulations. The methodology, applied here at a self-consistent tight-binding level but extendable to more sophisticated Hamiltonians such as density functional tight binding and time dependent density functional theory, promises to be helpful for quantifying the power and quantum efficiency of nanoscale electroluminescent devices. [ABSTRACT FROM AUTHOR]

We develop a quasiparticle approach to capture the dynamics of open quantum systems coupled to bosonic thermal baths of arbitrary complexity based on the Hierarchical Equations of Motion (HEOM). This is done by generalizing the HEOM dynamics and mapping it into that of the system in interaction with a few bosonic fictitious quasiparticles that we call bexcitons. Bexcitons arise from a decomposition of the bath correlation function into discrete features. Specifically, bexciton creation and annihilation couple the auxiliary density matrices in the HEOM. The approach provides a systematic strategy to construct exact quantum master equations that include the system–bath coupling to all orders even for non-Markovian environments. Specifically, by introducing different metrics and representations for the bexcitons it is possible to straightforwardly generate different variants of the HEOM, demonstrating that all these variants share a common underlying quasiparticle picture. Bexcitonic properties, while unphysical, offer a coarse-grained view of the correlated system–bath dynamics and its numerical convergence. For instance, we use it to analyze the instability of the HEOM when the bath is composed of underdamped oscillators and show that it leads to the creation of highly excited bexcitons. The bexcitonic picture can also be used to develop more efficient approaches to propagate the HEOM. As an example, we use the particle-like nature of the bexcitons to introduce mode-combination of bexcitons in both number and coordinate representation that uses the multi-configuration time-dependent Hartree to efficiently propagate the HEOM dynamics. [ABSTRACT FROM AUTHOR]

The time-dependent quantum transportation through a metal/polymer/metal system is theoretically investigated on the basis of a Su–Schrieffer–Heeger model combined with the hierarchical equations of motion formalism. Using a non-adiabatic dynamical method, the evolution of the electron subspace and lattice atoms with time can be obtained. It is found that the calculated transient currents vary with time and reach stable values after a response time under the bias voltages. However, the stable current as the system reaches its dynamical steady state exhibits a discrepancy between two sweep directions of the bias voltage, which results in pronounced electrical hysteresis loops in the current–voltage curve. By analyzing the evolution of instantaneous energy eigenstates, the occupation number of the instantaneous eigenstates, and the lattice of the polymer, we show that the formation of excitons and the delay of their annihilation are responsible for the hysteretic current–voltage characteristics, where electron–phonon interactions play the key factor. Furthermore, the hysteresis width and amplitude can also be modulated by the strength of the electron–phonon coupling, level-width broadening function, and temperature. We hope these results about past condition-dependent switching performance at a sweep voltage can provide further insight into some of the basic issues of interest in hysteresis processes in conducting polymers. [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]

A time-dependent vibrational electronic coupled-cluster (VECC) approach is proposed to simulate photo-electron/UV-VIS absorption spectra as well as time-dependent properties for non-adiabatic vibronic models, going beyond the Born–Oppenheimer approximation. A detailed derivation of the equations of motion and a motivation for the ansatz are presented. The VECC method employs second-quantized bosonic construction operators and a mixed linear and exponential ansatz to form a compact representation of the time-dependent wave-function. Importantly, the method does not require a basis set, has only a few user-defined inputs, and has a classical (polynomial) scaling with respect to the number of degrees of freedom (of the vibronic model), resulting in a favorable computational cost. In benchmark applications to small models and molecules, the VECC method provides accurate results compared to multi-configurational time-dependent Hartree calculations when predicting short-time dynamical properties (i.e., photo-electron/UV–VIS absorption spectra) for non-adiabatic vibronic models. To illustrate the capabilities, the VECC method is also successfully applied to a large vibronic model for hexahelicene with 14 electronic states and 63 normal modes, developed in the group by Aranda and Santoro [J. Chem. Theory Comput. 17, 1691, (2021)]. [ABSTRACT FROM AUTHOR]

Controlling molecular reactivity by shaped laser pulses is a long-standing goal in chemistry. Here, we suggest a direct optimal control approach that combines external pulse optimization with other control parameters arising in the upcoming field of vibro-polaritonic chemistry for enhanced controllability. The direct optimal control approach is characterized by a simultaneous simulation and optimization paradigm, meaning that the equations of motion are discretized and converted into a set of holonomic constraints for a nonlinear optimization problem given by the control functional. Compared with indirect optimal control, this procedure offers great flexibility, such as final time or Hamiltonian parameter optimization. A simultaneous direct optimal control theory will be applied to a model system describing H-atom transfer in a lossy Fabry–Pérot cavity under vibrational strong coupling conditions. Specifically, optimization of the cavity coupling strength and, thus, of the control landscape will be demonstrated. [ABSTRACT FROM AUTHOR]

The operations of current quantum computers are still significantly affected by decoherence caused by interaction with the environment. In this work, we employ the non-perturbative hierarchical equations of motion (HEOM) method to simulate the operation of model quantum computers and reveal the effects of dissipation on the entangled quantum states and on the performance of well-known quantum algorithms. Multi-qubit entangled states in Shor's factorizing algorithm are first generated and propagated using the HEOM. It is found that the failure of factorization is accompanied by a loss of fidelity and mutual information. An important challenge in using the HEOM to simulate quantum computers in a dissipative environment is how to efficiently treat systems with many qubits. We propose a two-dimensional tensor network scheme for this problem and demonstrate its capability by simulating a one-dimensional random circuit model with 21 qubits. [ABSTRACT FROM AUTHOR]

Scaffolded molecular networks are important building blocks in biological pigment–protein complexes, and DNA nanotechnology allows analogous systems to be designed and synthesized. System–environment interactions in these systems are responsible for important processes, such as the dissipation of heat and quantum information. This study investigates the role of nanoscale molecular parameters in tuning these vibronic system–environment dynamics. Here, genetic algorithm methods are used to obtain nanoscale parameters for a DNA-scaffolded chromophore network based on comparisons between its calculated and measured optical spectra. These parameters include the positions, orientations, and energy level characteristics within the network. This information is then used to compute the dynamics, including the vibronic population dynamics and system–environment heat currents, using the hierarchical equations of motion. The dissipation of quantum information is identified by the system's transient change in entropy, which is proportional to the heat currents according to the second law of thermodynamics. These results indicate that the dissipation of quantum information is highly dependent on the particular nanoscale characteristics of the molecular network, which is a necessary first step before gleaning the systematic optimization rules. Subsequently, the I-concurrence dynamics are calculated to understand the evolution of the vibronic system's quantum entanglement, which are found to be long-lived compared to these system–bath dissipation processes. [ABSTRACT FROM AUTHOR]

A first implementation of analytic gradients for spinor-based relativistic equation-of-motion coupled-cluster singles and doubles method using an exact two-component Hamiltonian augmented with atomic mean-field spin–orbit integrals is reported. To demonstrate its applicability, we present calculations of equilibrium structures and harmonic vibrational frequencies for the electronic ground and excited states of the radium mono-amide molecule (RaNH2) and the radium mono-methoxide molecule (RaOCH3). Spin–orbit coupling is shown to quench Jahn–Teller effects in the first excited state of RaOCH3, resulting in a C3v equilibrium structure. The calculations also show that the radium atoms in these molecules serve as efficient optical cycling centers. [ABSTRACT FROM AUTHOR]

In this work, we formulate the equations of motion corresponding to the Hermitian operator method in the framework of the doubly occupied configuration interaction space. The resulting algorithms turn out to be considerably simpler than the equations provided by that method in more conventional spaces, enabling the determination of excitation energies in N-electron systems under an affordable polynomial computational cost. The implementation of this technique only requires to know the elements of low-order reduced density matrices of an N-electron reference state, which can be obtained from any approximate method. We contrast our procedure against the reduced Bardeen–Cooper–Schrieffer and Richardson–Gaudin–Kitaev integrable models, pointing out the reliability of our proposal. [ABSTRACT FROM AUTHOR]

Interpreting high-dimensional data from molecular dynamics simulations is a persistent challenge. In this paper, we show that for a small peptide, deca-alanine, metastable states can be identified through a neural net based on structural information alone. While processing molecular dynamics data, dimensionality reduction is a necessary step that projects high-dimensional data onto a low-dimensional representation that, ideally, captures the conformational changes in the underlying data. Conventional methods make use of the temporal information contained in trajectories generated through integrating the equations of motion, which forgoes more efficient sampling schemes. We demonstrate that EncoderMap, an autoencoder architecture with an additional distance metric, can find a suitable low-dimensional representation to identify long-lived molecular conformations using exclusively structural information. For deca-alanine, which exhibits several helix-forming pathways, we show that this approach allows us to combine simulations with different biasing forces and yields representations comparable in quality to other established methods. Our results contribute to computational strategies for the rapid automatic exploration of the configuration space of peptides and proteins. [ABSTRACT FROM AUTHOR]

When an activated complex, as defined in transition state theory (TST), has a polyhedral shape, its kinetic energy is found to be diagonal in a system of spherical polar coordinates. If, in addition, the polyhedron is characterized by a high symmetry, then its dynamics considerably simplifies. An application of this approach to the most symmetrical TS known to date, i.e., that which controls the Cl− + CH3Cl → ClCH3 + Cl− SN2 nucleophilic substitution, is presented and an analytical expression of its potential energy surface is provided. In a substantial range around the saddle point, approximate equations of motion for the two components of the reaction coordinate, i.e., the antisymmetrical stretching motion of the ClCCl core and the wagging motion of the hydrogen triad, can be derived in an analytical form. During an extensive period of time, the main component of the reaction coordinate is governed by an unexpectedly simple equation of motion that depends on a single initial condition, irrespective of the other ones and of the internal energy. Reactive trajectories are observed to form a perfectly collimated bundle characterized by undetectable dispersion, thereby giving a spectacular example of regular dynamics in an anharmonic potential. Regularity and collimation are brought about by local symmetry, which is a widespread feature of potential energy surfaces. Anharmonicity is observed to influence the dynamics only at a late stage. As energy increases, trajectories tend to fan out and to deviate from the analytical equation. For the wagging motion, chaos sets in at much lower energies. [ABSTRACT FROM AUTHOR]

Numerically "exact" methods addressing the dynamics of coupled electron–phonon systems have been intensively developed. Nevertheless, the corresponding results for the electron mobility μdc are scarce, even for the one-dimensional (1d) Holstein model. Building on our recent progress on single-particle properties, here we develop the momentum-space hierarchical equations of motion (HEOM) method to evaluate real-time two-particle correlation functions of the 1d Holstein model at a finite temperature. We compute numerically "exact" dynamics of the current–current correlation function up to real times sufficiently long to capture the electron's diffusive motion and provide reliable results for μdc in a wide range of model parameters. In contrast to the smooth ballistic-to-diffusive crossover in the weak-coupling regime, we observe a temporally limited slow-down of the electron on intermediate time scales already in the intermediate-coupling regime, which translates to a finite-frequency peak in the optical response. Our momentum-space formulation lowers the numerical effort with respect to existing HEOM-method implementations, while we remove the numerical instabilities inherent to the undamped-mode HEOM by devising an appropriate hierarchy closing scheme. Still, our HEOM remains unstable at too low temperatures, for too strong electron–phonon coupling, and for too fast phonons. [ABSTRACT FROM AUTHOR]

We present a formulation of special relativistic dissipative hydrodynamics (SRDHD) derived from the well-established Müller–Israel–Stewart (MIS) formalism using an expansion in deviations from ideal behaviour. By re-summing the non-ideal terms, our approach extends the Euler equations of motion for an ideal fluid through a series of additional source terms that capture the effects of bulk viscosity, shear viscosity, and heat flux. For efficiency these additional terms are built from purely spatial derivatives of the primitive fluid variables. The series expansion is parametrized by the dissipation strength and time-scale coefficients, and is therefore rapidly convergent near the ideal limit. We show, using numerical simulations, that our model reproduces the dissipative fluid behaviour of other formulations. As our formulation is designed to avoid the numerical stiffness issues that arise in the traditional MIS formalism for fast relaxation time-scales, it is roughly an order of magnitude faster than standard methods near the ideal limit. [ABSTRACT FROM AUTHOR]

We introduce the fused sequential addition and migration (fSAM) algorithm for generating microstructures of fiber composites with long, flexible, nonoverlapping fibers and industrial volume fractions. The proposed algorithm is based on modeling the fibers as polygonal chains and enforcing, on the one hand, the nonoverlapping constraints by an optimization framework. The connectivity constraints, on the other hand, are treated via constrained mechanical systems of d'Alembert type. In case of straight, that is, nonflexible, fibers, the proposed algorithm reduces to the SAM (Comput. Mech., 59, 247–263, 2017) algorithm, a well‐established method for generating short fiber‐reinforced composites. We provide a detailed discussion of the equations governing the motion of a flexible fiber and discuss the efficient numerical treatment. We elaborate on the integration into an existing SAM code and explain the selection of the numerical parameters. To capture the fiber length distributions of long fiber reinforced composites, we sample the fiber lengths from the Gamma distribution and introduce a strategy to incorporate extremely long fibers. We study the microstructure generation capabilities of the proposed algorithm. The computational examples demonstrate the superiority of the novel microstructure‐generation technology over the state of the art, realizing large fiber aspect ratios (up to 2800) and high fiber volume fractions (up to 32%$$ 32\% $$ for an aspect ratio of 150) for experimentally measured fiber orientation tensors. [ABSTRACT FROM AUTHOR]

We take a fresh look at the relation between generalised Kähler geometry and N = (2 , 2) supersymmetric sigma models in two dimensions formulated in terms of (2, 2) superfields. Dual formulations in terms of different kinds of superfield are combined to give a formulation with a doubled target space and both the original superfield and the dual superfield. For Kähler geometry, we show that this doubled geometry is Donaldson's deformation of the holomorphic cotangent bundle of the original Kähler manifold. This doubled formulation gives an elegant geometric reformulation of the equations of motion. We interpret the equations of motion as the intersection of two Lagrangian submanifolds (or of a Lagrangian submanifold with an isotropic one) in the infinite dimensional symplectic supermanifold which is the analogue of phase space. We then consider further extensions of this formalism, including one in which the geometry is quadrupled, and discuss their geometry. [ABSTRACT FROM AUTHOR]

This study analyzes the Pazarcik and Elbistan earthquakes, which occurred on February 6, 2023 and are among the most destructive seismic events of the 21st century. Since the greatest damage was seen in Hatay in these earthquakes centered in Kahramanmaras, the study aims to contribute to the field of earthquake engineering by evaluating the seismic data obtained from these regions. In the first part of the analysis, peak ground accelerations (PGA) recorded at the stations in Kahramanmaras and Hatay were examined and these data were compared with the DD1 (maximum considered earthquake) and DD2 (design basis earthquake) design levels of the Turkish Building Earthquake Code (TBEC 2018). In addition, the effects of ground properties and proximity of faults on seismic records obtained from various stations were evaluated by examining the PGA distribution spatially. The impacts of factors such as the near-fault effect, directivity effect, ground amplification effect and impulse-like motions were determined by examining the peak ground accelations, peak ground velocities, peak ground displacements and spectral accelerations. The study uses NGA-West2 Ground Motion Prediction Equations (GMPEs) to evaluate peak ground accelerations in stiff soil and to consider impulse and directivity effects. In addition, the applicability of USGS Vs30 maps in Turkiye is evaluated by comparing with AFAD data. These comprehensive analysis provide critical insights from a structural safety perspective on how seismic characteristics change with ground properties and proximity to earthquake sources. [ABSTRACT FROM AUTHOR]

Recent (sub)millimetre polarization observations of protoplanetary discs reveal toroidally aligned, effectively prolate dust grains large enough (at least |$\sim 100$| |$\mu$| m) to efficiently scatter millimetre light. The alignment mechanism for these grains remains unclear. We explore the possibility that gas drag aligns grains through gas–dust relative motion when the grain's centre of mass is offset from its geometric centre, analogous to a badminton birdie's alignment in flight. A simple grain model of two non-identical spheres illustrates how a grain undergoes damped oscillations from flow-induced restoring torques which align its geometric centre in the flow direction relative to its centre of mass. Assuming specular reflection and subsonic flow, we derive an analytical equation of motion for spheroids where the centre of mass can be shifted away from the spheroid's geometric centre. We show that a prolate or an oblate grain can be aligned with the long axis parallel to the gas flow when the centre of mass is shifted along that axis. Both scenarios can explain the required effectively prolate grains inferred from observations. Application to a simple disc model shows that the alignment time-scales are shorter than or comparable to the orbital time. The grain alignment direction in a disc depends on the disc (sub-)structure and grain Stokes number (St) with azimuthal alignment for large St grains in sub-Keplerian smooth gas discs and for small St grains near the gas pressure extrema, such as rings and gaps. [ABSTRACT FROM AUTHOR]

The present study investigates the vibration analysis of cylindrical shells composed of fiber metal laminate (FML) with embedded piezoelectric layers, undergoing fluid-structure interaction (FSI) and resting on a Pasternak elastic foundation based on the principles of three-dimensional elasticity theory. Using the state space approach, the equations of motion were derived under simply supported boundary conditions. The natural frequencies of the FML cylindrical shell, accounting for the presence of a moving fluid, were computed by solving the eigenfrequency equations. The study examined the influence of various parameters, including boundary conditions, length-to-radius ratio, fluid type, fluid velocity, circumferential wave number, and radius-to-thickness ratio, on glass-reinforced aluminum laminate (GLARE), aramid-reinforced aluminum laminate (ARALL), and carbon-reinforced aluminum laminate (CARALL). A constant composite/metal volume ratio was assumed. The results obtained were validated by comparing with natural frequency values from the existing literature, confirming the agreement and convergence with previous studies. The results confirm that the highest natural frequency values are assigned to the CARALL, ARALL and GRALE structures in descending order. Furthermore, an increase in fluid flow velocity through the cylindrical shell correlates with a reduction in natural frequency. [ABSTRACT FROM AUTHOR]

This paper investigates the wave propagation in the graphene platelets (GPLs)-enhanced functionally graded porous plates. The governing equations of motion are obtained using the first-order shear deformation plate theory (FSDT). Subsequently, the equations are transformed into the state-space form. The wave dispersion relation is derived by solving the state space equation by means of the method of reverberation-ray matrix and the accuracy of this approach is validated through a comparative analysis with the results from relevant literature. In addition, parametric analyses are carried out, including boundary conditions, porosity coefficient, GPL mass fraction, porosity distribution, GPL distribution, and thickness-to-width ratio, on the dispersion behavior of functionally graded GPLs-reinforced porous plates. The use of GPLs in these composites is particularly promising, and the findings offer valuable insights into the design of composites with tailored properties for specific engineering applications. [ABSTRACT FROM AUTHOR]

The stability of the parametric resonance of a controlled stay cable with time delay is investigated. The in-plane nonlinear equations of motion are initially determined via the Hamilton principle. Then, utilizing the method of multiple scales, the modulation equations that govern the nonlinear dynamics are obtained. These equations are then utilized to investigate the effect of time delays on the amplitude and frequency-response behavior and, subsequently, on the stability of the parametric resonance of the controlled cable, that it is shown to depend on the excitation amplitude and the commensurability of the delayed-response frequency to the excitation frequency. The stability region of the parametric resonance is shifted, and the effects of control on the cable become worse by increasing time delay. The work plays a guiding role in the parametric design of the control system for stay cables. [ABSTRACT FROM AUTHOR]

Improving the performance of nutation blood pumps has long been a focus of medical research. In this study, we introduce the finite element method combined with sensor-based experimental measurements to explore the factors affecting the stability of the levitated magnetic rotor in the blood pump. Initially, we derived Euler angle and motion state equations to describe the motion of the magnetically levitated rotor inside the blood pump. Subsequently, employing the finite element method, we analyzed various factors impacting the dynamic levitation stability of the rotor. Our results indicate that the size of the permanent magnet in the rotor bearing, the nutation angle, and the levitation gap are critical factors affecting the rotor's levitation stability. We discuss these crucial parameters in detail, highlighting their correlation with the levitation stability of the rotor. Furthermore, comparing our theoretically obtained magnetic forces with experimental results, we found the maximum difference to be within 8%. Overall, this study provides valuable insights to enhance the functionality of blood pumps and understand the key parameters affecting their design. [ABSTRACT FROM AUTHOR]

In this work, a rectangular cellular microplate is taken into consideration which is embedded between two functionally graded carbon nanotube-reinforced composite layers that have the piezoelectric ability. Different patterns for the carbon nanotube dispersion are considered. Moreover, an external electrical voltage is applied to them. The displacements are described based on a higher-order shear deformation theory which is called hyperbolic theory and the size influence is captured via the modified couple stress formulations. The governing motion equations are then derived using the variational technique and Hamilton’s principle. Then, for simply supported edge conditions, a closed-form solution is provided. Next, it turns to validate the results in simpler states, and after that, the effect of the most prominent parameters on the results is investigated. It is observed that by increasing the external applied voltage to the face sheets, the frequencies are reduced. Also, the natural frequencies have a tendency to decrease as the dimensionless material length scale parameter increases. This study’s outcomes may help design and manufacture micro-/nano-electro systems, sensors and actuators, small-scale devices, and other engineering structures. [ABSTRACT FROM AUTHOR]