Search NASA⌕ Search

SEARCH · Search NASA

Results for “Full coupled”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 379 records · Page 21

Extension of the high-resolution thermal-hydraulics code ESCOT to hexagonal core geometries for multi-physics calculations

The extension of the capabilities of the pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT to analyze hexagonal fueled cores and its performance are presented. ESCOT is an accurate yet fast core thermal-hydraulics solution aiming at high-fidelity and high-resolution multi-physics core analysis in the framework of massively parallel computing platforms. Its algorithm solution is based on the four-equation drift-flux model for two-phase calculations, these are numerically solved by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models such as turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analysis. ESCOT is parallelized by a double (radial and axial) domain decomposition that enables its highly parallelized execution. The coupling of the code with the neutronics whole core solver for hexagonal geometries nTRACER is described. The newly implemented ESCOT features are validated by comparing single assembly and full core steady state nTRACER-ESCOT solutions with nTRACER standalone internal one-dimensional T/H solver results. The validation problems are based on the VVER 440 and VVER 1000 cores. ESCOT results show differences within an acceptable range with respect to the simple 1D nTRACER built-in solver. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

The new physics case for beam-dump experiments with accelerated muon beams

As the field examines a future muon collider as a possible successor to the LHC, we must consider how to fully utilize not only the high-energy particle collisions, but also any lower-energy staging facilities necessary in the R&D process. An economical and efficient possibility is to use the accelerated muon beam from either the full experiment or from cooling and acceleration tests in beam-dump experiments. Beam-dump experiments are complementary to the main collider as they achieve sensitivity to very small couplings with minimal instrumentation. We demonstrate the utility of muon beam-dump experiments for new physics searches at energies from 10 GeV to 5 TeV. We find that, even at low energies like those accessible at staging or demonstrator facilities, it is possible to probe new regions of parameter space for a variety of generic BSM models, including muonphilic, leptophilic, L μ − L τ , and dark photon scenarios. Such experiments could therefore provide opportunities for discovery of new physics well before the completion of the full multi-TeV collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Warm Arctic, Cold Siberia Pattern: Role of Full Arctic Amplification Versus Sea Ice Loss Alone

The effect of future Arctic amplification (AA) on the extratropical atmospheric circulation remains unclear in modeling studies. Using a collection of coordinated atmospheric and coupled global climate model perturbation experiments, we find an emergent relationship between the high-latitude 1,000–500 hPa thickness response and an enhancement of the Siberian High in winter. This wave number-1-like sea level pressure anomaly pattern is linked to an equatorward shift of the eddy-driven jet and a dynamical cooling response in eastern Asia. Additional simulations, where AA is imposed directly into the model domain by nudging, demonstrate how the sea ice forcing is insufficient by itself to capture the vertical extent of the warming and by extension the amplitude of the response in the Siberian High. Additionally, this study demonstrates the importance of the vertical extent of the tropospheric warming over the polar cap in revealing the “warm Arctic, cold Siberia” anomaly pattern in future projections.

54 ENVIRONMENTAL SCIENCES↗

Spin-phonon dispersion in magnetic materials

Microscopic coupling between the electron spin and the lattice vibration is responsible for an array of exotic properties from morphic effects in simple non magnets to magnetodielectric coupling in multiferroic spinels and hematites. Traditionally, a single spin–phonon coupling constant is used to characterize how effectively the lattice can affect the spin, but it is hardly enough to capture novel electromagnetic behaviors to the full extent. Here, we introduce a concept of spin–phonon dispersion to project the spin moment change along the phonon crystal momentum direction, so the entire spin change can be mapped out. Different from the phonon dispersion, the spin–phonon dispersion has both positive and negative frequency branches even in the equilibrium ground state, which correspond to the spin enhancement and spin reduction, respectively. Our study of bcc Fe and hcp Co reveals that the spin force matrix, that is, the second-order spatial derivative of spin moment, is similar to the vibrational force matrix, but its diagonal elements are smaller than the off-diagonal ones. This leads to the distinctive spin–phonon dispersion. Furthermore, the concept of spin–phonon dispersion expands the traditional Elliott-Yafet theory in nonmagnetic materials to the entire Brillouin zone in magnetic materials, thus opening the door to excited states in systems such as CoF 2 and NiO, where a strong spin-lattice coupling is detected in the THz regime.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Preliminary design overview and performance assessment of the ITER low-field side reflectometer

