Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “neutral model”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

2025 ASMS Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS

Title (20 words): Investigation of the Collision-Induced Dissociation Mechanism of Protonated TODGA with IRIS Introduction (120 words): One of the challenges facing wide-spread adoption of nuclear power is the development of efficient separation processes for used nuclear fuel. The molecules in separation processes are subjected to an extreme environment due to the high radiation fields from the used fuel and highly acidic media used for fuel dissolution, which results in significant molecular degradation, leading to reduced process efficiency. These degradation products must be identified and studied so mitigation strategies can be developed to maintain process efficiency. However, complex systems can have many degradation products, complicating identification. Untargeted analysis tools could be used to understand radiation chemistry in complex systems. However, this would necessitate improved understanding of the gas-phase fragmentation mechanisms of fuel cycle molecules like tetraoctyldiglycolamide (TODGA). Methods (120 words): The gas-phase fragmentation of protonated TODGA was investigated using collision-induced dissociation (CID), resonance ejection, and infrared ion spectroscopy (IRIS). CID and resonance ejection experiments were conducted using a Bruker Daltonics (Bremen, Gemany) SolariX XR fourier transform ion cyclotron resonance (FT-ICR) mass spectrometer. IRIS spectra of protonated TODGA and its two CID fragmentation products were measured using a modified Bruker amaZon Speed ETD 3D quadrupole ion trap mass spectrometer coupled to the Free Electron Lasers for Infrared eXperiments (FELIX) free electron laser. Measured spectra were compared with density functional theory (DFT) calculations using the Gaussian 16, Revision C.02 software package with the ?B97X-D functional and def2-TZVPP basis sets. Candidate structures were generated using the CREST 3.0 conformational sampling software tool. Preliminary Data (300 words): Collision-induced dissociation of protonated TODGA ([C36H73N2O3]+, m/z=581.562) results two fragment ions, one at m/z=340.285 assigned as [C20H38NO3]+ and the other at m/z=312.290, assigned as [C19H38NO2]+. Based on the assigned formula and the structure of protonated TODGA, the fragment at m/z=340.285 is likely formed from elimination of neutral dioctylamine. Comparison of the IRIS spectrum of m/z=340.285 with DFT predictions suggests it contains a ring structure, and is assigned as N-octyl-N-(6-oxo-1,4-dioxan-2-ylidene)octan-1-aminium. Based on this structure and the structure of protonated TODGA, we hypothesize this fragment formed from elimination of neutral dioctylamine followed by a ring closure mechanism. Comparison of the IRIS spectrum of the fragment at m/z=312.290 with DFT predictions also indicated the presence of a ring structure, assigned as N-(1,3-dioxolan-4-ylidene)-N-octyloctan-1-aminium. This product could be formed from elimination of carbon monoxide from the ring of m/z=340.285 as a sequential fragmentation or formed directly from protonated TODGA via elimination of neutral N,N-dioctylformamide followed by a ring closure. Resonance ejection experiments where m/z=340 was continuously ejected from the IRC cell showed no decrease in intensity of m/z=312.290 across several collision energies, suggesting that the later, direct formation mechanism, dominates. The location of the ionizing proton in protonated TODGA is important for modeling the fragmentation mechanisms. DFT calculations suggested that the position of bands involving the coupled vibrations of the amide C—N and C=O bonds in TODGA are the most sensitive to proton location. Evaluation of the IRIS spectrum of protonated TODGA suggests that the ionizing proton is located between the two amid oxygens. This protonation location was calculated to lie approximately 30 kJ/mol lower in energy than the next lowest energy location, with the proton located solely on one of the amide oxygens. Novel aspect (20 words): Infrared ion spectroscopy combined with resonance ejection experiments and density functional theory to probe the collision-induced dissociation mechanism of tetraoctyldiglycolamide.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2021

Dynamical phase transitions in models of collective neutrino oscillations

Collective neutrino oscillations can potentially play an important role in transporting lepton flavor in astrophysical scenarios where the neutrino density is large, typical examples are the early universe and supernova explosions. It has been argued in the past that simple models of the neutrino Hamiltonian designed to describe forward scattering can support substantial flavor evolution on very short timescales t ≈ log(N) / (G Fρν ), with N the number of neutrinos, G F the Fermi constant and ρ ν the neutrino density. This finding is in tension with results for similar but exactly solvable models for which t ≈ √N/(G Fρν ) instead. In this work we provide a coherent explanation of this tension in terms of dynamical phase transitions (DPT) and study the possible impact that a DPT could have in more realistic models of neutrino oscillations and their mean-field approximation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

