Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sampling methods”

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

Protocol for single-cell isolation and genome amplification of environmental microbial eukaryotes for genomic analysis

We describe environmental microbial eukaryotes (EMEs) sample collection, single-cell isolation, lysis, and genome amplification, followed by the rDNA amplification and OTU screening for recovery of high-quality species-specific genomes for de novo assembly. These protocols are part of our pipeline that also includes bioinformatic methods. The pipeline and its application on a wide range of phyla of different sample complexity are described in our complementary paper. In addition, this protocol describes optimized lysis, genome amplification, and OTU screening steps of the pipeline. For complete details on the use and execution of this protocol, please refer to Ciobanu et al. (2021).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Separation of Lipoproteins for Quantitative Analysis of 14 C-Labeled Lipid-Soluble Compounds by Accelerator Mass Spectrometry

To date, 14 C tracer studies using accelerator mass spectrometry (AMS) have not yet resolved lipid-soluble analytes into individual lipoprotein density subclasses. The objective of this work was to develop a reliable method for lipoprotein separation and quantitative recovery for biokinetic modeling purposes. The novel method developed provides the means for use of small volumes (10–200 µL) of frozen plasma as a starting material for continuous isopycnic lipoprotein separation within a carbon- and pH-stable analyte matrix, which, following post-separation fraction clean up, created samples suitable for highly accurate 14 C/ 12 C isotope ratio determinations by AMS. Manual aspiration achieved 99.2 ± 0.41% recovery of [5- 14 CH 3 ]-(2R, 4'R, 8'R)-α-tocopherol contained within 25 µL plasma recovered in triacylglycerol rich lipoproteins (TRL = Chylomicrons + VLDL), LDL, HDL, and infranatant (INF) from each of 10 different sampling times for one male and one female subject, n = 20 total samples. Small sample volumes of previously frozen plasma and high analyte recoveries make this an attractive method for AMS studies using newer, smaller footprint AMS equipment to develop genuine tracer analyses of lipophilic nutrients or compounds in all human age ranges.

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

Temporal Study 2021-2022: Sensor-Based Time Series of Surface Water Temperature, Specific Conductance, Total Dissolved Solids, Turbidity, pH, and Dissolved Oxygen from across Multiple Watersheds in the Yakima River Basin in Washington, USA (v2)

This dataset supports a broader study examining the drivers of temporal variability in sediment respiration rates in the Yakima River Basin. The dataset provides periodic (weekly or biweekly) in situ hydrological and water chemistry sensor data, handheld sensor water chemistry data, general environmental context photos, and field metadata collected at six sites within multiple rivers in the Yakima River Basin in Washington, USA. In addition to the sensor data, there are plots of continuous in situ sensor data and R scripts used to generate the plots. Related sample-based water chemistry data are published separately at https://data.ess-dive.lbl.gov/datasets/doi:10.15485/1898912.The data package was originally published in September 2022. It was updated in June 2025 (v2; modified files). See the change history section in the readme for more details. For details on how to navigate data packages generated by this project, see https://data.ess-dive.lbl.gov/portals/PNNLRiverCorridorSFA/About. This dataset is comprised of one main data folder containing two sensor-specific subfolders, and one photographs folder. The main data folder includes file-level metadata (flmd), data dictionary (dd), installation methods, field metadata, handheld sensor data, field data collection protocols, international generic sample number (IGSN) mapping file, and a readme file. Each sensor subfolder (BarotrollAtm and MantaRiverData) contains a subfolder containing sensor timeseries data and plots. The BarotrollAtm Data subfolder contains In Situ Rugged BaroTROLL sensor pressure and air temperature data. The MantaRiverData subfolder contains Eureka Manta+ 35B multisonde temperature, specific conductance, and turbidity. The FieldPhotos folder contains environmental context photographs and videos. All files are .csv, .pdf, .R, .jpg, .jpeg, .heic, .mov, or .mp4.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Simultaneous Determination of Halogens and Metals in Waste Plastic Pyrolysis Oil by Inductively Coupled Plasma Mass Spectrometry