The design of the ITER low-field side reflectometer (LFSR) has matured to a complete end-to-end preliminary design. LFSR will supply three key plasma measurements: (1) electron density profile, (2) electron density fluctuations, and (3) poloidal rotation. Simultaneous measurements of the three quantities are enabled by an array of six monostatic antennas which inject from an equatorial port on the outboard side of the ITER vessel. Low-loss transmission lines, consisting of corrugated, overmoded waveguide and miter bends, transmit the 30–165 GHz, O- and X-mode signals to and from the ITER plasma. Integrated transmission-line components serve a range of purposes, such as protection from high-power stray radiofrequency radiation, accommodation of transmission-line displacement, and simultaneous measurement of reference and plasma phases during the discharge. Broadband transmission signals are realized by full-band microwave transceivers combined with quasi-optical multiplexing. A field-programmable gate array (FPGA) processor demodulates the profile reflectometer signals, enabling real-time density profile measurements for plasma control system feedback. A full-scale transmission line test facility provides an integrated environment to assess the performance of critical LFSR components. Theoretical modeling together with insertion loss measurements provide the basis for a comprehensive power budget, which accounts for transmitted output power, transmission-line losses, antenna coupling, and plasma effects. Results indicate that high signal-to-noise ratios are achievable with the current design. A synthetic reflectometer model, using real design parameters and baseline ITER profiles, has been developed to estimate the return signal. With evolving microwave and data acquisition technologies, full-band, ultrafast sweeps (<1 μs) will be realizable for ITER.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

An Anthradithiophene Donor Polymer for Organic Solar Cells with a Good Balance between Efficiency and Synthetic Accessibility

Over the past few years, the organic photovoltaics technology (OPV) has reached remarkable power conversion efficiencies (PCEs), mostly thanks to the advent of nonfullerene acceptors as well as to a high level of materials engineering. Yet, the complex materials syntheses behind these results ultimately limit technological readiness. The quest for scalable organic compounds offering high PCE and reasonably low synthetic complexity is a must to close the gap between laboratory R&D and commercial products. The synthesis and full characterization of a new conjugated polymer called PATTD is reported herein, based on a novel anthradithiophene as an electron–rich building block coupled with a commercially available dithienylbenzodithiophenedione as an electron–withdrawing comonomer. Its photovoltaic properties are studied in blends with IT–4F and IDIC as acceptors. PATTD:IT–4F–based photovoltaic devices exhibit a PCE approaching 10% and over 2300 h shelf–life stability. Finally, the PATTD scalability factor (SF), together with the PATTD–based photovoltaic performances, lead to a PCE/SF value equal to 0.297, placing such devices into the innermost circle of OPV materials, achieving one of the best compromises between efficiency and synthetic complexity and at the same time offering promising industrial perspectives.

14 SOLAR ENERGY↗

Visualization of conformational changes and membrane remodeling leading to genome delivery by viral class-II fusion machinery

Chikungunya virus (CHIKV) is a human pathogen that delivers its genome to the host cell cytoplasm through endocytic low pH-activated membrane fusion mediated by class-II fusion proteins. Though structures of prefusion, icosahedral CHIKV are available, structural characterization of virion interaction with membranes has been limited. Here, we have used cryo-electron tomography to visualize CHIKV’s complete membrane fusion pathway, identifying key intermediary glycoprotein conformations coupled to membrane remodeling events. Using sub-tomogram averaging, we elucidate features of the low pH-exposed virion, nucleocapsid and full-length E1-glycoprotein’s post-fusion structure. Contrary to class-I fusion systems, CHIKV achieves membrane apposition by protrusion of extended E1-glycoprotein homotrimers into the target membrane. The fusion process also features a large hemifusion diaphragm that transitions to a wide pore for intact nucleocapsid delivery. Our analyses provide comprehensive ultrastructural insights into the class-II virus fusion system function and direct mechanistic characterization of the fundamental process of protein-mediated membrane fusion.

59 BASIC BIOLOGICAL SCIENCES↗

Topological excitonic insulator with tunable momentum order

