Search NASA⌕ Search

SEARCH · Search NASA

Results for “Rescaling”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Proton decay from quark and lepton compositeness

Within a chiral SU(15) gauge theory in which the Standard Model fermions are bound states of massless preons, we show that proton-decay operators are likely induced at the compositeness scale, Λ pre . Our estimate of the limit imposed by searches for proton decays is Λ pre > 10 4 TeV 1/2 $C^{1/4}_{8} m^{1/2}_{Q}$, where C 8 is a rescaled coefficient of an 8-prebaryon operator induced by SU(15) dynamics, and m Q is the mass of a composite vectorlike quark. The latter has a lower limit related to the mass of a composite vectorlike lepton, which in turn is required by LHC searches to be above 1 TeV. For C 8 in the 10 –5 – 1 range, the lower limit on Λ pre varies between 3 × 10 3 TeV and 5 × 10 4 TeV. We point out that exotic proton decay modes, into a π + and a heavy right-handed neutrino, could be observed using the Super-Kamiokande or DUNE detectors.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Discrete Superconvergence Analysis for Quantum Magnus Algorithms of Unbounded Hamiltonian Simulation

Motivated by various applications, unbounded Hamiltonian simulation has recently garnered great attention. Quantum Magnus algorithms, designed to achieve commutator scaling for time-dependent Hamiltonian simulation, have been found to be particularly efficient for such applications. When applied to unbounded Hamiltonian simulation in the interaction picture, they exhibit an unexpected superconvergence phenomenon. However, existing proofs are limited to the spatially continuous setting and do not extend to discrete spatial discretizations. Here, in this work, we provide the first superconvergence estimate in the fully discrete setting with a finite number of spatial discretization points N, and show that it holds with an error constant uniform in N. The proof is based on the two-parameter symbol class, which, to our knowledge, is applied for the first time in algorithm analysis. The key idea is to establish a semiclassical framework by identifying two parameters through the discretization number and the time step size rescaled by the operator norm, such that the semiclassical uniformity guarantees the uniformity of both. This approach may have broader applications in numerical analysis beyond the specific context of this work.

