Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic”

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 1,801 records · Page 100
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Real-Space Local Dynamics of Molten Inorganic Salts Using Van Hove Correlation Function

Molten inorganic salts are attracting resurgent attention because of their unique physicochemical properties, making them promising media for next-generation concentrating solar power systems and molten salt reactors. The dynamics of these highly disordered ionic media is largely studied by theoretical simulations, while the robust experimental techniques capable of observing local dynamics are not well-developed. To provide fundamental insights into the atomic-scale transport properties of molten salts, we report the real-space dynamics of molten magnesium chloride at high temperatures employing the Van Hove correlation function obtained by inelastic neutron scattering. Overall, our results directly depict the distance-dependent dynamics of a molten salt on the picosecond time scale. This study demonstrates the capability of the developed approach to describe the locally correlated- and self-dynamics in molten salts, significantly improving our understanding of the interplay between microscopic structural parameters and their dynamics that ultimately control physical properties of condensed matter in extreme environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Structure and Dynamics of Bimodal Colloidal Dispersions in a Low-Molecular-Weight Polymer Solution

We present an experimental study of the structural and dynamical properties of bimodal, micrometersized colloidal dispersions (size ratio ≈ 2) in an aqueous solution of low-molecular weight polymer (polyethylene glycol 2000) using synchrotron ultra-small angle X-ray scattering (USAXS) and USAXSbased X-ray photon correlation spectroscopy. We fixed the volume fraction of the large particles at 5 % and systematically increased the volume fraction of the small particles from 0 % to 5 % to evaluate its effect on the structure and dynamics. The bimodal dispersions were homogenous through the investigated parameter space. We found that the partial structure factors can be satisfactorily retrieved for the bimodal colloidal dispersions using a Percus-Yevick hard sphere potential when the particle size distributions of the particles were taken into account. We also found that the partial structure factor between the large particles does not exhibit significant variation with increasing volume fraction of small particles, whereas the isothermal compressibility of the binary mixture was found to decrease with increasing volume fraction of small particles. The dynamics of single-component large particle dispersion obey the principles of de Gennes narrowing, where the wave vector dependence of the interparticle diffusion coefficient is inversely proportional to the interparticle structure factor. The dynamics of the bimodal dispersions demonstrate strong dependence on the fraction of small particles. As a result, we also made a comparison between the experimental effective dynamic viscosity of the bimodal dispersion with theoretical predictions, which suggest that the complex mutual interactions between large and small particles have a strong effect on the dynamic behaviors of bimodal dispersions.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Untangleing the effects of chain rigidity on the structure and dynamics of strongly adsorbed polymer melts

Here, we present a detailed analysis of coarse-grained molecular dynamics simulations of semiflexible polymer melts in contact with a strongly adsorbing substrate. We have characterized the segments in the interfacial layer by counting the number of trains, loops, tails and unadsorbed segments. For more rigid chains, a tail and an adsorbed segment (a train) dominate while loops are more prevalent in more flexible chains. The tails exhibit a non-uniformly stretched conformation akin to the polydispersed pseudobrush envisioned by Guiselin. To probe the dynamics of the segments we computed the layer z-resolved intermediate coherent collective dynamics structure factor, S(q, t, z), mean-square displacement of segments, and the 2nd Legendre polynomial of the time-autocorrelation of unit bond vectors, 2 [n i (t,z)•n i (0,z)]>. Our results show that segmental dynamics is slower for stiffer chains and there is a strong correlation between the structure and dynamics in the interfacial layer. There is no glassy layer, and the slowing down in dynamics of stiffer chains in the adsorbed region can be attributed to the densification and the more persistent layering of segments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Structural and Dynamical Insights into the Formation Process of a Cross-Linked Polymer Network in Acrylic Adhesives During Thermal Curing