Correlated topological materials often maintain a delicate balance among physical symmetries. Many topological orders are symmetry protected, whereas most correlated phenomena arise from spontaneous symmetry breaking. Cases where symmetry breaking induces a non-trivial topological phase are rare. Here we demonstrate the presence of two such phases in Ta 2 Pd 3 Te 5 , where Coulomb interactions form excitons that condense below 100 K, one with zero and the other with finite momentum. We observed a full spectral bulk gap, which stems from exciton condensation. This topological excitonic insulator state spontaneously breaks mirror symmetries but involves a weak structural coupling. Scanning tunnelling microscopy shows gapless boundary modes in the bulk insulating phase. Their magnetic field response, together with theoretical modelling, indicates a topological origin. These observations establish Ta 2 Pd 3 Te 5 as a topological excitonic insulator in a three-dimensional crystal. Finally, our results manifest a unique sequence of topological exciton condensations in a bulk crystal, offering exciting opportunities to study critical behaviour and excitations.

36 MATERIALS SCIENCE↗

Rapidly convergent quantum Monte Carlo using a Chebyshev projector

The multireference coupled-cluster Monte Carlo (MR-CCMC) algorithm is a determinant-based quantum Monte Carlo (QMC) algorithm that is conceptually similar to Full Configuration Interaction QMC (FCIQMC). It has been shown to offer a balanced treatment of both static and dynamic correlation while retaining polynomial scaling, although application to large systems with significant strong correlation remained impractical. In this paper, we document recent algorithmic advances that enable rapid convergence and a more black-box approach to the multireference problem. These include a logarithmically scaling metric-tree-based excitation acceptance algorithm to search for determinants connected to the reference space at the desired excitation level and a symmetry-screening procedure for the reference space. We show that, for moderately sized reference spaces, the new search algorithm brings about an approximately 8-fold acceleration of one MR-CCMC iteration, while the symmetry screening procedure reduces the number of active reference space determinants with essentially no loss of accuracy. We also introduce a stochastic implementation of an approximate wall projector, which is the infinite imaginary time limit of the exponential projector, using a truncated expansion of the wall function in Chebyshev polynomials. Notably, this wall-Chebyshev projector can be used to accelerate any projector-based QMC algorithm. We show that it requires significantly fewer applications of the Hamiltonian to achieve the same statistical convergence. We benchmark these acceleration methods on the beryllium and carbon dimers, using initiator FCIQMC and MR-CCMC with basis sets up to cc-pVQZ quality.

Zhao, Zijun↗

Line Emission Mapper: an X-ray probe mission concept to study the cosmic ecosystems and the physics of galaxy formation

In the 2020 Astrophysics Decadal Survey, the National Academies identified cosmic feedback and structure formation as a key question that should drive research in the upcoming decade. In response to this recommendation, NASA released a call for X-ray and IR probe-class missions, with a $\$$1B cost cap. The line emission mapper (LEM) is a mission concept designed in response to this call. LEM is a single-instrument X-ray telescope that consists of a Wolter–Schwarzschild type I X-ray optic with a 4 m focal length, coupled with an X-ray microcalorimeter with a 30' field of view (FoV), 15" angular resolution, and 2.5 eV energy resolution [full-width half maximum (FWHM)], with a 1.3 eV FWHM energy resolution central subarray. The high throughput X-ray mirror combined with the large FoV and excellent energy resolution allows for efficient mapping of extended emission-line dominated astrophysical objects from megaparsecs to sub-pc scales to study cosmic ecosystems and unveil the physical drivers of galaxy formation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

$\mathrm{O}(a)$ improvement of the flavour singlet scalar density in a setup with Wilson fermions

We report on our Ward identity determination of the O(a) improvement coefficient for the flavour singlet scalar density, namely gS , from three-flavour lattice QCD with Wilson-clover fermions and the tree-level Symanzik improved gauge action. We employ five couplings, g20∈[1.5,1.77] , that cover the range used in large-volume CLS simulations. While gS itself is for instance relevant for the O(a) improvement of meson and baryon sigma terms, a relation to bg , the O(a) improvement parameter of the gauge coupling, can also be established, allowing for its non-perturbative extraction as well. With Wilson fermions, bg is in principle required for full O(a) improvement at non-vanishing sea quark masses. We outline our procedure for extracting bg

Petrak, Pia Jones↗

