Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “time dependent behavior”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Performance Evaluation and Optimal Design Analysis of Continuous-Operation Solar-Driven Cooling Absorption Systems With Thermal Energy Storage

A solar absorption cooling system consisting of a flat plate collector, thermal energy storage tank, and absorption chiller is analyzed in this work. A dimensionless model is developed from the energy balance on each component and the chiller's characteristic performance curves. The model is used to determine the interaction and influence of different parameters such as tank size, solar collector area, chiller size, cooling load, cooling temperature, heat loss, and mass flow rates on the performance. From the analysis, smaller solar collector areas are required for lower cooling loads and smaller tank volumes. A specific cooling load of 1 x 10^-5 will require a specific solar collector area between two and six times larger, depending on the initial tank temperature, than the area required for a baseline system that considers typical commercial design and operation parameters. A similar behavior was observed for the specific tank volume. For the baseline system, the minimum specific area of the collector of 9.57 is achieved for an initial tank temperature of 1.19. For a cooling load of 1 x 10^-5, the optimum initial tank temperature will be 1.11 that results in a minimum specific solar collector area of 25.26. A specific tank volume of 4 x 10^-4 will also have an optimum initial tank temperature of 1.11 that minimizes the specific solar collector area to a value of 28.18. The approach and analysis in this work can be used to determine design parameters for solar absorption cooling systems based on a proper relation among system's dimensions to achieve optimum operation.

cooling load↗
DOE OSTI2023

Photoionization of C 60 at high energies

Calculations of the photoionization cross section of C60 including all valence and core subshells have been carried out for photon energies up to 1 keV using density functional theory (DFT) and time-dependent DFT (TDDFT) within the framework of a fully molecular model. The high-energy valence subshell cross sections behave rather differently than the results from model potentials where the cross section falls far too rapidly with energy. This unphysical behavior of the model calculations is traced to smearing out of the carbon nuclei in the model potentials which makes it difficult for momentum to be conserved at the higher energies. Comparison with 60 times carbon atomic cross section sheds light on similarities and differences. The high energy behavior of individual valence cross sections reflects the amount of C 2s/2p mixing in the molecular orbitals.

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

Compound pulse characteristics of a heterogeneous composite scintillator in a gamma-ray field

An important aspect of radiation detection in scintillators is the time dependence of emitted fluorescence, which results in characteristic measured pulse shape. In some materials, a relatively strong dependence of the scintillation pulse shape on stopping power exists, which can provide reliable identification of the incident particle. An alternative method to realize particle identification is by combining materials with different fluorescence lifetimes into heterogeneous structures. Such composite scintillators derive their properties from both the constituent material characteristics and the geometry of the composite structure. In composite scintillators with a high degree of heterogeneous loading, a superposition of fluorescence contributions originating in multiple materials, which is observed as single waveform with an intermediate pulse shape, may affect its ability to reject rays. We measure the time dependence of light output from a scintillator composed of 6 Li-containing glass shards and scintillating polyvinyl toluene and identify events that exhibit such compound behavior when exposed to rays and fast neutrons. We develop a modeling and simulation framework that reproduces the pulse shapes in heterogeneous scintillators and use it study the effect that the weight percentage of 6 Li glass has on rejection. The modeling framework is applied to the experimentally studied scintillator, finding a good agreement. Here, the developed modeling and simulation approach will help optimize the design of heterogeneous scintillators to meet the desired trade-o between neutron capture efficiency and -ray rejection.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

In situ probing of structure and deagglomeration of SnO 2 colloids via small-angle X-ray scattering