Many modern adhesives, sealants, and coatings rely on the controlled transition from a liquid to a solid state by forming a three-dimensional cross-linked polymer network, often referred to as curing. The curing process, which is initiated by the mixing of reactive components or an external trigger, defines the structure of a network and further controls the final mechanical properties of cured materials. However, the curing mechanism is not fully understood yet due to the lack of experimental tools capable of directly probing the structure and dynamics of a network over relevant time- and length scales. Here, in this paper, we report the curing process of a commercial two-component methyl methacrylate (MMA) adhesive using in operando X-ray photon correlation spectroscopy (XPCS), a method that closely simulates the target manufacturing environment of the adhesive. The results are then integrated with those obtained by rheology, differential scanning calorimetry (DSC), and transmission electron microscopy to establish the structure–dynamics–process–property relationship. The XPCS results identify four distinct stages in the curing process after the mixing, extrusion, and deposition of the acrylic adhesive: (i) At a cure time (or “aging time”, t age ) of less than 1 min, nanodomains of polymerized MMA are formed within a liquid monomeric MMA matrix. The average size is several nm and remains constant over t age , while the dynamics of the nanodomains are slowed down with t age due to an increase in the viscosity of the MMA matrix. (ii) After t age > 1 min, the size of the nanodomains increases with t age until the gel point (= 6.3 min after mixing as determined by rheology). The dynamics of the nanodomains also increase due to the heat generated by the exothermic reaction. (iii) At the gel point, the nanodomains begin to interconnect each other, resulting in a network structure with a characteristic length of about 100 nm. This characteristic network size does not change for the rest of the curing process up to t age = 500 min. The dynamics of the network structure, however, show a rapid slowing down with t age up to t age ≈ 12 min, corresponding to the onset of vitrification (as determined by rheology). (iv) At t age > 12 min, when the DSC and rheology data can no longer provide meaningful information, the XPCS data show a further slowing down of the network dynamics associated with vitrification. Our results provide rich and complex insights into the physics and material design of thermosets in commercially relevant processes, which are essential for future industrial applications.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Chain and Ion Dynamics in Precise Polyethylene Ionomers

In this paper, we analyze the dynamics from microsecond-long, atomistic molecular dynamics (MD) simulations of a series of precise poly(ethylene- co -acrylic acid) ionomers neutralized with lithium, with three different spacer lengths between acid groups on the ionomers and at two temperatures. At short times, the intermediate structure factor calculated from the MD simulations is in reasonable agreement with quasi-elastic neutron scattering data for partially neutralized ionomers. For ionomers that are 100% neutralized with lithium, the simulations reveal three dynamic processes in the chain dynamics. The fast process corresponds to hydration librations, the medium-time process corresponds to local conformational motions of the portions of the chains between ionic aggregates, and the long-time process corresponds to relaxation of the ionic aggregates. At 600 K, the dynamics are sufficiently fast to observe the early stages of lithium-ion motion and ionic aggregate rearrangements. In the partially neutralized ionomers with isolated ionic aggregates, the Li-ion-containing aggregates rearrange by a process of merging and breaking up, similar to what has been observed in coarse-grained (CG) simulations. In the 100% neutralized ionomers that contain percolated ionic aggregates, the chains remain pinned by the percolated aggregate at long times, but the lithium ions are able to move along the percolated aggregate. In this work, the ion dynamics are also qualitatively similar to those seen in previous CG simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Microscopic Origins of the Nonlinear Behavior of Particle-Filled Rubber Probed with Dynamic Strain XPCS