Evaluation of the Planetary Boundary Layer Height From ERA5 Reanalysis With MOSAiC Observations Over the Arctic Ocean

The planetary boundary layer height (PBLH) is a crucial indicator reflecting the region of the atmosphere characterized by continuous turbulence. Here, we use radiosonde and surface meteorological observations (4–7 times per day, year-round measurements) during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition to derive the PBLH (PBLH MOSAiC ), and further evaluate the PBLH from the ERA5 reanalysis (PBLH ERA5 ). Comparisons between PBLH MOSAiC and PBLH ERA5 from different perspectives reveal that: (a) The overestimation of PBLH ERA5 when the sea ice concentration is >90% is significant with the centered root mean squared error reaching up to 201 m; (b) The difference between the two products is notably pronounced in cold seasons, while it is comparatively diminished in warm seasons; (c) In neutral boundary layers, differences in PBLH ERA5 are larger compared with stable and convective boundary layers. In addition, the analysis of error sources indicates that the bias of PBLH ERA5 is sensitive to the bias of vertical thermal structure and wind speed profiles in ERA5 data sets in all conditions. Finally, we find a Random Forest model effectively reduces the bias of PBLH ERA5 with the index of agreement reaching up to 0.71 in the test data set, while a multiple linear regression demonstrates comparable performance to the Random Forest model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Impact of crowding on the diversity of expanding populations

Crowding effects critically impact the self-organization of densely packed cellular assemblies, such as biofilms, solid tumors, and developing tissues. When cells grow and divide, they push each other apart, remodeling the structure and extent of the population’s range. Recent work has shown that crowding has a strong impact on the strength of natural selection. However, the impact of crowding on neutral processes, which controls the fate of new variants as long as they are rare, remains unclear. Here, we quantify the genetic diversity of expanding microbial colonies and uncover signatures of crowding in the site frequency spectrum. By combining Luria–Delbrück fluctuation tests, lineage tracing in a novel microfluidic incubator, cell-based simulations, and theoretical modeling, we find that the majority of mutations arise behind the expanding frontier, giving rise to clones that are mechanically “pushed out” of the growing region by the proliferating cells in front. These excluded-volume interactions result in a clone-size distribution that solely depends on where the mutation first arose relative to the front and is characterized by a simple power law for low-frequency clones. Our model predicts that the distribution depends on a single parameter—the characteristic growth layer thickness—and hence allows estimation of the mutation rate in a variety of crowded cellular populations. Combined with previous studies on high-frequency mutations, our finding provides a unified picture of the genetic diversity in expanding populations over the whole frequency range and suggests a practical method to assess growth dynamics by sequencing populations across spatial scales.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Initial TCV operation with a baffled divertor

The Tokamak à Configuration Variable (TCV) tokamak is in the midst of an upgrade to further its capability to investigate conventional and alternative divertor configurations. To that end, modular and removable gas baffles have been installed to decrease the coupling between the divertor and the plasma core. The baffles primarily seek to suppress the transit of recycling neutrals to closed flux surfaces. A first experimental campaign with the gas baffles has shown that the baffled divertor remains compatible with a wide range of configurations including snowflake and super-X divertors. Plasma density ramp experiments reveal an increase of the neutral pressure in the divertor by up to a factor ×5 compared to the unbaffled divertor and thereby qualitatively confirm simulations with the SOLPS-ITER code that were used to guide the baffle design. Finally, together with a range of new and upgraded divertor diagnostics, the baffled TCV divertor is now used to validate divertor models for ITER and next step devices with particular emphasis on geometric variations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Photon classification with Gradient Boosted Trees at CLAS12

Dihadron semi-inclusive deep inelastic scattering (SIDIS) of 10.6 GeV longitudinally polarized electrons off the proton has been measured using the CLAS12 detector at Jefferson Lab. Two separate channels, π + π 0 and π - π 0 , were analyzed, requiring the reconstruction of diphoton pairs. Here, in this analysis, we addressed the problem of false neutral particles being reconstructed by CLAS12's event builder, polluting the otherwise physical combinatorial background underneath the π 0 peak. A photon classifier using a Gradient Boosted Trees (GBTs) architecture was trained with Monte Carlo simulations to reduce the amount of background π 0 's. We show that the nearest-neighbor features learned by the model lead to a substantial increase in signal vs. background discrimination compared to previous CLAS12 π^0 analyses. The machine learning approach recovers several times more dihadron statistics for the dataset.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Effects-Based Monitoring of Geomagnetically-Induced Current Using a Convolutional Neural Network