Transforming dry nanopowders into stable colloidal dispersions remains challenging due to the cohesive forces between the nanoparticles (NPs) which promote agglomeration. Effective dispersion and deagglomeration of these agglomerates is a critical process in the formulation and preparation of nanoparticle-based functional materials via the colloidal route. Understanding the deagglomeration dynamics provides information to improve microstructural quality in various applications by enabling engineering of agglomerate size, structure and morphology. However, the deagglomeration process dynamics with respect to the evolution of the fractal agglomerate structures, particularly for very small NPs, is still poorly understood and requires further investigation. This study employs in situ small-angle X-ray scattering (SAXS) to investigate the sonication-induced deagglomeration of SnO 2 NPs. Electrostatically stabilized SnO 2 colloids with varying primary particle size (6–21 nm) are investigated in a specifically designed in situ cell using synchrotron-based SAXS to study the influence of sonication time and intensity on the nanoscaled agglomerates. Complete structural analysis via SAXS reveals a direct correlation between changes in agglomerate size and structure, size-dependent deagglomeration behavior and a dependence on the overall energy introduced during sonication into the dispersion regardless of the actual power as well as ultrasonic process parameters in case of SnO 2 NPs. The results suggest that control of the dispersion process during ultrasonic deagglomeration results in tailoring agglomerates with respect to size and structure.

Deagglomeration↗
DOE OSTI2020

Interplay between Intrinsic Thermal Stability and Expansion Properties of Functionalized UiO-67 Metal–Organic Frameworks

The UiO family of metal-organic frameworks (MOFs) has been extensively studied for several applications owing to their high stability rendered by their robust secondary building units. The efficient design and use of these materials require a fundamental understanding of their thermal stability and its impact on chemical and structural functionality. In this work, we provide a detailed characterization of the intrinsic thermal behavior of the UiO-67 and functional analogues, UiO-67-NH 2 and UiO-67-CH 3 . Using in situ temperature-programmed X-ray diffraction, we find that distortion of the carboxylate group on the organic linker leads to negative thermal expansion (NTE) of the UiO-67 MOFs during heating. This NTE behavior is correlated with rich and reversible thermal changes observed in the MOF infrared spectral signature as samples are heated to the sample activation temperature (473 K). We find that in the absence of oxygen, activated UiO-67 samples show higher thermal stability compared to ambient or inert environments, with temperature-programmed desorption revealing an overall stability trend: UiO-67 > UiO-67- CH 3 > UiO-67-NH 2 . Two stages of change are observed during thermal treatment above 473 K, which are directly related to deformation of the inorganic node and the isotropic NTE behavior of these materials. Ultimately, these results provide a real-time interpretation of the fundamental thermoresponsive behavior of UiO-67 MOFs and offer a foundation for accurate interpretation of MOF interactions with guest molecules and their temperature dependence.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Systematic Study of Solid-State U(VI) Photoreactivity: Long-Lived Radicalization and Electron Transfer in Uranyl Tetrachloride

Reported are the syntheses, structural characterizations and luminescence properties of three novel [UO 2 Cl 4 ] 2- bearing compounds containing substituted 1,1’-dialkyl-4,4’-bipyridinum dications (i.e. viologens). These compounds undergo photoinduced luminescence quenching upon exposure to UV radiation. Further, this reactivity is concurrent with two phenomena: radicalization of the uranyl tetrachloride anion and photoelectron transfer to the viologen which constitutes the formal transfer of one electron from the [UO 2 Cl 4 ] 2- to the viologen species. This behavior is elucidated using electron paramagnetic resonance (EPR) spectroscopy and further probed through a series of characterization and computational techniques including Rehm-Weller analysis, time-dependent density functional theory (TD-DFT), and density of states (DOS). This work provides a systematic study of the photoreactivity of the uranyl unit in the solid state, an under-described aspect of fundamental uranyl chemistry.

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

Pulse-shape analysis in Cd 0.9 Zn 0.1 Te 0.98 Se 0.02 room-temperature radiation detectors

The effect of trapping of charge carriers on the shape of the charge pulses from room temperature nuclear detectors, based on a recently discovered quaternary semiconductor Cd 0.9 Zn 0.1 Te 0.98 Se 0.02 (CZTS), has been studied. Here, we present a charge-trapping model that explains the pulse shapes for its entire duration from radiation detectors containing multiple defect types. A piecewise continuous model, based on the movement of a large number of charge carriers rather than a single charge, has been proposed to explain the pulse shapes from planar room temperature nuclear detectors. An experimental method based on the model has been described to extract charge-trapping information. This model can be applied to similar semiconductor detectors with electrically active trap centers, in general, to extract information like pre-trapping drift times. Charge pulses from a CZTS planar detector, exposed to 5486 keV alpha particles emitted by a 241Am source, were acquired using a digitizer, and the pulse shapes were explained based on the proposed model. The pre-trapping drift times for electrons as well as holes in CZTS were calculated and studied as a function of operating bias voltage, and a contrast between their bias dependent behaviors was noticed, indicating the behavioral difference of electron and hole traps.