Ch3MS-RF: a random forest model for chemical characterization and improved quantification of unidentified atmospheric organics detected by chromatography–mass spectrometry techniques

Abstract. The chemical composition of ambient organic aerosols plays a critical role in driving their climate and health-relevant properties and holds important clues to the sources and formation mechanisms of secondary aerosol material. In most ambient atmospheric environments, this composition remains incompletely characterized, with the number of identifiable species consistently outnumbered by those that have no mass spectral matches in the literature or the National Institute of Standards and Technology/National Institutes of Health/Environmental Protection Agency (NIST/NIH/EPA) mass spectral databases, making them nearly impossible to definitively identify. This creates significant challenges in utilizing the full analytical capabilities of techniques which separate and generate spectra for complex environmental samples. In this work, we develop the use of machine learning techniques to quantify and characterize novel, or unidentifiable, organic material. This work introduces Ch3MS-RF (Chemical Characterization by Chromatography–Mass Spectrometry Random Forest Modeling), an open-source, R-based software tool, for efficient machine-learning-enabled characterization of compounds separated in chromatography–mass spectrometry applications but not identifiable by comparison to mass spectral databases. A random forest model is trained and tested on a known 130 component representative external standard to predict the response factors of novel environmental organics based on position in volatility–polarity space and mass spectrum, enabling the reproducible, efficient, and optimized quantification of novel environmental species. Quantification accuracy on a reserved 20 % test set randomly split from the external standard compound list indicates that random forest modeling significantly outperforms the commonly used methods in both precision and accuracy, with a median response factor percent error of −2 %, for modeled response factors, compared to > 15 %, for typically used proxy assignment-based methods. Chemical properties modeling, evaluated on the same reserved 20 % test set and an extrapolation set of species identified in ambient organic aerosol samples collected in the Amazon rainforest, also demonstrate robust performance. Extrapolation set property prediction mean absolute errors for carbon number, oxygen to carbon ratio (O : C), average carbon oxidation state (OSc‾), and vapor pressure are 1.8, 0.15, 0.25, and 1.0 (log(atm)), respectively. Extrapolation set out-of-sample R2 for all properties modeled are above 0.75, with the exception of vapor pressure. While predictive performance for vapor pressure is less robust compared to the other chemical properties modeled, random-forest-based modeling was significantly more accurate than other commonly used methods of vapor pressure prediction, decreasing the mean vapor pressure prediction error to 0.24 (log(atm)) from 0.55 (log(atm)) (chromatography-based vapor pressure prediction) and 1.2 (log(atm)) (chemical formula-based vapor pressure prediction). The random forest model significantly advances an untargeted analysis of the full scope of chemical speciation yielded by two-dimensional gas chromatography (GCxGC-MS) techniques and can be applied to gas chromatography coupled with electron ionization mass spectrometry (GC-MS) as well. It enables the accurate estimation of key chemical properties commonly utilized in the atmospheric chemistry community, which may be used to more efficiently identify important tracers for further individual analysis and to characterize compound populations uniquely formed under specific ambient conditions.

54 ENVIRONMENTAL SCIENCES↗

Multiphysics Simulations of MSRE with NEAMS Thermal Hydraulics Tools

This report documents the benchmarks being developed and simulations performed using tools and codes developed under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, utilizing MSRE experimental data. In FY23, three main work scopes were investigated under the NEAMS MSR work package at ANL. The first scope investigated the Griffin-SAM coupling model for simulating the pump startup transient experiment of MSRE. The analyses start with a simple model (single-channel, single-lattice), gradually adding more details (multi-channel, full-core) into the model. The results show that the reactivity loss curve is very sensitive to the axial boundary conditions and the radial core discretization. The simple model can predict a similar reactivity trend as that of the more sophisticated model, which is likely due to error cancellation. Accurately modeling the axial boundary condition may further improve the reactivity trend but would require significant efforts to generate the mesh of the MSRE inlet and upper plenum. The core channel radial discretization for the Griffin-SAM coupled model also depends on the flow distribution. Given the complex geometry in the inlet plenum, the flow distribution needed to be calculated from CFD analysis, which was performed using the NekRS code. This analysis employed a MSRE CAD model developed by Copenhagen Atomics. The CAD model was disassembled to keep the inlet plenum region only, which was subsequently cleaned and modified so that the mesh generated is under the memory limit. The results are merged to a few radial regions to show that the flow rate is highest in the central region. This would be useful for future improvement of the Griffin-SAM coupling model of the MSRE core. The last task investigated is tritium transport modeling using the standalone SAM code. This task aimed to initiate the effort to demonstrate and validate the tritium transport model implemented in SAM. The preliminary investigation employed an MSRE model consisting of the primary loop. Three tritium transport pathways were examined including the retention in the graphite, the permeation through the HX tube wall, and the removal from the off-gas system. The results compare well with the MSRE data, but improvements are still needed on the initial conditions (i.e., the present state may not have reached equilibrium), the boundary conditions, the off-gas system modeling, and a better numerical strategy to reach the equilibrium state.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Design and Engineering of LUPIN: A Test-Bed Radio-Frequency Ion Source for Enhanced Neutral Beam Injection on DIII-D