Geomagnetically-induced current (GIC) due to space weather can flow in the power grid causing undesirable effects such as transformer overheating, misoperation of protection devices, and potential blackouts. It is therefore important to monitor GIC in the power grid to improve online situational awareness and decision-making of system operators during a geomagnetic disturbance. To avoid the costly installation of GIC monitors at transformers’ neutrals, it is desirable to find correlations between GIC and already-monitored parameters. Hence, this work proposed the use of a convolutional neural network (CNN) to compute GIC amplitudes from learned patterns in the time-series data of transformer even harmonic currents. Using an electromagnetic transient program, GIC injection simulations were performed for a modeled Dominion Energy Virginia (DEV) substation with two 504 MVA, 500/230 kV transformers. Data collected from these offline simulations were used to train the CNN to provide online GIC monitoring. Testing the CNN performance involved using real GIC measurements from published literature and from a physical GIC monitor in the DEV area. Finally, the results showed that the proposed method was able to provide GIC readings with a root mean squared error of 1.56 A/phase (equivalent to an average accuracy of 94%) for these real GIC waveforms.

42 ENGINEERING↗
DOE OSTI2025

Nitrate-To-Ammonia Electroconversion at Neutral pH on Polycrystalline Vanadium Sulfide Derived from Vanadium Disulfide

The electrochemical nitrate reduction reaction (NO3RR) offers a pathway to produce NH3 for fuel and fertilizer from waste NO3-. In this work, a polycrystalline vanadium sulfide (VSx), which is derived from solvothermally grown and annealed VS2, is shown to exhibit excellent NO3RR activity (2.3 +- 0.6 mg.cm-2 geo..h-1 @ -0.92 VRHE) and Faradaic efficiency to NH4+ (69 +- 6% at -0.69 VRHE) in buffered neutral pH electrolyte containing 0.1 M NO3-. A variety of characterization techniques are leveraged to support the VSx assignment, including X-ray photoelectron spectroscopy, near-edge X-ray absorption fine structure spectroscopy, selected area electron diffraction, and X-ray diffraction measurements. The VS2 annealing step reduces the oxide character and generates VSx, which, based on the improved NO3RR activity, results in the creation of active sites for NO3- binding. To help shed light on NO3RR on VSx, VS2 is used as a model system, and a grand-canonical density functional theory (GC-DFT) investigation of VS2 shows strong evidence that S vacancies are active sites for NO3RR, where NO3- outcompetes H+ for adsorption at the S-vacancy sites. Moreover, GC-DFT results highlight a thermodynamically favorable reaction to generate NH4+ in an aqueous electrolyte at relevant cathodic potentials. As an annealed material, VSx may contain undersaturated V sites, which show an electronic structure similar to the theoretically calculated S-vacancy site of VS2, and these sites may contribute to the observed increase in NO3RR activity and selectivity for NH4+ on VSx versus unannealed VS2. Finally, kinetic isotope effect measurements suggest that the kinetic rate-limiting step of the NO3RR on VSx is not proton-coupled, indicating it may be the first electron transfer to adsorbed NO3*.

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

Overview of ASDEX upgrade results in view of ITER and DEMO