42 ENGINEERING↗
DOE OSTI2020

Modular Linear Thermoviscoelastic Model

Time-dependent, viscoelastic responses of materials like polymers and glasses have long been studied. As such, a variety of models have been put forth to describe the behavior including simple rheological models (e.g. Maxwell, Kelvin), linear "fading memory" theories, and hereditary integral based linear thermal viscoelastic approaches as well as more recent nonlinear theories that are either integral, fictive temperature, or differential internal state variable based. The current work details a new LINEAR_THERMOVISCOELASTIC model that has been added to LAME. This formulation represents a viscoelastic theory that neglects some of the phenomenological details of the PEC/SPEC models in favor of efficiency and simplicity. Furthermore, this new model is a first step towards developing "modular" viscoelastic capabilities akin to those available with hardening descriptions for plasticity models in LAME. Specifically, multiple different (including user-defined) shift-factor forms are implemented with each being easily selected via parameter specification rather than requiring distinct material models.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Spatiotemporal mapping of mesoscopic liquid dynamics

The study of liquid dynamics at mesoscopic scales is still strewn with difficulty due to limitations in theory and experiment. Historically, significant attention has been given to the analysis of space-time correlation functions and their frequency-Fourier transforms at a few discrete wave numbers. The massive computing power afforded by modern high performance computing clusters and the advent of a wide-angle neutron spin-echo spectrometer, however, have unlocked a more intuitive and fruitful approach to this problem. Using molecular dynamics simulations, here we demonstrate the benefits of spatiotemporally mapping intermediate scattering functions on a dense grid of correlation times and wave numbers. Four model systems are investigated: a Lennard-Jones liquid, a coarse-grained bead-spring polymer, a molten sodium chloride, and a poly(ethylene oxide) melt. We show that the spatiotemporal mapping approach is particularly useful for elucidating the mesoscopic dynamics in these liquids, where several underlying mechanisms, such as molecular relaxations, hydrodynamic modes, and nonhydrodynamic excitations, are potentially at play. Therefore, compared to the traditional method, direct visualization of density space-time correlation functions on two-dimensional color maps permits appraisals of complicated dynamical behavior at mesoscales in a global manner. For example, the scaling relations between space and time for different types of molecular motions can be straightforwardly identified on these plots, without any model-dependent analysis. Additionally, we show how theoretical ideas regarding collective mesoscopic dynamics, such as the classical hydrodynamic theory, the convolution approximation, and a recently proposed phenomenological model, can be discussed in terms of the global features of spatiotemporal maps of intermediate scattering functions. The new perspective offered by the spatiotemporal mapping method should prove useful for the study of liquid dynamics in general.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Driven responses of periodically patterned superconducting films

Here, we simulate the motion of a commensurate vortex lattice in a periodic lattice of artificial circular pinning sites having different diameters, pinning strengths, and spacings using the time-dependent Ginzburg-Landau formalism. Above some critical DC current density J c , the vortices depin, and the resulting steady-state motion then induces an oscillatory electric field E (t) with a defect "hopping" frequency f 0 , which depends on the applied current density and the pinning landscape characteristics. The frequency generated can be locked to an applied AC current density over some range of frequencies, which depends on the amplitude of the DC as well as the AC current densities. Both synchronous and asynchronous collective hopping behaviors are studied as a function of the supercell size of the simulated system and the (asymptotic) synchronization threshold current densities determined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Effects of flow on uranium speciation in soils impacted by acidic waste fluids

