Quantum time correlation functions (TCFs) involving two states are important for describing nonadiabatic dynamical processes such as charge transfer (CT). Based on a previous single-state method, we propose an imaginary-time open-chain path-integral (OCPI) approach for evaluating the two-state symmetrized TCFs. Expressing the forward and backward propagation on different electronic potential energy surfaces as a complex-time path integral, we then transform the path variables to average and difference variables such that the integration over the difference variables up to the second order can be performed analytically. The resulting expression for the symmetrized TCF is equivalent to sampling the open-chain configurations in an effective potential that corresponds to the average surface. Using importance sampling over the extended OCPI space via open path-integral molecular dynamics, we tested the resulting path-integral approximation by calculating the Fermi’s golden rule CT rate constant within a widely used spin-boson model. Comparing with the real-time linearized semiclassical method and analytical result, we show that the imaginary-time OCPI provides an accurate two-state symmetrized TCF and rate constant in the typical turnover region. It is shown that the first bead of the open chain corresponds to physical zero-time and that the endpoint bead corresponds to final time t; oscillations of the end-to-end distance perfectly match the nuclear mode frequency. The two-state OCPI scheme is seen to capture the tested model’s electronic quantum coherence and nuclear quantum effects accurately.
more »
« less
This content will become publicly available on August 1, 2024
Matrix Approach to Accelerate Spin‐Up of CLM5
Key Points A new semi‐analytical spin‐up (SASU) framework combines the default accelerated spin‐up method and matrix analytical algorithm SASU accelerates CLIM5 spin‐up by tens of times, becoming the fastest method to our knowledge SASU is applicable to most biogeochemical models and enables computationally costly study, for example, sensitivity analysis
more »
« less
- NSF-PAR ID:
- 10446677
- Date Published:
- Journal Name:
- Journal of Advances in Modeling Earth Systems
- Volume:
- 15
- Issue:
- 8
- ISSN:
- 1942-2466
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Abstract Exosomes, a subset of extracellular vesicles (EVs, 30–200‐nm diameter), serve as biomolecular snapshots of their cell of origin and vehicles for intercellular communication, playing roles in biological processes, including homeostasis maintenance and immune modulation. The large‐scale processing of exosomes for use as therapeutic vectors has been proposed, but these applications are limited by impure, low‐yield recoveries from cell culture milieu (CCM). Current isolation methods are also limited by tedious and laborious workflows, especially toward an isolation of EVs from CCM for therapeutic applications. Employed is a rapid (<10 min) EV isolation method on a capillary‐channeled polymer fiber spin‐down tip format. EVs are isolated from the CCM of suspension‐adapted human embryonic kidney cells (HEK293), one of the candidate cell lines for commercial EV production. This batch solid‐phase extraction technique allows 1012EVs to be obtained from only 100‐µl aliquots of milieu, processed using a benchtop centrifuge. The tip‐isolated EVs were characterized using transmission electron microscopy, multi‐angle light scattering, absorbance quantification, an enzyme‐linked immunosorbent assay to tetraspanin marker proteins, and a protein purity assay. It is believed that the demonstrated approach has immediate relevance in research and analytical laboratories, with opportunities for production‐level scale‐up projected.
-
more » « less
Abstract Numerical techniques to efficiently model out-of-equilibrium dynamics in interacting quantum many-body systems are key for advancing our capability to harness and understand complex quantum matter. Here we propose a new numerical approach which we refer to as generalized discrete truncated Wigner approximation (GDTWA). It is based on a discrete semi-classical phase space sampling and allows to investigate quantum dynamics in lattice spin systems with arbitrary
S ≥ 1/2. We show that the GDTWA can accurately simulate dynamics of large ensembles in arbitrary dimensions. We apply it forS > 1/2 spin-models with dipolar long-range interactions, a scenario arising in recent experiments with magnetic atoms. We show that the method can capture beyond mean-field effects, not only at short times, but it also can correctly reproduce long time quantum-thermalization dynamics. We benchmark the method with exact diagonalization in small systems, with perturbation theory for short times, and with analytical predictions made for models which feature quantum-thermalization at long times. We apply our method to study dynamics in largeS > 1/2 spin-models and compute experimentally accessible observables such as Zeeman level populations, contrast of spin coherence, spin squeezing, and entanglement quantified by single-spin Renyi entropies. We reveal that largeS systems can feature larger entanglement than correspondingS = 1/2 systems. Our analyses demonstrate that the GDTWA can be a powerful tool for modeling complex spin dynamics in regimes where other state-of-the art numerical methods fail. -
Dimerized valence bond solids appear naturally in spin-1/2 systems on bipartite lattices, with the geometric frustrations playing a key role both in their stability and the eventual `melting' due to quantum fluctuations. Here, we ask the question of the stability of such dimerized solids in spin-1 systems, taking the anisotropic square lattice with bilinear and biquadratic spin-spin interactions as a paradigmatic model. The lattice can be viewed as a set of coupled spin-1 chains, which in the limit of vanishing inter-chain coupling are known to possess a stable dimer phase. We study this model using the density matrix renormalization group (DMRG) and infinite projected entangled-pair states (iPEPS) techniques, supplemented by the analytical mean-field and linear flavor wave theory calculations. While the latter predicts the dimer phase to remain stable up to a reasonably large interchain-to-intrachain coupling ratio r≲0.6, the DMRG and iPEPS find that the dimer solid melts for much weaker interchain coupling not exceeding r≲0.15. We find the transition into a magnetically ordered state to be first order, manifested by a hysteresis and order parameter jump, precluding the deconfined quantum critical scenario. The apparent lack of stability of dimerized phases in 2D spin-1 systems is indicative of strong quantum fluctuations that melt the dimer solid.more » « less
-
An analytical implementation of static dipole polarizabilities within the generalized Kohn–Sham semicanonical projected random phase approximation (GKS-spRPA) method for spin-restricted closed-shell and spin-unrestricted open-shell references is presented. General second-order analytical derivatives of the GKS-spRPA energy functional are derived using a Lagrangian approach. By resolution-of-the-identity and complex frequency integration methods, an asymptotic [Formula: see text] scaling of operation count and [Formula: see text] scaling of storage is realized, i.e., the computational requirements are comparable to those for GKS-spRPA ground state energies. GKS-spRPA polarizabilities are assessed for small molecules, conjugated long-chain hydrocarbons, metallocenes, and metal clusters, by comparison against Hartree–Fock (HF), semilocal density functional approximations (DFAs), second-order Møller–Plesset perturbation theory, range-separated hybrids, and experimental data. For conjugated polydiacetylene and polybutatriene oligomers, GKS-spRPA effectively addresses the “overpolarization” problem of semilocal DFAs and the somewhat erratic behavior of post-PBE RPA polarizabilities without empirical adjustments. The ensemble averaged GKS-spRPA polarizabilities of sodium clusters (Na n for n = 2, 3, …, 10) exhibit a mean absolute deviation comparable to PBE with significantly fewer outliers than HF. In conclusion, analytical second-order derivatives of GKS-spRPA energies provide a computationally viable and consistent approach to molecular polarizabilities, including systems prohibitive for other methods due to their size and/or electronic structure.more » « less