A major barrier to integrating pyrolysis-derived oil into conventional refinery technology is the presence of impurities, particularly halogens and metals, that can deactivate catalysts. This study presents a novel, cost-effective approach for the simultaneous analysis of a subset of halogens, metals, nonmetals, and metalloids in complex, industrially relevant distilled pyrolysis oil samples, this was previously achievable only through multiple techniques. Excellent results were obtained using a widely accessible inductively coupled plasma mass spectrometry (ICP-MS) method with helium gas mode and standard laboratory consumables, enabling high-throughput analysis and efficient evaluation of adsorbent performance. Sample preparation is straightforward, requiring only dilution in a compatible matrix, and provides accurate and precise quantification, with 75%–137% spike recovery for Be, Ti, V, Cr, Fe, Ni, Co, Cu, As, Se, Mo, Cd, Sb, Tl, and Pb, and 62%–65% spike recovery for Cl and 109%–133% spike recovery for Br. Additionally, an enhanced version of the method using ICP-MS/MS and hydrogen gas is described, which has higher accuracy with 69%–112% spike recovery for Be, B, Ti, V, Fe, Co, Ni, Cu, As, Se, Mo, Cd, Sb, Tl, and Pb, with 85%–102% spike recovery for Cl and 63%–93% spike recovery for Br. Helium mode detection limits for industrially relevant elements (V, Fe, Ni, Cu, As, and Pb) are less than 1.7 µg/kg, and less than 0.2 mg/kg for Br and Cl. Furthermore, this methodology facilitates rapid systematic evaluation of adsorption capacities of materials under time-on-stream conditions and supports robust comparisons across diverse operating environments.