Radioactive acidic liquid waste is a common byproduct of uranium (U) and plutonium (Pu) enrichment and recycling processes whose accidental and planned release has led to a significant input of U into soils and sediments across the world, including at the U.S. DOE's Hanford site (WA, USA). Because of the particularly hazardous nature of U, it is important to predict its speciation when introduced into soils and sediments by acidic waste fluids. Of fundamental importance are the coupled effects of acid-driven mineral transformation and reactive transport on U speciation. To evaluate the effect of waste-fluid residence time and co-associated dissolved phosphate concentrations on U speciation in impacted soils and sediments, uncontaminated surface materials (from the Hanford Site) were reacted with U-containing synthetic acidic waste fluids (pH 2) amended with dissolved phosphate concentrations in both batch (no flow) and flow-through column systems for 7–365 days. By comparing dissolved U behavior and solid phase speciation as a function of flow regimen, we found that the availability of proton-promoted dissolution products (such as Si) to sequester U into uranyl silicates was dependent on waste fluid-sediment contact time as uranyl silicates were not detected in short contact time flow-through systems but were detected in no-flow, long contact time, reactors. Moreover, the dominance of uranyl phosphate as neoprecipitate U scavenger (principally in the form of meta-ankoleite) in phosphate amended systems confirmed the importance of phosphate amendments for an efficient sequestration of U in the soils and sediments. Overall, our experiments suggest that the formation of uranyl silicates in soils impacted by acidic waste fluids is likely to be limited unless reaction products are allowed to accumulate in soil pores, highlighting the importance of investigating soil U speciation in flow-through, transport-driven systems as opposed to no-flow, batch systems. Finally, this study provides insights into uranium speciation and its potential changes under acidic conditions for better prediction of risks and subsequent development of efficient remediation strategies.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Three-dimensional simulations of reshocked inclined Richtmyer-Meshkov instability: Effects of initial perturbations

The effect of initial perturbations on the evolution of the inclined Richtmyer-Meshkov turbulent mixing layer before and after reshock initiated by a shock wave with Mach number 1.55 is investigated through three-dimensional (3D) simulations using the flash code. The 3D simulations aim to reproduce both predominantly single-mode and multimode interfaces between light and heavy gases (N 2 -CO 2 , Atwood number, A≈0.22; amplitude to wavelength ratio of 0.088) which were created in an inclined shock tube facility to analyze the effects of initial conditions on mixing development in the entire flow field. The two-dimensional center slices of 3D simulations are compared with the experimental results to validate the computational code. Mixing width, mixed mass, mixed-mass thickness, and circulation in addition to concentration fields are shown to be in good agreement with the experimental data. The three-dimensional density and vorticity fields are first presented to qualitatively describe the flow behavior before and after reshock. Several measured density/velocity-related quantities indicate that the growth of the mixing material is strongly dependent on initial conditions. Before reshock and at early times after reshock, flow is clearly maintaining the memory of initial perturbations. However, at late time after reshock, although the large wavelength feature still dominates the flow motion, and the morphology of the two different interfaces indicates several differences, by breakdown of large-scale coherent structures to much finer scales, the memory of small scales of the multimode initial perturbation is not as clear as pre-reshock. Regarding three-dimensionality of the flow, before reshock in the multimode case, the baroclinic vorticity production, circulation, turbulent kinetic energy, and turbulent mass flux suggest that the small-scale roll-up features along the large inclined wavelength quickly evolves in all three dimensions. The coherent vortex tubes break down to smaller wormlike vortex structures, and turbulent fluctuations in the out-of-plane dimension are comparable to the spanwise direction. After reshock, this three-dimensionality of mixing growth was observed in the flow for both initial conditions. The results of this work represent a significant extension of previous computational studies performed on this specific topic. A different code with a different numerical method is validated through comparison with the experimental data. The initial perturbations are directly measured from the experimental results. Moreover, the entire three-dimensional experimental shock tube domain is simulated, and more quantities are investigated to understand the mixing mechanism and instability evolution in all three dimensions.

42 ENGINEERING↗
DOE OSTI2020

Quasilinear dispersion in electronic band structure and high Seebeck coefficient in CuFeS 2 -based thermoelectric materials