Experiments on ASDEX Upgrade (AUG) in 2021 and 2022 have addressed a number of critical issues for ITER and EU DEMO. A major objective of the AUG programme is to shed light on the underlying physics of confinement, stability, and plasma exhaust in order to allow reliable extrapolation of results obtained on present day machines to these reactor-grade devices. Concerning pedestal physics, the mitigation of edge localised modes (ELMs) using resonant magnetic perturbations (RMPs) was found to be consistent with a reduction of the linear peeling-ballooning stability threshold due to the helical deformation of the plasma. Conversely, ELM suppression by RMPs is ascribed to an increased pedestal transport that keeps the plasma away from this boundary. Candidates for this increased transport are locally enhanced turbulence and a locked magnetic island in the pedestal. The enhanced D-alpha (EDA) and quasi-continuous exhaust (QCE) regimes have been established as promising ELM-free scenarios. Here, the pressure gradient at the foot of the H-mode pedestal is reduced by a quasi-coherent mode, consistent with violation of the high-n ballooning mode stability limit there. This is suggestive that the EDA and QCE regimes have a common underlying physics origin. In the area of transport physics, full radius models for both L- and H-modes have been developed. These models predict energy confinement in AUG better than the commonly used global scaling laws, representing a large step towards the goal of predictive capability. A new momentum transport analysis framework has been developed that provides access to the intrinsic torque in the plasma core. In the field of exhaust, the X-Point Radiator (XPR), a cold and dense plasma region on closed flux surfaces close to the X-point, was described by an analytical model that provides an understanding of its formation as well as its stability, i.e., the conditions under which it transitions into a deleterious MARFE with the potential to result in a disruptive termination. With the XPR close to the divertor target, a new detached divertor concept, the compact radiative divertor, was developed. Here, the exhaust power is radiated before reaching the target, allowing close proximity of the X-point to the target. No limitations by the shallow field line angle due to the large flux expansion were observed, and sufficient compression of neutral density was demonstrated. With respect to the pumping of non-recycling impurities, the divertor enrichment was found to mainly depend on the ionisation energy of the impurity under consideration. In the area of MHD physics, analysis of the hot plasma core motion in sawtooth crashes showed good agreement with nonlinear 2-fluid simulations. This indicates that the fast reconnection observed in these events is adequately described including the pressure gradient and the electron inertia in the parallel Ohm’s law. Concerning disruption physics, a shattered pellet injection system was installed in collaboration with the ITER International Organisation. Thanks to the ability to vary the shard size distribution independently of the injection velocity, as well as its impurity admixture, it was possible to tailor the current quench rate, which is an important requirement for future large devices such as ITER. Progress was also made modelling the force reduction of VDEs induced by massive gas injection on AUG. The H-mode density limit was characterised in terms of safe operational space with a newly developed active feedback control method that allowed the stability boundary to be probed several times within a single discharge without inducing a disruptive termination. Regarding integrated operation scenarios, the role of density peaking in the confinement of the ITER baseline scenario (high plasma current) was clarified. The usual energy confinement scaling ITER98(p,y) does not capture this effect, but the more recent H20 scaling does, highlighting again the importance of developing adequate physics based models. Advanced tokamak scenarios, aiming at large non-inductive current fraction due to non-standard profiles of the safety factor in combination with high normalised plasma pressure were studied with a focus on their access conditions. A method to guide the approach of the targeted safety factor profiles was developed, and the conditions for achieving good confinement were clarified. Based on this, two types of advanced scenarios (‘hybrid’ and ‘elevated’

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Photoenzymatic Hydrogenation of Heteroaromatic Olefins Using ‘Ene’‐Reductases with Photoredox Catalysts

Abstract Flavin‐dependent ‘ene’‐reductases (EREDs) are highly selective catalysts for the asymmetric reduction of activated alkenes. This function is, however, limited to enones, enoates, and nitroalkenes using the native hydride transfer mechanism. Here we demonstrate that EREDs can reduce vinyl pyridines when irradiated with visible light in the presence of a photoredox catalyst. Experimental evidence suggests the reaction proceeds via a radical mechanism where the vinyl pyridine is reduced to the corresponding neutral benzylic radical in solution. DFT calculations reveal this radical to be “dynamically stable”, suggesting it is sufficiently long‐lived to diffuse into the enzyme active site for stereoselective hydrogen atom transfer. This reduction mechanism is distinct from the native one, highlighting the opportunity to expand the synthetic capabilities of existing enzyme platforms by exploiting new mechanistic models.

Nakano, Yuji↗
DOE OSTI2025

Unraveling the solvation phenomena of Ln 3+ cations in room-temperature ionic liquids: A computational study

Rare-earth elements (REEs), classified as critical materials, are difficult to separate due to their similar chemical properties and slight differences in ionic radius. Room-temperature ionic liquids (RTILs) have gained significant attention for REE separation because of their unique physicochemical properties and environmental advantages. In this study, we have investigated the solvation mechanisms of two RE cations (Ln 3+ ), Nd 3+ (light REE) and Yb 3+ (heavy REE), in two 1-butyl-3-methylimidazolium ([BMIM] + )-based RTILs using bis(trifluoromethylsulfonyl)imide ([NTf 2 ] − ) and acetate ([OAc] − ) as the respective anions, employing classical molecular dynamics (MD) simulations and density functional theory (DFT) techniques. From MD simulations, it is evident that Ln 3+ cations are primarily solvated by RTIL anions in the first solvation shell, while [BMIM] + cations form the second solvation shell with their numbers depending on the size and composition of the first solvation shell. The neutralizing [NO 3 − ] counterions are mainly solvated by [BMIM] + cations in their first solvation shell. Relative free energy change of solvation calculations using thermodynamic integration (TI) method indicate that Yb 3+ is more strongly solvated than Nd 3+ in both RTILs, a trend further supported by DFT calculations. While both methods predict consistent qualitative behavior, differences in models and energy evaluations lead to variations in absolute values. Overall, this study provides a comprehensive understanding of the solvation behavior of Ln 3+ cations in RTILs, demonstrating a stronger solvation preference for the heavy Ln 3+ cation (Yb 3+ ). In conclusion, these findings have implications for the design of RTIL-based separation processes for REEs.

Ash, Tamalika [Ames Lab., and Iowa State Univ., Am↗
DOE OSTI2025

Rarefied xenon flow in orificed hollow cathodes

A parametric study is conducted to quantify the effect of the keeper electrode geometry on the xenon neutral flow quantities within orificed hollow cathodes, prior to ignition. The keeper impinges directly on the flow out of the cathode orifice and its geometry influences the product between the pressure in the orifice–keeper region and the cathode-to-keeper distance. A representative cathode is simulated using the Direct Simulation Monte Carlo method. The numerical model is first validated with computational results from the literature. A parametric study is then conducted. Parameters include the cathode pressure–diameter in the range of 1–5 Torr cm and the following geometric ratios (and ranges): cathode orifice-to-inner radii (0.1–0.7), keeper orifice-to-cathode orifice radii (1–5), and keeper distance-to-cathode-orifice diameter (0.5–10). It is found that, if both keeper and cathode have identical orifice radii, the flow remains subsonic in the orifice-to-keeper region. In most cases, however, the flow becomes underexpanded and supersonic, and the static pressure within the orifice-to-keeper region is, on average, 4% that of the upstream pressure value. The orifice–keeper region pressure increases with either a decrease in the keeper orifice diameter or an increase in the distance between cathode and keeper, in agreement with literature data. Both trends are explained through conservation laws. A statistical study of numerical results reveals that the ratio of ignition-to-nominal mass flow rates has a most probable value of 50, which suggests that heaterless cathode ignition at a minimum DC voltage may be achieved by increasing the input mass flow rate by a factor of 50.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Predicting soil carbon changes in switchgrass grown on marginal lands under climate change and adaptation strategies

Abstract The United States Great Lakes Region (USGLR) is a critical geographic area for future bioenergy production. Switchgrass ( Panicum virgatum ) is widely considered a carbon (C)‐neutral or C‐negative bioenergy production system, but projected increases in air temperature and precipitation due to climate change might substantially alter soil organic C (SOC) dynamics and storage in soils. This study examined long‐term SOC changes in switchgrass grown on marginal land in the USGLR under current and projected climate, predicted using a process‐based model (Systems Approach to Land‐Use Sustainability) extensively calibrated with a wealth of plant and soil measurements at nine experimental sites. Simulations indicate that these soils are likely a net C sink under switchgrass (average gain 0.87 Mg C ha −1 year −1 ), although substantial variation in the rate of SOC accumulation was predicted (range: 0.2–1.3 Mg C ha −1 year −1 ). Principal component analysis revealed that the predicted intersite variability in SOC sequestration was related in part to differences in climatic characteristics, and to a lesser extent, to heterogeneous soils. Although climate change impacts on switchgrass plant growth were predicted to be small (4%–6% decrease on average), the increased soil respiration was predicted to partially negate SOC accumulations down to 70% below historical rates in the most extreme scenarios. Increasing N fertilizer rate and decreasing harvest intensity both had modest SOC sequestration benefits under projected climate, whereas introducing genotypes better adapted to the longer growing seasons was a much more effective strategy. Best‐performing adaptation scenarios were able to offset >60% of the climate change impacts, leading to SOC sequestration 0.7 Mg C ha −1 year −1 under projected climate. On average, this was 0.3 Mg C ha −1 year −1 more C sequestered than the no adaptation baseline. These findings provide crucial knowledge needed to guide policy and operational management for maximizing SOC sequestration of future bioenergy production on marginal lands in the USGLR.

09 BIOMASS FUELS↗
DOE OSTI2021

Dynamic life-cycle carbon analysis for fast pyrolysis biofuel produced from pine residues: implications of carbon temporal effects

Abstract Background Woody biomass has been considered as a promising feedstock for biofuel production via thermochemical conversion technologies such as fast pyrolysis. Extensive Life Cycle Assessment studies have been completed to evaluate the carbon intensity of woody biomass-derived biofuels via fast pyrolysis. However, most studies assumed that woody biomass such as forest residues is a carbon–neutral feedstock like annual crops, despite a distinctive timeframe it takes to grow woody biomass. Besides, few studies have investigated the impacts of forest dynamics and the temporal effects of carbon on the overall carbon intensity of woody-derived biofuels. This study addressed such gaps by developing a life-cycle carbon analysis framework integrating dynamic modeling for forest and biorefinery systems with a time-based discounted Global Warming Potential (GWP) method developed in this work. The framework analyzed dynamic carbon and energy flows of a supply chain for biofuel production from pine residues via fast pyrolysis. Results The mean carbon intensity of biofuel given by Monte Carlo simulation across three pine growth cases ranges from 40.8–41.2 g CO 2 e MJ −1 (static method) to 51.0–65.2 g CO 2 e MJ −1 (using the time-based discounted GWP method) when combusting biochar for energy recovery. If biochar is utilized as soil amendment, the carbon intensity reduces to 19.0–19.7 g CO 2 e MJ −1 (static method) and 29.6–43.4 g CO 2 e MJ −1 in the time-based method. Forest growth and yields (controlled by forest management strategies) show more significant impacts on biofuel carbon intensity when the temporal effect of carbon is taken into consideration. Variation in forest operations and management (e.g., energy consumption of thinning and harvesting), on the other hand, has little impact on the biofuel carbon intensity. Conclusions The carbon temporal effect, particularly the time lag of carbon sequestration during pine growth, has direct impacts on the carbon intensity of biofuels produced from pine residues from a stand-level pine growth and management point of view. The carbon implications are also significantly impacted by the assumptions of biochar end-of-life cases and forest management strategies.

09 BIOMASS FUELS↗
DOE OSTI2021

Multiwavelength Observation Campaign of the TeV Gamma-Ray Binary HESS J0632 + 057 with NuSTAR, VERITAS, MDM, and Swift

HESS J0632+057 belongs to a rare subclass of binary systems that emit gamma rays above 100 GeV. It stands out for its distinctive high-energy light curve, which features a sharp "primary" peak and broader "secondary" peak. We present the results of contemporaneous observations by NuSTAR and VERITAS during the secondary peak between 2019 December and 2020 February, when the orbital phase (Φ) is between 0.55 and 0.75. NuSTAR detected X-ray spectral evolution, while VERITAS detected TeV emission. We fit a leptonic wind-collision model to the multiwavelength spectra data obtained over the four NuSTAR and VERITAS observations, constraining the pulsar spin-down luminosity and the magnetization parameter at the shock. Despite long-term monitoring of the source from 2019 October to 2020 March, the MDM observatory did not detect significant variation in Hα and Hβ line equivalent widths, an expected signature of Be-disk interaction with the pulsar. Furthermore, fitting folded Swift-XRT light-curve data with an intrabinary shock model constrained the orbital parameters, suggesting two orbital phases (at Φ D = 0.13 and 0.37), where the pulsar crosses the Be-disk, as well as phases for the periastron (Φ 0 = 0.30) and inferior conjunction (Φ IFC = 0.75). The broadband X-ray spectra with Swift-XRT and NuSTAR allowed us to measure a higher neutral hydrogen column density at one of the predicted disk-passing phases.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Symmetrized Drude Oscillator Force Fields Improve Numerical Performance of Polarizable Molecular Dynamics

Drude oscillator potentials are a popular and computationally efficient class of polarizable models that represent each polarizable atom as a positively charged Drude core harmonically bound to a negatively charged Drude shell. Here we show that existing force fields that place all non-Coulomb forces on the Drude core and none on the shell inadvertently couple the dipole to non-Coulombic forces. This introduces errors where interactions with neutral particles can erroneously induce atomic polarization, leading to spurious polarizations in the absence of an electric field, exacerbating violations of equipartition in the employed Carr–Parinello scheme. A suitable symmetrization of the interaction potential that correctly splits the force between the Drude core and shell can correct this shortcoming, improving the stability and numerical performance of Drude oscillator-based simulations. The symmetrization procedure is straightforward and only requires the rescaling of a few force field parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