Lazarcik, James [University of Wisconsin-Madison, ↗
DOE OSTI2023

Yakima River Basin Temporal Study: Sensor and Sample Data from Wenas Creek following the Evans Canyon Fire in Washington, USA

The Evans Canyon Fire occurred in Washington in August 2020 and burned 76,000 acres of the semi-arid shrub-steppe landscape at low to moderate severities. This fire burned across the Wenas Creek watershed, allowing for a sampling design that included a portion of the stream within the burned area and a reference site upstream of the burn. Monthly data was collected at an unburned (W10) and burned (W20) site from November 2020 to August 2022 for in situ hydrological and water chemistry sensor data, handheld sensor water chemistry data, general environmental context photos, and field metadata. This dataset is comprised of one folder with field photos and one main data folder containing (1) sensor and sample field protocols; (2) file-level metadata (flmd); (3) data dictionary (dd); (4) international geo-sample number (IGSN) mapping file; (5) field metadata; (6) readme; (7) methods codes; (8) dissolved organic carbon (DOC; measured as non-purgeable organic carbon; NPOC); (9) total nitrogen (TN); (10) total suspended solids (TSS); (11) sensor data (specific conductance, pH, total dissolved solids, temperature, pressure, dissolved oxygen, and turbidity) averages; and (12) sensor installation methods. All files are .csv, .jpg, .jpeg, or .pdf.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Determination of acoustic nonlinearity parameters using thermal modulation of ultrasonic waves

This study presents a test method and its theoretical framework to determine the acoustic nonlinearity parameters (α, β, δ) of material using thermal modulation of ultrasonic waves. Temperature change-induced thermal strain excites the nonlinear response of the material and modulates the ultrasonic wave propagating in it. Experimental results showed a strong correlation between the relative wave velocity change and the temperature change. With a quadratic polynomial model, the acoustic nonlinearity parameters were obtained from the polynomial coefficients by curve fitting the experimental data. Their effects on thermal-induced velocity change were discussed. The parameters α, β, and δ govern the hysteretic gap, average slope, and curvature of the correlation curve, respectively. The proposed theory was validated on aluminum, steel, intact and damaged concrete samples. The obtained nonlinear parameters show reasonable agreement with values reported in the literature. Compared to other nonlinear acoustic methods using vibration or acoustic excitation, the thermal modulation method generates more uniform, slow changing, and larger strain field in the test sample. Employing the thermal effect as the driving force for nonlinearity instead of an undesired influencing factor, this method can measure the absolute values of α, β, and δ with good accuracy using a simple ultrasonic test setup.

Physics↗
DOE OSTI2021

Pulsar Based Timing for Grid Synchronization

Existing synchronization systems in the power grid, such as the global positioning system, are susceptible to temporary or permanent failures due to various unpredictable and uncontrollable factors such as cyber-attack and electromagnetic interferences, thus affecting the accuracy and reliability of generated timing signal. In this article, a pulsar astronomy-based timing system is proposed to provide an alternative synchronization signal. Further, this clock will offer significant security improvements to power grid applications, such as a wide-area monitoring system, which depends on a precise timing signal. The hardware and software frameworks are described in detail. First, a high-speed sampling hardware platform is designed to collect signals from radio telescopes. Then a periodic pulse extraction method with three steps is proposed to process the pulsar signal, including polyphase filterbanks, incoherent de-dispersion, and sliding window folding. Lastly, three experiments are conducted to verify the effectiveness of the frameworks. The generated pulsar timing pulse is presented, and the factors affecting its accuracy are also discussed. The analysis results demonstrate that the pulsar signals can provide high-accurate timing pulses for grid synchronization.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Challenges in studying water fluxes within the soil-plant-atmosphere continuum: A tracer-based perspective on pathways to progress

Tracing and quantifying water fluxes in the hydrological cycle is crucial for understanding the current state of ecohydrological systems and their vulnerability to environmental change. Especially the interface between ecosystems and the atmosphere that is strongly mediated by plants is important to meaningfully describe ecohydrological system functioning. Many of the dynamic interactions generated by water fluxes between soil, plant and the atmosphere are not well understood, which is partly due to a lack of interdisciplinary research. This opinion paper reflects the outcome of a discussion among hydrologists, plant ecophysiologists and soil scientists on open questions and new opportunities for collaborative research on the topic "water fluxes in the soil-plant-atmosphere continuum" especially focusing on environmental and artificial tracers. We emphasize the need for a multi-scale experimental approach, where a hypothesis is tested at multiple spatial scales and under diverse environmental conditions to better describe the small-scale processes (i.e., causes) that lead to large-scale patterns of ecosystem functioning (i.e., consequences). Novel in-situ, high-frequency measurement techniques offer the opportunity to sample data at a high spatial and temporal resolution needed to understand the underlying processes. Here we advocate for a combination of long-term natural abundance measurements and event-based approaches. Multiple environmental and artificial tracers, such as stable isotopes, and a suite of experimental and analytical approaches should be combined to complement information gained by different methods. Virtual experiments using process-based models should be used to inform sampling campaigns and field experiments, e.g., to improve experimental designs and to simulate experimental outcomes. On the other hand, experimental data are a pre-requisite to improve our currently incomplete models. Interdisciplinary collaboration will help to overcome research gaps that overlap across different earth system science fields and help to generate a more holistic view of water fluxes between soil, plant and atmosphere in diverse ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Preventing Bio-Bloopers and XFEL Follies: Best Practices from your Friendly Instrument Staff

Serial Femtosecond Crystallography (SFX) at X-ray Free electron Lasers (XFELs) is a relatively new field promising to deliver unparalleled spatial and temporal resolution on biological systems and there dynamics. Over the past decade, though, there have been a handful of results that have truly delivered on these promises. Why? SFX has many paradigm shifting techniques not seen in typical structural biology arenas, such as creating a concentrated slurry of microcrystals rather than a handful of loopable prisms worthy of a catalog photo. Then taking that slurry and high speed jetting them towards the vacuum or helium interation region to destroy less than 1% of your sample and waste the other 99. The literature is full of techniques and methods promising to cure what ails your experiment, yet as an instrument scientist will tell you –and a first author might admit after a few drinks at the conference happy hour—is that there are a lot more failures than the success we published, results may vary. We will walk through a best practices on how to prepare your sample and chose a sample delivery technique that will amerliorate your studies rather than undermine your hardwork and hopefully lead to better experimental planning and execution, inching you closer to that scientific goal and that call from Stockholm. This will be written in a more editorialized fashion and is meant to give the reader an idea of what to try or how they should be thinking. Welcome to SFX, now what?

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

NIR Dataset and Models (Near-Infrared Spectroscopy Calibration Dataset for Butanediol Fermentation and Multivariate Calibration Models) [SWR-22-60]

2,3-butanediol (2,3-BDO) is an economically important platform chemical that can be used in a variety of chemical feedstocks, liquid fuels, and biosynthetic building blocks. While 2,3-BDO can be efficiently produced by fermentation, the fermentation requires continuous monitoring and control to maximize 2,3-BDO yields and minimize inhibitory coproducts. Because of the time required for sampling and at-line measurement of fermentation samples with high pressure liquid chromatography (HPLC), the ability for operators to perform real-time modification to fermentation conditions is limited. To overcome this challenge, researchers from the National Renewable Energy Laboratory (NREL) have developed a calibration model which can predict the concentration of several analytes in real-time using near-infrared (NIR) spectra of the filtered fermentation broth. While significantly reducing the need for off-line sampling, NREL expects this technology to play a critical role in maximizing 2,3-BDO production. 2,3-butanediol (2,3-BDO) is a useful chemical platform that can be used to create a variety of products. For instance, 2,3-BDO can be (1) dehydrated and converted into methyl ethyl ketone, a liquid fuel additive or (2) deoxydehydrated into 1,3-butadiene for synthetic rubber, which can also be oligomerized in high yields to gasoline, diesel, and jet fuel. In order to maximize 2,3-BDO production, frequent measurement of fermentation samples is needed, as small changes in oxygen concentration can drive the fermentation to undesired products. For example, oxygen-deficient conditions result in glycerol production, while excess oxygen concentrations result in acetoin production. This results in the need for measuring dissolved oxygen, glucose, and xylose concentrations in order to optimize the aeration rate of the fermentation. Traditional monitoring methods occurs off-line and can take up to 30 minutes per sample. With multiple fermenters and high-pressure liquid chromatography (HPLC) injectors, resulting in the need for multiple samples, the sampling process can take hours to complete. NREL’s calibration model can predict the glucose, xylose, 2,3-BDO, acetoin, and glycerol concentrations from NIR spectra of filtered fermentation liquor samples. Using a partial least-squares (PLS) calibration model, NREL’s model can monitor the concentration of these analytes during subsequent fermentations at bench- and pilot-scale, demonstrating the utility of NIR spectroscopy combined with chemometrics for real-time, at-line monitoring of 2,3-BDO fermentations.

Wolfrum, Edward↗
DOE OSTI2023

Comprehensive analysis of common polymers using hyphenated TGA-FTIR-GC/MS and Raman spectroscopy towards a database for micro- and nanoplastics identification, characterization, and quantitation

Environmental contamination by micro- and nanoplastics (MNPs) is well documented with potential for their increased accumulation globally. Growing public concern over environmental, ecological, and human exposure to MNPs has led to exponential increase in publications, news articles, and reports. Significant knowledge gap exists in standardized analytical methods for the identification and quantification of MNPs from real world environmental samples. Here, in this study, we report comprehensive datasets utilizing thermogravimetric analyzer (TGA) coupled to a Fourier transformed infrared spectrometer (FTIR) and a gas chromatography/mass spectrometer (GC/MS) with corresponding Raman spectral data for the most common polymers documented to be present in the environment (35 plastics of 12 polymer types), to serve as a base line reference for the identification and quantitation of MNPs. Various parameters for TGA-FTIR-GC/MS data acquisition were optimized. Commercial consumer plastic product compositions were identified using this analytical database. Case studies to showcase the utility of the method for polymer mixtures analysis is included. This dataset would serve towards the development of a collaborative, global, comprehensive, and curated public database for the identification of various MNPs and mixtures.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2023

Predicting rare events using neural networks and short-trajectory data

Estimating the likelihood, timing, and nature of events is a major goal of modeling stochastic dynamical systems. When the event is rare in comparison with the timescales of simulation and/or measurement needed to resolve the elemental dynamics, accurate prediction from direct observations becomes challenging. In such cases a more effective approach is to cast statistics of interest as solutions to Feynman-Kac equations (partial differential equations). Here, we develop an approach to solve Feynman-Kac equations by training neural networks on short-trajectory data. Our approach is based on a Markov approximation but otherwise avoids assumptions about the underlying model and dynamics. This makes it applicable to treating complex computational models and observational data. Additionally, we illustrate the advantages of our method using a low-dimensional model that facilitates visualization, and this analysis motivates an adaptive sampling strategy that allows on-the-fly identification of and addition of data to regions important for predicting the statistics of interest. Finally, we demonstrate that we can compute accurate statistics for a 75-dimensional model of sudden stratospheric warming. This system provides a stringent test bed for our method.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

The information content of projected galaxy fields

ABSTRACT The power spectrum of the non-linearly evolved large-scale mass distribution recovers only a minority of the information available on the mass fluctuation amplitude. We investigate the recovery of this information in 2D ‘slabs’ of the mass distribution averaged over ≈100 h−1 Mpc along the line of sight, as might be obtained from photometric redshift surveys. We demonstrate a Hamiltonian Monte Carlo method to reconstruct the non-Gaussian mass distribution in slabs, under the assumption that the projected field is a point-transformed Gaussian random field, Poisson-sampled by galaxies. When applied to the Quijote N-body suite at z = 0.5 and at a transverse resolution of 2 h−1 Mpc, the method recovers ∼30 times more information than the 2D power spectrum in the well-sampled limit, recovering the Gaussian limit on information. At a more realistic galaxy sampling density of 0.01 h3 Mpc−3, shot noise reduces the information gain to a factor of 5 improvement over the power spectrum at resolutions of 4 h−1 Mpc or smaller.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Barium stars as tracers of s -process nucleosynthesis in AGB stars

Barium (Ba) stars help to verify asymptotic giant branch (AGB) star nucleosynthesis models since they experienced pollution from an AGB binary companion and thus their spectra carry the signatures of the slow neutron capture process (s process). For a large number (180) of Ba stars, we searched for AGB stellar models that match the observed abundance patterns. We aim to uncover any systematic deviations of the sample abundances from the predictions of the nucleosynthesis models. We employed three machine learning algorithms as classifiers: a Random Forest method, developed for this work, and the two classifiers used in our previous study. Compared to that work, we also expanded our observational sample with 11 Ba stars available in the supersolar metallicity range. We studied the statistical behaviour of the different s-process elements in the observational sample to investigate if the AGB models systematically under- or overpredict the abundances observed in the Ba stars and show the results in the form of violin plots of the residuals between spectroscopic abundances and model predictions. We inspected the correlations between the observed [Fe/H], the s-process elemental abundances, and the residuals. We employed the [Zr/Fe] and [Nb/Fe] abundances as a thermometer to constrain the operational temperature that rules the production of these elements in the sample stars, assuming a steady-state s process. We also investigated the mass distribution of the identified polluter AGB stars and the behaviour of the δ parameter, which describes the fraction of accreted AGB material relative to the Ba star envelope. We find a significant trend in the residuals that implies an underproduction of the elements just after the first s-process peak (Nb, Mo, and Ru) in the models relative to the observations. This may originate from a neutron-capture process (e.g. the intermediate neutron-capture process, i process) not yet included in the AGB models of metallicity from solar to roughly 1/5 solar, corresponding to the range of the Ba stars. Correlations are found between the residuals of these peculiar elements, suggesting a common origin for the deviations from the models. In addition, there is a weak metallicity dependence of the residuals of these elements. The s-process temperatures derived with the [Zr/Fe] – [Nb/Fe] thermometer have an unrealistic value for the majority of our stars. The most likely explanation is that at least a fraction of these elements are not produced in a steady-state s process, and instead may be due to processes not included in the AGB models. The mass distribution of the identified models confirms that our sample of Ba stars was polluted by low-mass AGB stars (< 4 M ⊙ ). Most of the matching AGB models require low accreted mass, but a few systems with high accreted mass are needed to explain the observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Machine Learning for First Principles Calculations of Material Properties for Ferromagnetic Materials

The investigation of finite temperature properties using Monte-Carlo (MC) methods requires a large number of evaluations of the system’s Hamiltonian to sample the phase space needed to obtain physical observables as function of temperature. DFT calculations can provide accurate evaluations of the energies, but they are too computationally expensive for routine simulations. To circumvent this problem, machine-learning (ML) based surrogate models have been developed and implemented on high-performance computing (HPC) architectures. In this paper, we describe two ML methods (linear mixing model and HydraGNN) as surrogates for first principles density functional theory (DFT) calculations with classical MC simulations. These two surrogate models are used to learn the dependence of target physical properties from complex compositions and interactions of their constituents. We present the predictive performance of these two surrogate models with respect to their complexity while avoiding the danger of overfitting the model. An important aspect of our approach is the periodic retraining with newly generated first principles data based on the progressive exploration of the system’s phase space by the MC simulation. The numerical results show that HydraGNN model attains superior predictive performance compared to the linear mixing model for magnetic alloy materials.

Eisenbach, Markus↗
DOE OSTI2026

Determination of proton PDF uncertainties with Markov chain Monte Carlo

We present an analysis of parton distribution functions (PDFs) of the proton using Markov chain Monte Carlo (MCMC) methods. The MCMC approach naturally implements Bayes’ theorem and, thus, provides a means to directly sample the underlying probability distribution—in this case, the probability distribution of the PDF parameters. This allows for a straightforward propagation of the resulting uncertainties into any PDF-dependent observable, preserving their simple probabilistic interpretation. In our analysis we include a broad set of deep inelastic scattering data from HERA, BCDMS and NMC experiments along with the Drell-Yan, 𝑊 and 𝑍 boson data from LHC and Tevatron experiments, which combined with theoretical calculations at next-to-next-to-leading order in QCD allow for realistic determination of PDFs. The main focus of this analysis is to explore alternative methods for PDF uncertainty estimation that are more firmly grounded in statistical principles. We show that the flexibility of the Bayes framework, allowing one, e.g., to account for non-Gaussianity or inconsistencies of datasets, is crucial to extract realistic uncertainties when such assumptions are not fulfilled. We also demonstrate that MCMC allows one to determine the Δ⁢𝜒 2 value corresponding to a given confidence level in the sample, which can, in turn, be used as a statistically well-founded tolerance criterion used in the Hessian method, thus addressing one of its main long-standing drawbacks.

Risse, Peter Clemens [Universität Münster (Germany↗
DOE OSTI2025

Emerging Per- and Polyfluoroalkyl Substances in Tap Water from the American Healthy Homes Survey II

Humans experience widespread exposure to anthropogenic per- and polyfluoroalkyl substances (PFAS) through various media, which can lead to a wide range of negative health impacts. Tap water is an important source of exposure in communities with any degree of contamination but routine or large-scale PFAS monitoring often depends on targeted analytical methods limited to measuring specific PFAS. We analyzed 680 tap water samples from the American Healthy Homes Survey II for PFAS using non-targeted analysis (NTA) to expand the range of detectable PFAS. Based on detection frequency and relative abundance, about half of the identified PFAS were found only by NTA. We identified (with varying degrees of confidence) 75 distinct PFAS, including 57 exclusively detected by NTA. The identified PFAS are members of seven structural subclasses differentiated by their head groups and degree of fluorination. Clustering analysis categorized the PFAS into four coabundance groups dominated by specific PFAS subclasses. One group uniquely identified by NTA contains zwitterionic PFAS and other PFAS transformation products which are likely associated with aqueous firefighting foam contaminants in a small number of spatially correlated samples. These results help further characterize the scope of exposure to emerging PFAS experienced by the U.S. population via tap water and augment nationwide targeted-PFAS monitoring programs.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2021

Bringing discrete-time Langevin splitting methods into agreement with thermodynamics

In light of the recently published complete set of statistically correct Grønbech–Jensen (GJ) methods for discrete-time thermodynamics, we revise a differential operator splitting method for the Langevin equation in order to comply with the basic GJ thermodynamic sampling features, namely, the Boltzmann distribution and Einstein diffusion, in linear systems. This revision, which is based on the introduction of time scaling along with flexibility of a discrete-time velocity attenuation parameter, provides a direct link between the ABO splitting formalism and the GJ methods. This link brings about the conclusion that any GJ method has at least weak second order accuracy in the applied time step. It further helps identify a novel half-step velocity, which simultaneously produces both correct kinetic statistics and correct transport measures for any of the statistically sound GJ methods. Explicit algorithmic expressions are given for the integration of the new half-step velocity into the GJ set of methods. Finally, numerical simulations, including quantum-based molecular dynamics (QMD) using the QMD suite Los Alamos Transferable Tight-Binding for Energetics, highlight the discussed properties of the algorithms as well as exhibit the direct application of robust, time-step-independent stochastic integrators to QMD.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Refine your search to explore more results.