Borns-Weil, Yonah [University of California, Berke↗
DOE OSTI2025

Optimizing the optimizer for physics-informed neural networks and Kolmogorov-Arnold networks

Physics-Informed Neural Networks (PINNs) have revolutionized the computation of PDE solutions by integrating partial differential equations (PDEs) into the neural network’s training process as soft constraints, becoming an important component of the scientific machine learning (SciML) ecosystem. More recently, physics-informed Kolmogorv-Arnold networks (PIKANs) have also shown to be effective and comparable in accuracy with PINNs. In their current implementation, both PINNs and PIKANs are mainly optimized using first-order methods like Adam, as well as quasi-Newton methods such as BFGS and its low-memory variant, L-BFGS. However, these optimizers often struggle with highly nonlinear and non-convex loss landscapes, leading to challenges such as slow convergence, local minima entrapment, and (non)degenerate saddle points. In this study, we investigate the performance of Self- Scaled BFGS (SSBFGS), Self-Scaled Broyden (SSBroyden) methods and other advanced quasi-Newton schemes, including BFGS and L-BFGS with different line search strategies. These methods dynamically rescale updates based on historical gradient information, thus enhancing training efficiency and accuracy. We systematically compare these optimizers – using both PINNs and PIKANs – on key challenging PDEs, including the Burgers, Allen-Cahn, Kuramoto-Sivashinsky, Ginzburg-Landau, and Stokes equations. Additionally, we evaluate the performance of SSBFGS and SSBroyden for Deep Operator Network (DeepONet) architectures, demonstrating their effectiveness for data-driven operator learning. Our findings provide state-of-the-art results with orders-of-magnitude accuracy improvements without the use of adaptive weights or any other enhancements typically employed in PINNs. More broadly, our work reveal insights into the effectiveness of quasi-Newton optimization strategies in significantly improving the convergence and accurate generalization of PINNs and PIKANs.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A deep learning model for automatic analysis of cavities in irradiated materials

Transmission electron microscopy (TEM) is a commonly used technique in materials science for defect investigation. Quantitative analysis of defects is important for understanding the properties of a material, but manual analysis of TEM micrographs can be time-consuming and prone to error, especially when the defects have irregular shapes rather than spherical shapes. Many existing methods or deep learning models do not handle a wide range of sizes for the same object type within a single image. In this work, we present a framework that enables users to train an instance segmentation model called Mask R- CNN on any microstructure dataset, perform multi-detection on the same image at different scales, and obtain properties (e.g., size, area) of the objects based on the desired shape (e.g., circle, ellipse, rectangle). Additionally, we have developed a parallel detection module that uses multiple GPUs to increase the efficiency of the object detection process. We demonstrate the capabilities of our framework using a set of TEM images of cavities with different shapes, size distributions, and background contrasts. Finally, we show that the performance of our model in terms of density, size, and swelling of the cavities is comparable to the human average and that our model achieves the highest recall value compared to existing methods due to the use of image multi-rescaling.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Adaptive elasticity policies for staging-based in situ visualization

In situ processing aims to alleviate the growing gap between computation and I/O capabilities by performing data processing close to the data source. In situ processing is widely used to process data generated by multiple data sources, including observation data from edge devices or scientific observational facilities and the simulation data generated by scientific computation on a high-performance computing (HPC) platform. For a scientific workflow that is run on an HPC platform and composed of a simulation program and an in situ data analytics or visualization (abbreviated as ana/vis) task, there is an implicit assumption that the computing resources assigned to the workflow keep static during the workflow execution. However, with the converging trend between the HPC and cloud computing platform, running the in situ ana/vis task in an elastic way is promising to decrease its overhead and improve its resource utilization rate. Resource elasticity represents the ability to change resource configurations such as the number of computing nodes/processes during workflow execution. An elastic job may dynamically adjust resource configurations; it may use a few resources at the beginning and more resources toward the end of the job when interesting data appear. However, it is hard to predict a priori how many computing nodes/processes need to be added/removed during the workflow execution to adapt to changing workflow needs. How to efficiently guide elasticity operations, such as growing or shrinking the number of processes used for in situ analysis during workflow execution, is an open-ended research question. In this article, we present adaptive elasticity policies that adopt workflow runtime information collected during workflow execution to predict how to trigger the addition/removal of processes in order to minimize in situ processing overhead. Taking in situ visualization tasks as an example, we integrate the presented elasticity policies into a staging-based elastic workflow and evaluate its efficiency in multiple elasticity scenarios. Compared with the situation without elasticity or with a static elasticity policy that uses a fixed number of processes for each rescaling operation, the adaptive elasticity policy can save overhead in finding a proper resource configuration and improve resource utilization efficiency. Furthermore, one experiment illustrates that the adaptive elasticity policy saves 41% of core-hours compared with the situation without the resource elasticity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Quantitative SANS and multi-model analysis of spacer-dependent micellization of urea-based gemini surfactants

The micellization behavior of urea-based cationic gemini surfactants was investigated using small-angle neutron scattering (SANS) with multi-model form factor analysis. A homologous series of surfactants with urea group included in the hydrophobic tail and polymethylene spacers consisting of two to ten methylene units was analyzed using three form factor models: a core–shell ellipsoid and two variants of homogeneous ellipsoids. The results from all models show a consistent trend of the micelle structures, confirming that the spacer length critically influences micellar geometry, aggregation number, and hydration. The surfactant with four CH 2 groups in the spacer formed the largest micelles with the highest aggregation number, while longer spacers led to progressively smaller, more compact aggregates. The shell hydration—quantified as the volume fraction of heavy water within the hydrophilic region—decreased systematically with increasing spacer length due to enhanced hydrophobicity of the headgroup-spacer region. Intermicellar interactions, modeled as screened Coulomb interaction using the rescaled mean spherical approximation (RMSA), revealed the strongest electrostatic repulsion for the case of four methylene groups in the spacer, corresponding to the highest micellar charge and largest interparticle spacing. The observed spacer-dependent trends were robust across all modeling approaches, demonstrating that the spacer length serves as a key structural determinant of self-assembly in this type of urea-based gemini systems. These findings provide insight into the design of gemini surfactants with tailored aggregation behavior for applications in drug delivery, nanostructure templating, and solubilization technologies.

Core–shell ellipsoid model↗
DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

The HydroBio Dataset: a new data resource for evaluating existing and potential hydropower capacity and freshwater biodiversity in the conterminous United States

Hydropower is a critical source of affordable and reliable electricity and energy system stability services in the United States. Opportunities to expand US hydropower production include retrofitting existing non-powered dams to produce power, retrofitting existing hydropower dams to improve efficiency or increase capacity, or constructing new hydropower infrastructure on currently unregulated river reaches. We created the HydroBio Dataset, which summarizes existing and potential hydropower capacity and freshwater biodiversity at the sub-basin scale in the conterminous US to contextualize existing and potential grid contributions with the freshwater ecosystems in which dams are situated. We demonstrate a use-case of this dataset by rescaling and comparing potential non-powered dam nominal capacity to rarity-threat-weighted freshwater species richness for sub-basins where both types of data exist. On average, normalized freshwater biodiversity exceeded normalized potential non-powered dam nominal capacity in these sub-basins. Potential non-powered dam nominal capacity was concentrated in sub-basins in the Upper Mississippi and Ohio hydrologic regions while freshwater biodiversity was concentrated in the South Atlantic-Gulf, Ohio, and Tennessee hydrologic regions. Additionally, non-powered dams and existing hydropower dams are located in sub-basins with similar indices of freshwater biodiversity. The HydroBio Dataset adds an additional ecological dimension of context to our understanding of current and potential future US hydropower capabilities and is a valuable decision support tool for stakeholders tasked with balancing gains in services to the US power grid with the public and environmental benefits of freshwater ecosystems.

Biodiversity↗
DOE OSTI2023

Towards elastic in situ analysis for high-performance computing simulations

In situ analysis and visualization have grown increasingly popular for enabling direct access to data from high-performance computing (HPC) simulations. As a simulation progresses and interesting physical phenomena emerge, however, the data produced may become increasingly complex, and users may need to dynamically change the type and scale of in situ analysis tasks being carried out and consequently adapt the amount of resources allocated to such tasks. To date, none of the production in situ analysis frameworks offer such an elasticity feature, and for good reason: the assumption that the number of processes could vary during run time would force developers to rethink software and algorithms at every level of the in situ analysis stack. In this paper we present Colza, a data staging service with elastic in situ visualization capabilities. We demonstrate the use of Colza with the Deep Water Impact and the AMR-Wind simulations, coupling them with the ParaView Catalyst and Ascent in situ libraries, and show that Colza enables dynamic rescaling of these widely-used frameworks with no interruption to the simulation or staging service. Here, we highlight the challenges of enabling such elasticity, which requires overcoming these frameworks' reliance on MPI, using distinct engineering approaches, namely dependency injection and dependency overload. To the best of our knowledge, this work is the first to enable elastic in situ visualization capabilities for HPC applications on top of existing production analysis tools.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Under-capacitated and over-powered? Rural austerity and asymmetrical negotiating relationships in US wind energy development

Though rural local governments are central actors in renewable energy development, local governments in the United States (US) remain systematically under-funded. This paper considers what the manifestations of austerity in local governments broadly and rural localities specifically mean for renewable energy development and for energy transitions. Drawing on a survey of 262 elected county officials with experience with wind energy in eight US states, this paper asks how local officials understand the impacts of wind development, how local governments are involved in wind energy negotiations, how the resources and expertise needed to navigate negotiations are distributed among counties, and analyze the relationship between local capacity, access to resources, and involvement in negotiations. We find that local officials express simultaneously affective and material concerns with the impacts of wind development and see negotiations with the developer as central to realizing local benefits. However, the expertise and staffing needed to negotiate with developers is less accessible to poorer or sparsely populated counties, and counties with lower overall revenues have narrower scopes of negotiation, and counties incre. Our results suggest that uneven rural capacity heightens an already asymmetrical relationship between localities and developers. In analyzing how infrastructure developments are shaped by relationships between localities and developers that are conditioned by austerity and (under)capacity, this paper contributes to and bridges scholarly discussions on rural austerity, rescaling, and renewable energy transitions. These results challenge conventional wisdoms around centralizing energy siting processes, contextualize popular and academic debates about opposition to renewable energy development, and highlight the need for rural reinvestment to realize meaningfully participatory energy developments.

Elmallah, Salma↗
DOE OSTI2025

Corrigendum to ‘Under-capacitated and over-powered? Rural austerity and asymmetrical negotiating relationships in US wind energy development’ [J. Rural Stud., 119 (2025) 1–14]

The authors regret that there is an incomplete sentence in the abstract of the article, and request that the portion “, and counties incre” be deleted from the abstract (found at the end of the sentence beginning with “However …”). The portion to be deleted is underlined and bolded below. The authors would like to apologise for any inconvenience caused. Current abstract: Though rural local governments are central actors in renewable energy development, local governments in the United States (US) remain systematically under-funded. This paper considers what the manifestations of austerity in local governments broadly and rural localities specifically mean for renewable energy development and for energy transitions. Drawing on a survey of 262 elected county officials with experience with wind energy in eight US states, this paper asks how local officials understand the impacts of wind development, how local governments are involved in wind energy negotiations, how the resources and expertise needed to navigate negotiations are distributed among counties, and analyze the relationship between local capacity, access to resources, and involvement in negotiations. We find that local officials express simultaneously affective and material concerns with the impacts of wind development and see negotiations with the developer as central to realizing local benefits. However, the expertise and staffing needed to negotiate with developers is less accessible to poorer or sparsely populated counties, and counties with lower overall revenues have narrower scopes of negotiation, and counties incre. Our results suggest that uneven rural capacity heightens an already asymmetrical relationship between localities and developers. In analyzing how infrastructure developments are shaped by relationships between localities and developers that are conditioned by austerity and (under)capacity, this paper contributes to and bridges scholarly discussions on rural austerity, rescaling, and renewable energy transitions. These results challenge conventional wisdoms around centralizing energy siting processes, contextualize popular and academic debates about opposition to renewable energy development, and highlight the need for rural reinvestment to realize meaningfully participatory energy developments.

Elmallah, Salma↗
DOE OSTI2022

How to correct relative voxel scale factors for calculations of vector-difference Fourier maps in cryo-EM

The atomic coordinates derived from cryo-electron microscopy (cryo-EM) maps can be inaccurate when the voxel scaling factors are not properly calibrated. Here, we describe a method for correcting relative voxel scaling factors between pairs of cryo-EM maps for the same or similar structures that are expanded or contracted relative to each other. We find that the correction of scaling factors reduces the amplitude differences of Fourier-inverted structure factors from voxel-rescaled maps by up to 20–30%, as shown by two cryo-EM maps of the SARS-CoV-2 spike protein measured at pH 4.0 and pH 8.0. This allows for the calculation of the difference map after properly scaling, revealing differences between the two structures for individual amino acid residues. Unexpectedly, the analysis uncovers two previously overlooked differences of amino acid residues in structures and their local structural changes. Furthermore, we demonstrate the method as applied to two cryo-EM maps of monomeric apo-photosystem II from the cyanobacteria Synechocystis sp. PCC 6803 and Thermosynechococcus elongatus. The resulting difference maps reveal many changes in the peripheral transmembrane PsbX subunit between the two species.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

A semiempirical transparency model for dual energy cargo radiography applications

Cargo containers passing through ports are scanned by non-intrusive inspection systems to search for concealed illicit materials. By using two photon beams with different energy spectra, dual energy inspection systems are sensitive to both the area density and the atomic number of cargo contents. Most literature on the subject assumes a simple exponential attenuation model for photon intensity in which only free streaming photons are detected. However, this approximation neglects second order effects such as scattering, leading to a biased model and thus incorrect material predictions. This work studies the accuracy of the free streaming model by comparing it to simulation outputs, finding that the model shows poor atomic number reconstruction accuracy at high- Z and suffers significantly if the source energy spectra and detector response function are not known exactly. To address these challenges, this work introduces a semiempirical transparency model which modifies the free streaming model by rescaling different components of the mass attenuation coefficient, allowing the model to capture secondary effects ignored by the free streaming model. The semiempirical model displays improvement agreement with simulated results at high- Z and shows excellent extrapolation to materials and thicknesses which were not included during the calibration step. Furthermore, this work demonstrates that the semiempirical model yields accurate atomic number predictions even when the source spectra and detector response are not known exactly. In conclusion, using the semiempirical model, manufacturers can perform a simple calibration to enable more precise Z reconstruction capabilities, which has the potential to significantly improve the performance of existing dual energy radiographic systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Thermal diffusion, exhaust gas recirculation and blending effects on lean premixed hydrogen flames

Thermodiffusively-unstable lean premixed hydrogen flames are investigated using two-dimensional direct numerical simulation employing finite-rate chemical kinetics. Three databases are generated focussing on the inclusion of the Soret effect, the recirculation of exhaust gas, and blending with methane. A simple rescaling of a classic thermal diffusion model is presented and shown to mimic multicomponent diffusion with very low computational cost and little-to-no loss in accuracy. It is also shown that a previously developed model for mean local flame speeds in lean premixed hydrogen flames can still be used provided Soret effects are taken into account in one-dimensional calculations. The addition of exhaust gas to the unburned mixture is found to enhance thermodiffusive instability; the primary mechanism for this was shown to be the highly-efficient third-body nature of water, with the reduction of adiabatic flame temperature a second-order effect. Again, the existing mean local flame speed model proved sufficient. Finally, blending with methane was found to reduce the thermodiffusive response of the flame, more so than the existing model suggests, despite adjustment of the fuel Lewis number; an adapted model is presented to account for this.

08 HYDROGEN↗
DOE OSTI2020

Dissolution of a cylindrical disk in Hele-Shaw flow: a conformal-mapping approach

Here we apply conformal mapping to find the evolving shapes of a dissolving cylinder in a potential flow. Similar equations can be used to describe melting in a flowing liquid phase. Results are compared with microfluidic experiments and numerical simulations. Shapes predicted by conformal mapping agree almost perfectly with experimental observations, after a modest (20 %) rescaling of the time. Finite-volume simulations show that the differences with experiment are connected to the underlying assumptions of the analytical model: potential flow and diffusion-limited dissolution. Approximate solutions of the equations describing the evolution of the shape of the undissolved solid can be derived from a Laurent expansion of the mapping function from the unit circle. Asymptotic expressions for the evolution of the area of the disk and the shift in its centre of mass have been derived at low and high Péclet number. Analytic approximations to the leading-order Laurent coefficients provide additional insight into the mechanisms underlying pore-scale dissolution.

58 GEOSCIENCES↗
DOE OSTI2022

Verifying raytracing/Fokker–Planck lower-hybrid current drive predictions with self-consistent full-wave/Fokker–Planck simulations

Raytracing/Fokker–Planck (FP) simulations used to model lower-hybrid current drive (LHCD) often fail to reproduce experimental results, particularly when LHCD is weakly damped. A proposed reason for this discrepancy is the lack of ‘full-wave’ effects, such as diffraction and interference, in raytracing simulations and the breakdown of the raytracing approximation. Previous studies of LHCD using non-Maxwellian full-wave/FP simulations have been performed, but these simulations were not self-consistent and enforced power conservation between the FP and full-wave code using a numerical rescaling factor. Here, we have created a fully self-consistent full-wave/FP model for LHCD that is automatically power conserving. This was accomplished by coupling an overhauled version of the non-Maxwellian TORLH full-wave solver and the CQL3D FP code using the Integrated Plasma Simulator. We performed converged full-wave/FP simulations of Alcator C-Mod discharges and compared them with raytracing. We found that excellent agreement in the power deposition profiles from raytracing and TORLH could be obtained, however, TORLH had somewhat lower current drive efficiency and broader power deposition profiles in some cases. This discrepancy appears to be a result of numerical limitations present in the TORLH model and a small amount of diffractional broadening of the TORLH wave spectrum. Our results suggest full-wave simulation of LHCD is likely not necessary as diffraction and interference represented only a small correction that could not account for the differences between simulations and experiment.

Physics↗
DOE OSTI2024

Simulating nonlinear optical processes on a superconducting quantum device

Simulating plasma physics on quantum computers is difficult because most problems of interest are nonlinear, but quantum computers are not naturally suitable for nonlinear operations. In weakly nonlinear regimes, plasma problems can be modelled as wave–wave interactions. In this paper, we develop a quantization approach to convert nonlinear wave–wave interaction problems to Hamiltonian simulation problems. We demonstrate our approach using two qubits on a superconducting device. Unlike a photonic device, a superconducting device does not naturally have the desired interactions in its native Hamiltonian. Nevertheless, Hamiltonian simulations can still be performed by decomposing required unitary operations into native gates. To improve experimental results, we employ a range of error-mitigation techniques. Apart from readout error mitigation, we use randomized compilation to transform undiagnosed coherent errors into well-behaved stochastic Pauli channels. Moreover, to compensate for stochastic noise, we rescale exponentially decaying probability amplitudes using rates measured from cycle benchmarking. We carefully consider how different choices of product-formula algorithms affect the overall error and show how a trade-off can be made to best utilize limited quantum resources. This study provides an example of how plasma problems may be solved on near-term quantum computing platforms.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
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↗