The earth-abundant natural mineral chalcopyrite CuFeS 2 is a potential n-type thermoelectric material because of its large Seebeck coefficient at high carrier concentrations. For a long time, the large Seebeck coefficient of CuFeS 2 has been attributed to a large electron effective mass, but the reasons for this and the unusual carrier concentration dependent behavior have rarely been discussed. Here, we systematically investigated the special transport behavior of CuFeS 2 and found the classical parabolic band model to be inadequate in explaining it. Our experimental and theoretical studies indicate that there are two flat electronic pockets at the and Z points of the Brillouin zone near the conduction band edge of CuFeS 2 that dominate the charge transport. These electronic pockets result from nonparabolic quasilinearly dispersing bands that give rise to a linear wave vector dependent energy (E ~ k) and a carrier density dependent effective mass (m* ~ m 0 n 1/3 ). Such a strong carrier concentration dependent carrier effective mass results in the high Seebeck coefficient of CuFeS 2 compound under a large carrier density. Here, the work demonstrates that quasilinearly dispersing bands can give strongly enhanced Seebeck coefficient, and could be useful in optimizing the properties of thermoelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Beam Energy Dependence of Triton Production and Yield Ratio $N_t \times N_p / N^2_d$ in Au+Au Collisions at RHIC

We report the triton (t) production in midrapidity (|y| < 0.5) Au + Au collisions at $\sqrt{s_{NN}} = 7.7-200$ GeV measured by the STAR experiment from the first phase of the beam energy scan at the Relativistic Heavy Ion Collider. The nuclear compound yield ratio ($N_t \times N_p / N^2_d$), which is predicted to be sensitive to the fluctuation of local neutron density, is observed to decrease monotonically with increasing charged-particle multiplicity ($dN_{ch}/dη$) and follows a scaling behavior. The $dN_{ch}/dη$ dependence of the yield ratio is compared to calculations from coalescence and thermal models. Enhancements in the yield ratios relative to the coalescence baseline are observed in the 0%-10% most central collisions at 19.6 and 27 GeV, with a significance of 2.3σ and 3.4σ, respectively, giving a combined significance of 4.1σ. The enhancements are not observed in peripheral collisions or model calculations without critical fluctuation, and decreases with a smaller $p_T$ acceptance. The physics implications of these results on the QCD phase structure and the production mechanism of light nuclei in heavy-ion collisions are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Thixotropy, antithixotropy, and viscoelasticity in hysteresis

Thixotropy, antithixotropy, and viscoelasticity are three types of time-dependent dynamics that involve fundamentally different underlying physical processes. Here, we show that the three dynamics exhibit different signatures in hysteresis by examining the fingerprints of the simplest thixotropic kinetic model, a new antithixotropic model that we introduce here, and the Giesekus model. We start by showing that a consistent protocol to generate hysteresis loops is a discrete shear-rate controlled ramp that begins and ends at high shear rates, rather than at low shear rates. Using this protocol, we identify two distinguishing features in the resulting stress versus shear rate loops. The first is the direction of the hysteresis loops: clockwise for thixotropy, but counterclockwise for viscoelasticity and antithixotropy. A second feature is achieved at high ramping rates where all responses lose hysteresis: the viscoelastic response shows a stress plateau at low shear rates due to lack of stress relaxation, whereas the thixotropic and antithixotropic responses are purely viscous with minimal shear thinning or thickening. We establish further evidence for these signatures by experimentally measuring the hysteresis of Laponite suspensions, carbon black suspensions, and poly(ethylene oxide) solutions, each representing a historically accepted example of each class of material behavior. The signatures measured in experiments are consistent with those predicted by the three models. This study reveals different fingerprints in hysteresis loops associated with thixotropy, antithixotropy, and viscoelasticity, which may be helpful in distinguishing the three time-dependent responses.

Mechanics↗
DOE OSTI2025

Model-based analysis of solute transport and potential carbon mineralization in the active layer of a hillslope underlain by permafrost with seasonal variability and climate change