The Large, Uniform Plasma for Ionizing Neutrals (LUPIN) is a radio-frequency (RF) inductively coupled plasma (ICP) chamber for demonstrating plasma performance of an RF ICP positive ion source upgrade for the DIII-D neutral beam injection (NBI) system. LUPIN will be used to investigate ion source physics, including neutral gas dynamics, plasma density uniformity, interactions with Faraday shields, and power coupling to novel RF antenna designs. LUPIN has an RF generator capable of delivering 20 kW of power at 2 MHz, which is coupled into a cylindrical quartz vessel measuring 20 cm in length and 10 cm in radius. This configuration matches the power density requirements for a full-scale ion source. Target hydrogen and deuterium plasma densities exceeding 10 18 m -3 would relate to extracted ion current densities of 2100 A/m 2 for 10s. Vacuum conductance and gas flow calculations predict a maximum achievable neutral gas flow rate of 1675 Pa ⋅ L/s at 5 Pa of He, which mimics the gas flow of the DIII-D NBI system. Designs have been developed for an internal Faraday shield to mitigate heat flux and ion sputtering on the dielectric vessel. Thermomechanical finite element simulations demonstrated the Faraday shield design to be capable of withstanding anticipated heat loads from worst-case operation scenarios. Finally, results of upcoming experimental investigations on LUPIN will guide the design of a full-scale prototype for DIII-D integration.

Faraday shield↗

Uncertainty-informed selection of CMIP6 Earth System Model subsets for use in multisectoral and impact models

Earth system models (ESMs) and general circulation models (GCMs) are heavily used to provide inputs to sectoral impact and multisector dynamic models, which include representations of energy, water, land, economics, and their interactions. Therefore, representing the full range of model uncertainty, scenario uncertainty, and interannual variability that ensembles of these models capture is critical to the exploration of the future co-evolution of the integrated human–Earth system. The pre-eminent source of these ensembles has been the Coupled Model Intercomparison Project (CMIP). With more modeling centers participating in each new CMIP phase, the size of the model archive is rapidly increasing, which can be intractable for impact modelers to effectively utilize due to computational constraints and the challenges of analyzing large datasets. In this work, we present a method to select a subset of the latest phase, CMIP6, featuring models for use as inputs to a sectoral impact or multisector dynamics models, while prioritizing preservation of the range of model uncertainty, scenario uncertainty, and interannual variability in the full CMIP6 ensemble results. This method is intended to help impact modelers select climate information from the CMIP archive efficiently for use in downstream models that require global coverage of climate information. This is particularly critical for large-ensemble experiments of multisector dynamic models that may be varying additional features beyond climate inputs in a factorial design, thus putting constraints on the number of climate simulations that can be used. We focus on temperature and precipitation outputs of CMIP6 models, as these are two of the most used variables among impact models, and many other key input variables for impacts are at least correlated with one or both of temperature and precipitation (e.g., relative humidity). Besides preserving the multi-model ensemble variance characteristics, we prioritize selecting CMIP6 models in the subset that preserve the very likely distribution of equilibrium climate sensitivity values as assessed by the latest Intergovernmental Panel on Climate Change (IPCC) report. This approach could be applied to other output variables of climate models and, possibly when combined with emulators, offers a flexible framework for designing more efficient experiments on human-relevant climate impacts. It can also provide greater insight into the properties of existing CMIP6 models.