The underlying microscopic response of filler networks in reinforced rubber to dynamic strain is not well understood due to the experimental difficulty of directly measuring filler network behavior in samples undergoing dynamic strain. This difficulty can be overcome with in situ X-ray photon correlation spectroscopy (XPCS) measurements. The contrast between the silica filler and the rubber matrix for X-ray scattering allows us to isolate the filler network behavior from the overall response of the rubber. This in situ XPCS technique probes the microscopic breakdown and reforming of the filler network structure, which are responsible for the nonlinear dependence of modulus on strain, known in the rubber science community as the Payne effect. These microscopic changes in the filler network structure have consequences for the macroscopic material performance, especially for the fuel efficiency of tire tread compounds. In this work, we elucidate the behavior with in situ dynamic strain XPCS experiments on industrially relevant, vulcanized rubbers filled (13 vol %) with novel air-milled silica of ultrahigh-surface area (UHSA) (250 m 2 /g). The addition of a silane coupling agent to rubber containing this silica causes an unexpected and counterintuitive increase in the Payne effect and decrease in energy dissipation. For this rubber, we observe a nearly two-fold enhancement of the storage modulus and virtually equivalent loss tangent compared to a rubber containing a coupling agent and conventional silica. Interpretation of our in situ XPCS results simultaneously with interpretation of traditional dynamic mechanical analysis (DMA) strain sweep experiments reveals that the debonding or yielding of bridged bound rubber layers is key to understanding the behavior of rubber formulations containing the silane coupling agent and high-surface area silica. These results demonstrate that the combination of XPCS and DMA is a powerful method for unraveling the microscale filler response to strain which dictates the dynamic mechanical properties of reinforced soft matter composites. With this combination of techniques, we have elucidated the great promise of UHSA silica when used in concert with a silane coupling agent in filled rubber. Such composites simultaneously exhibit large moduli and low hysteresis under dynamic strain.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Improving the MJO Forecast of S2S Operation Models by Correcting Their Biases in Linear Dynamics

The operational dynamic subseasonal to seasonal (S2S) models for Madden-Julian oscillation (MJO) forecasting mostly still suffer from systematic errors in capturing the MJO's key dynamic features, such as its growth rate and propagation speed. By deriving the linear dynamic operators using the linear inverse modeling (LIM) approach, we propose a method to partly correct the errors in MJO linear dynamic operators to improve the MJO predictions of three operational dynamic S2S models. Correcting the deficiencies of the too-fast decay rates and the unrealistic propagating phase speeds lead to MJO prediction skills being extended by approximately 2–4 days. The improvements are more significant for the models with larger biases in MJO amplitude and propagation. This approach in principle may be extendable to predictions of other types of climate variability such as ENSO on one hand, and possible inclusions of nonlinear dynamics effects on the other hand.

58 GEOSCIENCES↗
DOE OSTI2022

Using Ensemble Data Assimilation to Estimate Transient Hydrologic Exchange Flow Under Highly Dynamic Flow Conditions

Abstract Quantifying dynamic hydrologic exchange flows (HEFs) within river corridors that experience high‐frequency flow variations caused by dam regulations is important for understanding the biogeochemical processes at the river water and groundwater interfaces. Heat has been widely used as a tracer to infer steady‐state flow velocities through analytical solutions of heat transport defined by the diurnal temperature signals. Under sub‐daily dynamic flow conditions, however, such analytical solutions are not applicable due to the violation of their fundamental assumptions. In this study, we developed a data assimilation‐based approach to estimate the sub‐daily flux under highly dynamic flow conditions using multi‐depth temperature observations at a 5‐min resolution. If the hydraulic gradient is measured, Darcy's law was used to calculate the flux with permeability estimated from temperature responses below the riverbed. Otherwise, flux was estimated directly by assimilating multi‐depth temperature data at 1‐ or 2‐hr time intervals assuming one‐dimensional flow and heat transport governing equation. By comparing estimated fluxes with model‐generated synthetic truth, we demonstrated that both schemes have robust performance in estimating fluxes under highly dynamic flow conditions. This data assimilation‐based flux estimation method was able to capture the vertical sub‐daily fluxes using multi‐depth high‐resolution temperature data alone, even in the presence of multi‐dimensional flow. This approach has been successfully applied to real field temperature data collected at the Hanford site, which experiences highly dynamic HEFs. Our study shows the promise of adopting distributed 1‐D temperature monitoring to capture spatial and temporal exchange dynamics in river corridors at a watershed scale or beyond.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Mapping in situ the assembly and dynamics in aqueous supramolecular polymers