Permafrost carbon, stored in frozen organic matter across vast Arctic and sub-Arctic regions, represents a substantial and increasingly vulnerable carbon reservoir. As global temperatures rise, the accelerated thawing of permafrost releases greenhouse gases, exacerbating climate change. However, freshly thawed permafrost carbon may also experience lateral transport by groundwater flow to surface water recipients such as rivers and lakes, increasing the terrestrial-to-aquatic transfer of permafrost carbon. Mobilization and subsurface transport mechanisms are poorly understood and not accounted for in global climate models, leading to high uncertainties in the predictions of the permafrost carbon feedback. Here, we focus on a hillslope in Endalen Valley, Svalbard, as a representative example of a high-Arctic hillslope underlain by continuous permafrost. We analyze solute transport in the form of a non-reactive tracer representing dissolved organic carbon (DOC) using a physics-based numerical model with the objective to study governing cryotic and hydrodynamic transport mechanisms relevant for warming permafrost regions. We first analyze transport times for DOC pools at different locations within the active layer under present-day climatic conditions and proceed to study susceptibility for deeper ancient carbon release in the upper permafrost due to thaw under different warming scenarios. Results suggest that DOC in the active layer near the permafrost table experiences rapid lateral transport upon thaw due to saturated conditions and lateral flow, while DOC close to the ground surface experiences slower transport due to flow in unsaturated soil. Deeper permafrost carbon release exhibits vastly different transport behaviors depending on warming and thaw rate. Gradual warming leads to small fractions of DOC being mobilized every year, while the majority moves vertically through percolation and cryosuction. Abrupt thaw resulting from a single very warm year leads to faster lateral transport times, similar to active layer DOC released in saturated conditions. Lastly, we analyze the potential susceptibility of DOC to mineralization to CO 2 prior to export due to soil moisture and temperature conditions. We find that high liquid saturation during transport coincides with very low mineralization rates and potentially inhibits mineralization into CO 2 before export. Overall, the results highlight the importance of subsurface hydrologic and thermal conditions for the retention and lateral export of permafrost carbon by subsurface flow.

Hamm, Alexandra [Stockholm Univ. (Sweden)] (ORCID:↗
DOE OSTI2025

Comparison of the Distribution of the In Situ Ring Current Density Based on Magnetic Field Observations From THEMIS, MMS and Cluster

Abstract We review and extend previous work on the ring current density distribution by applying the curlometer technique to multi‐spacecraft magnetic field data on different spatial scales, directly comparing the results of three space missions: Time History of Events and Macroscale Interactions during Substorms (THEMIS), The Magnetospheric Multiscale (MMS) and Cluster, that have so far accrued decades of data on a range of spatial separation scales. We suggest that the statistical results of Cluster data over two decades can illustrate well the overall form of the ring current (in the ring plane), but can underestimate the current density, while THEMIS (partially) and MMS are able to observe the detailed structure of the ring current and indeed show different behavior on some (smaller) spatial scales. These results suggest that, in addition to the local time and radial form, reported previously, the morphology of the ring current may also depend on magnetic latitude, as influenced by the tilt angle between the inner and outer magnetic fields.

Tan, X. [School of Space and Earth Sciences Beihan↗
DOE OSTI2024

Dynamic embedding of effective harmonic normal mode vibrations in all-atomistic energy gap fluctuations: Case study of light harvesting 2 complex

Environmental effects in excitation energy transfer have mostly been modeled by baths of harmonic oscillators, but to what extent such modeling provides a reliable description of actual interactions between molecular systems and environments remains an open issue. Here, we address this issue by investigating fluctuations in the excitation energies of the light harvesting 2 complex using a realistic all-atomistic simulation of the potential energy surface. Our analyses reveal that molecular motions exhibit significant anharmonic features, even for underdamped intramolecular vibrations. In particular, we find that the anharmonicity contributes to the broadening of spectral densities and substantial overlaps between neighboring peaks, which complicates the meaning of mode frequencies constituting a bath model. Thus, we develop a strategy to construct a minimally underdamped harmonic bath that has a clear connection to all-atomistic dynamics by utilizing actual normal modes of molecules but optimizing their frequencies such that the resulting bath model can best reproduce the all-atomistic simulation results. By subtracting the underdamped contribution from the entire fluctuations, we also show that identifying a residual spectral density representing all other contributions with overdamped behavior is possible. We find that this can be fitted well with a well-established analytic form of a spectral density function or, alternatively, modeled as explicit time dependent fluctuations with muti-exponential or power law type correlation functions. We provide an assessment and the implications of these possibilities. The approach presented here can also serve as a general strategy to construct a simplified bath model that can effectively represent the underlying all-atomistic bath dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