Snyder, Abigail C.↗

Report of the Atmosphere to Electrons Land-Based Mesoscale-to-Microscale Coupling Project (FY2020)

The overall goal of the Mesoscale-to-Microscale Coupling (MMC) project is to improve coupling between mesoscale and microscale simulations via improved guidance and new strategies for setting up simulations and for the development of new tools that can be used across the community. Including the mesoscale forcing is critical to modeling the full energy transfer across scales in the atmosphere. The project-specific objectives include: 1) To apply rigorous verification and validation (V&V) techniques to the new modeling tools that are developed as part of the project to ensure the accuracy of our codes and results and to develop estimates of the relative uncertainty. 2) To improve computational performance of the coupled MMC models through the development of methods that can be used to reduce turbulence spin-up time and hence the size of computational domains. 3) To improve representation of the surface layer in microscale models to enhance simulations of hub-height wind speed. 4) To develop guidance for the community describing the best ways to couple mesoscale and microscale models, including specific spatial scales at which the handoff to the microscale model should occur. 5) To prepare documentation and a suite of software tools that can be used across the community. And 6) To transition MMC research to the offshore environment. Major progress was made in each of these areas during FY20. The land-based portion of the project is reported herein, while the offshore portion will be reported separately. The team continued to advance the MMC tools and methodologies as well as to document their performance in journal papers and conference presentations, although several planned conferences were cancelled due to the COVID-19 pandemic.

54 ENVIRONMENTAL SCIENCES↗

Axial H-Bonding Solvent Controls Inhomogeneous Spectral Broadening, While Peripheral H-Bonding Solvent Controls Vibronic Broadening: Cresyl Violet in Methanol

The dynamics of the nuclei of both a chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore–solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both chromophore and its condensed-phase environment control many spectral features, including the vibronic and inhomogeneous broadening present in spectral line shapes. For the cresyl violet chromophore in methanol, we here analyze and isolate the effect of specific chromophore–solvent interactions on simulated spectral densities, reorganization energies, and linear absorption spectra. Employing both force field and ab initio molecular dynamics trajectories along with the inclusion of only certain solvent molecules in the excited-state calculations, we determine that the methanol molecules axial to the chromophore are responsible for the majority of inhomogeneous broadening, with a single methanol molecule that forms an axial hydrogen bond dominating the response. Furthermore, the strong peripheral hydrogen bonds do not contribute to spectral broadening, as they are very stable throughout the dynamics and do not lead to increased energy-gap fluctuations. We also find that treating the strong peripheral hydrogen bonds as molecular mechanical point charges during the molecular dynamics simulation underestimates the vibronic coupling. Including these peripheral hydrogen bonding methanol molecules in the quantum-mechanical region in a geometry optimization increases the vibronic coupling, suggesting that a more advanced treatment of these strongly interacting solvent molecules during the molecular dynamics trajectory may be necessary to capture the full vibronic spectral broadening.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Capturing electronic decoherence in quantum-classical dynamics using the ring-polymer-surface-hopping–density-matrix approach

Simulations of coupled electronic and nuclear dynamics in molecules can be quite challenging due to the involved interplay of the many degrees of freedom. Because a full quantum treatment of both electrons and nuclei is computationally very demanding, it is generally restricted to model systems or rather small molecules and short timescales. Mixed quantum-classical dynamics methods such as Tully's fewest switches surface hopping (FSSH) can be used to overcome this limitation. However, FSSH is known to poorly describe electronic coherences and decoherence phenomena. Here, we present an approach that combines FSSH with ring-polymer molecular dynamics (RPMD) in a specific way that aims to alleviate the coherence problem. Termed the ring-polymer-surface-hopping–density-matrix approach, this method uses an electronic density-matrix formulation to calculate surface hopping rates. Additionally, this incorporates decoherence effects into FSSH in a natural way by taking into account the spatial spreading of the ring polymer that mimics the width of a nuclear wave packet in each RPMD trajectory. By applying our method to Tully's one-dimensional model system, we demonstrate that this method captures a crucial decoherence mechanism that is missing in FSSH. Furthermore, our method turns out to be superior at describing electronic coherences compared with earlier attempts at combining RPMD and FSSH.

74 ATOMIC AND MOLECULAR PHYSICS↗