Supramolecular polymers, bonded through directional non-covalent interactions, closely mimic dynamic behaviors of biological nanofibers. However, the complexity of assembly pathways makes it highly challenging to unravel the nature of supramolecular dynamics in aqueous environments. Here we introduce a precise combinatorial titration methodology to probe in situ the assembly of peptide amphiphiles (PAs). This approach reveals a binary assembly mechanism governed by equilibrium between spheroidal micelles and β-sheet polymers. Weakening hydrogen bonding shifts the equilibrium towards micelles and decreases the internal structural order of filamentous polymers, promoting supramolecular dynamics. Extending this methodology to two-component copolymerization systems, we find a surprising tendency to form blocky nanostructures with reduced internal phase separation as the mismatch in peptide sequence decreases. Interestingly, while well-mixed copolymers acquire different dynamics, mismatched ones retain the characteristic supramolecular motion of their homopolymer counterparts. These critical insights into supramolecular dynamics offer strategies to tailor the dynamic functions of supramolecular nanomaterials.

Molecular self-assembly↗
DOE OSTI2025

Advancing simulations of coupled electron and phonon nonequilibrium dynamics using adaptive and multirate time integration

Electronic structure calculations in the time domain provide a deeper understanding of nonequilibrium dynamics in materials. The real-time Boltzmann equation (rt-BTE), used in conjunction with accurate interactions computed from first principles, has enabled reliable predictions of coupled electron and lattice dynamics. However, the timescales and system sizes accessible with this approach are still limited, with two main challenges being the different timescales of electron and phonon interactions and the cost of computing collision integrals. As a result, only a few examples of these calculations exist, mainly for two-dimensional (2D) materials. Here we leverage adaptive and multirate time integration methods to achieve a major step forward in solving the coupled rt-BTEs for electrons and phonons. Relative to conventional (non-adaptive) time-stepping, our approach achieves a 10x speedup for a target accuracy, or greater accuracy by 3–6 orders of magnitude for the same computational cost, enabling efficient calculations in both 2D and bulk materials. This efficiency is showcased by computing the coupled electron and lattice dynamics in graphene up to ~100 ps, as well as modeling ultrafast lattice dynamics and thermal diffuse scattering maps in bulk materials (silicon and gallium arsenide). In addition to improved efficiency, our adaptive method can resolve the characteristic rates of different physical processes, thus naturally bridging different timescales. This enables simulations of longer timescales and provides a framework for modeling multiscale dynamics of coupled degrees of freedom in matter. Our work opens new opportunities for quantitative studies of nonequilibrium physics in materials, including driven lattice dynamics with phonons coupled to electrons, spin, and other degrees of freedom.

Yao, Jia [California Institute of Technology (CalT↗
DOE OSTI2025

Dynamic nanodomains dictate macroscopic properties in lead halide perovskites

Lead halide perovskites have emerged as promising materials for solar energy conversion and X-ray detection owing to their remarkable optoelectronic properties. However, the microscopic origins of their superior performance remain unclear. Here we show that low-symmetry dynamic nanodomains present in the high-symmetry average cubic phases, whose characteristics are dictated by the A-site cation, govern the macroscopic behaviour. We combine X-ray diffuse scattering, inelastic neutron spectroscopy, hyperspectral photoluminescence microscopy and machine-learning-assisted molecular dynamics simulations to directly correlate local nanoscale dynamics with macroscopic optoelectronic response. Our approach reveals that methylammonium-based perovskites form densely packed, anisotropic dynamic nanodomains with out-of-phase octahedral tilting, whereas formamidinium-based systems develop sparse, isotropic, spherical nanodomains with in-phase tilting, even when crystallography reveals cubic symmetry on average. We demonstrate that these sparsely distributed isotropic nanodomains present in formamidinium-based systems reduce electronic dynamic disorder, resulting in a beneficial optoelectronic response, thereby enhancing the performance of formamidinium-based lead halide perovskite devices. By elucidating the influence of the A-site cation on local dynamic nanodomains, and consequently, on the macroscopic properties, we propose leveraging this relationship to engineer the optoelectronic response of these materials, propelling further advancements in perovskite-based photovoltaics, optoelectronics and X-ray imaging.

Materials Science↗
DOE OSTI2025

Liquid-like dynamics in a solid-state lithium electrolyte

Superionic materials represent a regime intermediate between the crystalline and liquid states of matter. Despite the considerable interest in potential applications for solid-state batteries or thermoelectric devices, it remains unclear whether the fast ionic diffusion observed in superionic materials reflects liquid-like dynamics or whether the hops of mobile ions are inherently coupled to more conventional lattice phonons. Here we reveal a crossover from crystalline vibrations to relaxational dynamics of ionic diffusion in the superionic compound Li 6 PS 5 Cl, a candidate solid-state electrolyte. By combining inelastic and quasi-elastic neutron-scattering measurements with first-principles-based machine-learned molecular dynamics simulations, we found that the vibrational density of states in the superionic state strongly deviates from the quadratic behaviour expected from the Debye law of lattice dynamics. The superionic dynamics emerges from overdamped phonon quasiparticles to give rise to a linear density of states characteristic of instantaneous normal modes in the liquid state. Further, we showed that the coupling of lattice phonons with a dynamic breathing of the Li+ diffusion bottleneck enables an order-of-magnitude increase in diffusivity. Thus, our results shed insights into superionics for future energy storage and conversion technologies.

25 ENERGY STORAGE↗
DOE OSTI2018

The ultrafast X-ray spectroscopic revolution in chemical dynamics

There have been rapid developments in short-pulse X-ray sources in the past two decades, which have enabled the study of nuclear and electronic dynamics by ultrafast X-ray spectroscopies with unprecedented time resolution ranging from nanoseconds to attoseconds. In this Perspective, we discuss some of the major achievements in the study of nuclear and electronic dynamics with X-ray pulses produced by high-harmonic, free-electron-laser and synchrotron sources. The particular dynamic processes probed by X-ray radiation highlighted in this Perspective are electronic coherences on attosecond to femtosecond timescales, chemical reactions, such as dissociations, and pericyclic ring-openings, spin-crossover dynamics, ligand-exchange dynamics and structural deformations in excited states. X-ray spectroscopic probing of chemical dynamics holds promising use for the future owing to the ongoing developments of new spectroscopies, such as four-wave mixing, and the continuous improvements in emerging laboratory-based, high-harmonic sources and large-scale, facility-based, free-electron lasers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Subterahertz collective dynamics of polar vortices

The collective dynamics of topological structures are of interest from both fundamental and applied perspectives. For example, studies of dynamical properties of magnetic vortices and skyrmions have not only deepened our understanding of many-body physics but also offered potential applications in data processing and storage. Topological structures constructed from electrical polarization, rather than electron spin, have recently been realized in ferroelectric superlattices, and these are promising for ultrafast electric-field control of topological orders. However, little is known about the dynamics underlying the functionality of such complex extended nanostructures. In this study, using terahertz-field excitation and femtosecond X-ray diffraction measurements, we observe ultrafast collective polarization dynamics that are unique to polar vortices, with orders-of-magnitude higher frequencies and smaller lateral size than those of experimentally realized magnetic vortices. A previously unseen tunable mode, hereafter referred to as a vortexon, emerges in the form of transient arrays of nanoscale circular patterns of atomic displacements, which reverse their vorticity on picosecond timescales. Its frequency is considerably reduced (softened) at a critical strain, indicating a condensation (freezing) of structural dynamics. We use first-principles-based atomistic calculations and phase-field modelling to reveal the microscopic atomic arrangements and corroborate the frequencies of the vortex modes. The discovery of subterahertz collective dynamics in polar vortices opens opportunities for electric-field-driven data processing in topological structures with ultrahigh speed and density.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Dynamical topological phase realized in a trapped-ion quantum simulator

Nascent platforms for programmable quantum simulation offer unprecedented access to new regimes of far-from-equilibrium quantum many-body dynamics in almost isolated systems. Here achieving precise control over quantum many-body entanglement is an essential task for quantum sensing and computation. Extensive theoretical work indicates that these capabilities can enable dynamical phases and critical phenomena that show topologically robust methods to create, protect and manipulate quantum entanglement that self-correct against large classes of errors. However, so far, experimental realizations have been confined to classical (non-entangled) symmetry-breaking orders. In this work, we demonstrate an emergent dynamical symmetry-protected topological phase, in a quasiperiodically driven array of ten 171 Yb + hyperfine qubits in Quantinuum’s System Model H1 trapped-ion quantum processor. This phase shows edge qubits that are dynamically protected from control errors, cross-talk and stray fields. Crucially, this edge protection relies purely on emergent dynamical symmetries that are absolutely stable to generic coherent perturbations. Here, this property is special to quasiperiodically driven systems: as we demonstrate, the analogous edge states of a periodically driven qubit array are vulnerable to symmetry-breaking errors and quickly decohere. Our work paves the way for implementation of more complex dynamical topological orders that would enable error-resilient manipulation of quantum information.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

Label-free functional analysis of root-associated microbes with dynamic quantitative oblique back-illumination microscopy

The increasing global demand for food, coupled with concerns about the environmental impact of synthetic fertilizers, underscores the urgency of developing sustainable agricultural practices. Nitrogen-fixing bacteria, known as diazotrophs, offer a potential solution by converting atmospheric nitrogen into bioavailable forms, reducing the reliance on synthetic fertilizers. However, a deeper understanding of their interactions with plants and other microbes is needed. In this study, we introduce a recently developed label-free 3D quantitative phase imaging technology called dynamic quantitative oblique back-illumination microscopy (DqOBM) to assess the functional dynamic activity of diazotrophs in vitro and in situ. Our experiments involved three different diazotrophs (Sinorhizobium meliloti, Azotobacter vinelandii, and Rahnella aquatilis) cultured on media with amendments of carbon and nitrogen sources. Over 5 days, we observed increased dynamics in nutrient-amended media. These results suggest that the observed bacterial dynamics correlate with their metabolic activity. Furthermore, we applied qOBM to visualize microbial dynamics within the root cap and elongation zone of Arabidopsis thaliana primary roots. This allowed us to identify distinct areas of microbial infiltration in plant roots without the need for fluorescent markers. Our findings demonstrate that DqOBM can effectively characterize microbial dynamics and provide insights into plant-microbe interactions in situ, offering a valuable tool for advancing our understanding of sustainable agriculture.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2017

Atomistic non-adiabatic dynamics of the LH2 complex with a GPU-accelerated ab initio exciton model

We recently outlined an efficient multi-tiered parallel ab initio excitonic framework that utilizes time dependent density functional theory (TDDFT) to calculate ground and excited state energies and gradients of large supramolecular complexes in atomistic detail – enabling us to undertake non-adiabatic simulations which explicitly account for the coupled anharmonic vibrational motion of all the constituent atoms in a supramolecular system. Here we apply that framework to the 27 coupled bacterio-chlorophyll-a chromophores which make up the LH2 complex, using it to compute an on-the-fly nonadiabatic surface-hopping (SH) trajectory of electronically excited LH2. Part one of this article is focussed on calibrating our ab initio exciton Hamiltonian using two key parameters: a shift δ, which corrects for the error in TDDFT vertical excitation energies; and an effective dielectric constant ε, which describes the average screening of the transition-dipole coupling between chromophores. Using snapshots obtained from equilibrium molecular dynamics simulations (MD) of LH2, we tune the values of both δ and ε through fitting to the thermally broadened experimental absorption spectrum, giving a linear absorption spectrum that agrees reasonably well with experiment. In part two of this article, we construct a time-resolved picture of the coupled vibrational and excitation energy transfer (EET) dynamics in the sub-picosecond regime following photo-excitation. Assuming Franck–Condon excitation of a narrow eigenstate band centred at 800 nm, we use surface hopping to follow a single nonadiabatic dynamics trajectory within the full eigenstate manifold. Consistent with experimental data, this trajectory gives timescales for B800→B850 population transfer (τ B800→B850 ) between 650–1050 fs, and B800 population decay (τ 800→ ) between 10–50 fs. The dynamical picture that emerges is one of rapidly fluctuating LH2 eigenstates that are delocalized over multiple chromophores and undergo frequent crossing on a femtosecond timescale as a result of the atomic vibrations of the constituent chromophores. The eigenstate fluctuations arise from disorder that is driven by vibrational dynamics with multiple characteristic timescales. The scalability of our ab initio excitonic computational framework across massively parallel architectures opens up the possibility of addressing a wide range of questions, including how specific dynamical motions impact both the pathways and efficiency of electronic energy-transfer within large supramolecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Refine your search to explore more results.