Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Structural error”

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 901 records · Page 50

Parameter splitting in dark energy: is dark energy the same in the background and in the cosmic structures?

We perform an empirical consistency test of General Relativity/dark energy by disentangling expansion history and growth of structure constraints. We replace each late-universe parameter that describes the behavior of dark energy with two meta-parameters: one describing geometrical information in cosmological probes, and the other controlling the growth of structure. If the underlying model (a standard wCDM cosmology with General Relativity) is correct, that is under the null hypothesis, the two meta-parameters coincide. If they do not, it could indicate a failure of the model or systematics in the data. We present a global analysis using state-of-the-art cosmological data sets which points in the direction that cosmic structures prefer a weaker growth than that inferred by background probes. This result could signify inconsistencies of the model, the necessity of extensions to it or the presence of systematic errors in the data. We examine all these possibilities. The fact that the result is mostly driven by a specific sub-set of galaxy clusters abundance data, points to the need of a better understanding of this probe.

79 ASTRONOMY AND ASTROPHYSICS↗

SU-F-T-285: Evaluation of a Patient DVH-Based IMRT QA System

Purpose: To evaluate the clinical performance of a patient DVH-based QA system for prostate VMAT QA. Methods: Mobius3D(M3D) is a QA software with an independent beam model and dose engine. The MobiusFX(MFX) add-on predicts patient dose using treatment machine log files. We commissioned the Mobius beam model in two steps. First, the stock beam model was customized using machine commissioning data, then verified against the TPS with 12 simple phantom plans and 7 clinical 3D plans. Secondly, the Dosimetric Leaf Gap(DLG) in the Mobius model was fine-tuned for VMAT treatment based on ion chamber measurements for 6 clinical VMAT plans. Upon successful commissioning, we retrospectively performed IMRT QA for 12 VMAT plans with the Mobius system as well as the ArcCHECK-3DVH system. Selected patient DVH values (PTV D95, D50; Bladder D2cc, Dmean; Rectum D2cc) were compared between TPS, M3D, MFX, and 3DVH. Results: During the first commissioning step, TPS and M3D calculated target Dmean for 3D plans agree within 0.7%±0.7%, with 3D gamma passing rates of 98%±2%. In the second commissioning step, the Mobius DLG was adjusted by 1.2mm from the stock value, reducing the average difference between MFX calculation and ion chamber measurement from 3.2% to 0.1%. In retrospective prostate VMAT QA, 5 of 60 MFX calculated DVH values have a deviation greater than 5% compared to TPS. One large deviation at high dose level was identified as a potential QA failure. This echoes the 3DVH QA result, which identified 2 instances of large DVH deviation on the same structure. For all DVH’s evaluated, M3D and MFX show high level of agreement (0.1%±0.2%), indicating that the observed deviation is likely from beam modelling differences rather than delivery errors. Conclusion: Mobius system provides a viable solution for DVH based VMAT QA, with the capability of separating TPS and delivery errors.

60 APPLIED LIFE SCIENCES↗

The Model Assessment Wizard (MAW): A Visualization System for Ontology Constraint Violations

Validation and verification of engineering models is important to understand potential weaknesses and issues in the model. This is accomplished through the application of constraint logic to the model. These models and the constraints put upon them can be represented through a graph structure. Here we give a visualization system to aid users understanding, locating, and fixing constraint violations in their systems. We give users several ways to narrow down on the specific errors and parts of the graph they’re interested in. Users have the opportunity to choose the types of errors that will be shown in the graph. Clustering is applied to the graph to help users narrow down their searches. Several other graph interactions are given to support discovery of constraint violations.

97 MATHEMATICS AND COMPUTING↗

SSTDR and FDR Detection of Un-Energized and Energized Cable Anomalies Including Thermal Degradation Using Machine Learning

Historically, cables are initially qualified for nuclear power plant use for 40 years. As plants extend their operating license to 60 and 80 years, continued use of these cables must shift to a performance-based approach since it is cost prohibitive to completely replace cables that are likely still capable of performing their design function. A variety of cable tests are available and are commonly applied during outages when the cables can be taken out of service. Frequency domain reflectometry (FDR) is one of these test methods that is being more broadly accepted and used because it not only detects anomalies along the cable with a low-voltage signal that does not stress the cable insulation, but the technique also locates the anomalies. This supports follow-up local inspection and local repair or partial replacement of a damaged cable segment. Currently, FDR testing is only applied to cables that are taken out of service since the test instrument would be damaged by operational voltages. A related technology that has found some acceptance in the aircraft and rail industry is spread spectrum time domain reflectometry (SSTDR). This technology has been implemented with a custom commercial instrument by LiveWire Innovation that is designed to operate on live cables up to 1000 volts and with a bandwidth of 48 MHz. Initial evaluation by the Pacific Northwest National Laboratory (PNNL) of the Live Wire system indicated that a broader bandwidth (BW) SSTDR may be better for many kinds of flaws. This led PNNL to develop an SSTDR laboratory instrument suitable for tests up to 500 MHz bandwidth. Testing on energized cables is also desirable for online monitoring systems so an inductive clamshell coupler was developed that allows energized cables to be tested up to at least 5 kV and likely higher voltage levels. Dielectric spectroscopy and tan delta testing plus various laboratory destructive tests were included in this data acquisition campaign directed to feed a machine learning (ML) study. With these kinds of developments, online energized cable tests may be possible with industrial adoption of such hardware advances but it will be completely impractical to have highly skilled data analysts continually examine these complex signals for indications of damage or compromised conditions. If online testing is to be implemented in new test hardware, it must be accompanied by software that can interpret the signals and alert plant operators of changing or degraded conditions. The thermally aged, shielded cable investigated here was separately treated for ML analysis. Visual analysis of electrical data showed generally increasing peaks where the cable entered and exited the oven. These peaks were not exactly aligned with expected locations, but these differences were attributed to velocity of propagation calibration errors. Only supervised ML was applied to the thermally aged data as this data was only available shortly before the committed publication date of this report. The supervised ML was structured to divide the 0 to 70-day responses as ‘normal’ from 0 to 35 days or ‘anomalous’ from 36 to 70 days, based on cable tensile elongation at break (EAB) insulation characterization. Using 80% of the data for training and 20% for testing, the supervised ML predicted normal versus anomalous was 70% accurate. Important conclusions include: • Accuracy to predict the presence of cable damage is improved from the 2023 effort by more training data. Weighted accuracies for comparisons among the instruments ranged from 67 to 89 % for unsupervised ML and 71 to 99% for supervised ML. • Based on the synthetic data tests, the unsupervised models are more generalizable to unseen anomalies. The Multi-Layer Perceptron classifier (MLP) model reported as high as 99.7% accuracy on the test data, but this dropped to 58.3% when tested on the synthetic data. In contrast, the unsupervised Pointwise model only achieved 89.7% accuracy on the experimental data but reported 78.3% accuracy on the synthetic data. • The best anomaly indicators are higher frequency (400 MHz BW) FDR data. Other tests may be interesting but for this study, this was the best predicter.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Multi-scale Material Parameter Identification Using LS-DYNA® and LS-OPT®

Ever-tightening regulations on fuel economy, and the likely future regulation of carbon emissions, demand persistent innovation in vehicle design to reduce vehicle mass. Classical methods for computational mass reduction include sizing, shape and topology optimization. One of the few remaining options for weight reduction can be found in materials engineering and material design optimization. Apart from considering different types of materials, by adding material diversity and composite materials, an appealing option in automotive design is to engineer steel alloys for the purpose of reducing plate thickness while retaining sufficient strength and ductility required for durability and safety. A project to develop computational material models for advanced high strength steel is currently being executed under the auspices of the United States Automotive Materials Partnership (USAMP) funded by the US Department of Energy. Under this program, new Third Generation Advanced High Strength Steel (i.e., 3GAHSS) are being designed, tested and integrated with the remaining design variables of a benchmark vehicle Finite Element model. The objectives of the project are to integrate atomistic, microstructural, forming and performance models to create an integrated computational materials engineering (ICME) toolkit for 3GAHSS. The mechanical properties of Advanced High Strength Steels (AHSS) are controlled by many factors, including phase composition and distribution in the overall microstructure, volume fraction, size and morphology of phase constituents as well as stability of the metastable retained austenite phase. The complex phase transformation and deformation mechanisms in these steels make the well-established traditional techniques obsolete, and a multi-scale microstructure-based modeling approach following the ICME [0]strategy was therefore chosen in this project. Multi-scale modeling as a major area of research and development is an outgrowth of the Comprehensive Test Ban Treaty of 1996 which banned surface testing of nuclear devices [1]. This had the effect that experimental work was reduced from large scale tests to multiscale experiments to provide material models with validation at different length scales. In the subsequent years industry realized that multi-scale modeling and simulation-based design were transferable to the design optimization of any structural system. Horstemeyer [1] lists a number of advantages of the use of multiscale modeling. Among these are: the reduction of product development time by alleviating costly trial-and-error iterations as well as the reduction of product costs through innovations in material, product and process designs. Multi-scale modeling can reduce the number of costly large scale experiments and can increase product quality by providing more accurate predictions. Research tends to be focussed on each particular length scale, which enhances accuracy in the long term. This paper serves as an introduction to the LS-OPT and LS-DYNA methodology for multi-scale modeling. It mainly focuses on an approach to integrate material identification using material models of different length scales. As an example, a multi-scale material identification strategy, consisting of a Crystal Plasticity (CP) material model and a homogenized State Variable (SV) model, is discussed and the parameter identification of the individual material models of different length scales is demonstrated. The paper concludes with thoughts on integrating the multi-scale methodology into the overall vehicle design.

Stander, Nielen↗

Efficient Precoding for Single Carrier Modulation in Multi-User Massive MIMO Networks

Abstract—Frequency domain (FD) multi-user detection (MUD) and precoding has been shown to be an effective means of approaching the theoretical per-user capacity for single carrier modulation (SCM) schemes in massive MIMO scenarios with highly dispersive channels. When a cyclic prefix is added to the SCM waveform, the circulant structures of the resulting convolutional channel matrix allows for relatively simple expressions for the FD detection and precoding. In this paper, we further refine the computationally efficient minimum mean squared error (MMSE) FD-MUD technique in a time-division duplexing (TDD) massive MIMO setup by pre-computing the scale factor needed to produce an unbiased estimate. The resulting scale factor and the matrix inverse that was calculated for the uplink is then reused for the downlink to form a computationally efficient precoding scheme with superior performance compared to zero-forcing (ZF) precoding. Performance analysis is provided and confirmed through simulation results.

5G and Beyond Communications↗

Development of a series of Korean pediatric and adult head computational phantoms and application to nuclear medicine dosimetry - Paper 26

Over the past decades, advances in Single Photon Emission Computed Tomography (SPECT) and Positron Emission Tomography (PET) imaging techniques and an increasing number of new radiopharmaceuticals have resulted in a surge of neuroimaging applications in nuclear medicine on pediatric and adult patients. It is not feasible to directly measure the absorbed dose to the organs of the human body from nuclear medicine procedures. The formalism of the Medical Internal Radiation Dose (MIRD) committee of the Society of Nuclear Medicine has been adopted as the calculation method to estimate radiation dose to target organs from radionuclides distributed in source organs. The absorbed dose to a target organ from different types of radionuclides distributed within source organs can be estimated by the products of the time-integrated activity in each source region and organ dose conversion coefficients, called S values (Gy/MBq-s). S values are computed from a comprehensive set of Specific Absorbed Fraction (SAF), which are calculated from computational human phantoms with realistic anatomical structures coupled with Monte Carlo transport technique. The current SAF data for pediatric and adult head and brain dosimetry is based on a series of stylized computational phantoms. However, dosimetric errors caused by the unrealistic anatomy of stylized phantoms in internal dosimetry calculations have been reported by several studies. In the current study, we developed a series of realistic computational head phantoms representing 6, 9, 12, 15, and 22-year-old standard Korean males, and computed a comprehensive set of SAFs using the phantoms coupled with Monte Carlo transport technique. (authors)

61 RADIATION PROTECTION AND DOSIMETRY↗

Binary pseudo-random array standards for calibration of 3D optical surface profilers used for metrology with significantly curved x-ray optics

High-accuracy metrology is vitally important in manufacturing ultra-high-quality free-form mirrors designed to manipulate X-ray light with nanometer-scale wavelengths. The current capabilities and possibility for improvements in X-ray mirror manufacturing are limited by inherent imperfections of the integrated metrology tools. In the case of Fizeau interferometry, metrology tools are currently calibrated with super-polished flat test-standard/reference mirrors. This is acceptable for measuring slightly curved X-ray optics. However, for even moderately curved aspherical X-ray mirrors the flat-reference calibration is not sufficiently accurate and stitching Fizeau interferometer-based surface metrology is used to mitigate the problem. But still, the retrace and aberration errors, as well as the limited spatial resolution, described with the instrument transfer function (ITF), can be transferred into the optical surface topography of X-ray mirrors obtained in stitching metrology. For ITF calibration, we have developed an original technique, based on test standards structured as two-dimensional (2D) highly-randomized (HR) binary pseudo-random arrays (BPRAs). The technique employs the unique properties of the HR BPRA patterns in the spatial frequency domain., i.e. the inherent 2D power spectral density of the HR BPRA pattern has a deterministic white-noise-like character that allows direct determination of the ITF with uniform sensitivity over the entire spatial frequency range and field-of-view of an instrument. Here, we explore technological, metrological, and analytical aspects essential for calibration of the retrace and aberration errors of Fizeau interferometers using different types of tilted test samples, including a super-polished reference mirror for the re-trace calibration and the uniformly redundant array (URA) BPRA standards for the geometrical distortion (aberration) calibration. While the first method was previously demonstrated by researchers at DIAMOND Light Source, a method based on the URA BPRA is described and demonstrated here for the first time. We outline the design and fabrication process used in fabrication of URA BPRA test standards, and present the results of application of the URA BPRA standards demonstrating the high efficacy of our approach to geometrical distortion calibration of Fizeau interferometers. We also discuss the possible sources of unexpected peculiarities of the systematic errors, including an astigmatic character of the retrace error, observed with Fizeau interferometers at the Advanced Light Source X-Ray Optics Laboratory.

Munechika, Keiko↗

Mapping of the HIDRA stellarator magnetic flux surfaces

The Hybrid Illinois Device for Research and Applications (HIDRA) is a classical stellarator designed for conducting plasma material interaction experiments and developing novel Plasma Facing Components (PFCs). Notably, the testing of two open-channel liquid lithium PFCs is imminent. Determining the shape of the plasma and its magnetic structure inside HIDRA is essential to carry out these tests. For this, electron traces were captured to build up the images of the HIDRA magnetic flux surfaces for several magnetic configurations, following the same procedure previously employed on the WEGA stellarator coupling an electron gun with a fluorescent detector. The FIELDLINES code has then been used to generate computational surfaces matching the experimental results. The obtained surfaces were found to be subject to a similar n = 1 error field as the one observed on WEGA, suggesting that the origin of this error field is inherent to the HIDRA vacuum vessel. Also, the effect of adding a vertical field was investigated, demonstrating the ability to radially shift the magnetic axis and move to a regime free of low-order rational resonances. Finally, this additional control over the HIDRA plasma and magnetic structure allows more freedom in setting up the PFC tests in the limiter and divertor regions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Data-driven material models for atomistic simulation

The central approximation made in classical molecular dynamics simulation of materials is the interatomic potential used to calculate the forces on the atoms. Great effort and ingenuity is required to construct viable functional forms and find accurate parametrizations for potentials using traditional approaches. Machine learning has emerged as an effective alternative approach to develop accurate and robust interatomic potentials. Beginning with a very general model form, the potential is learned directly from a database of electronic structure calculations and therefore can be viewed as a multiscale link between quantum and classical atomistic simulations. Risk of inaccurate extrapolation exists outside the narrow range of time and length scales where the two methods can be directly compared. Here, we use the spectral neighbor analysis potential (SNAP) and show how a fit can be produced with minimal interpolation errors which is also robust in extrapolating beyond training. To demonstrate the approach, we have developed a tungsten-beryllium potential suitable for the full range of binary compositions. Subsequently, large-scale molecular dynamics simulations were performed of high energy Be atom implantation onto the (001) surface of solid tungsten. The machine learned W-Be potential generates a population of implantation structures consistent with quantum calculations of defect formation energies. A very shallow ( < 2 nm ) average Be implantation depth is predicted which may explain ITER diverter degradation in the presence of beryllium.

36 MATERIALS SCIENCE↗

Cosmic voids and void lensing in the Dark Energy Survey science verification data

Galaxies and their dark matter halos populate a complicated filamentary network around large, nearly empty regions known as cosmic voids. Cosmic voids are usually identified in spectroscopic galaxy surveys, where 3D information about the large-scale structure of the Universe is available. Although an increasing amount of photometric data is being produced, its potential for void studies is limited since photometric redshifts induce line-of-sight position errors of ~50 Mpc/h or more that can render many voids undetectable. In this paper we present a new void finder designed for photometric surveys, validate it using simulations, and apply it to the high-quality photo-z redMaGiC galaxy sample of the Dark Energy Survey Science Verification (DES-SV) data. The algorithm works by projecting galaxies into 2D slices and finding voids in the smoothed 2D galaxy density field of the slice. Fixing the line-of-sight size of the slices to be at least twice the photo- z scatter, the number of voids found in these projected slices of simulated spectroscopic and photometric galaxy catalogs is within 20% for all transverse void sizes, and indistinguishable for the largest voids of radius ~70 Mpc/h and larger. The positions, radii, and projected galaxy profiles of photometric voids also accurately match the spectroscopic void sample. Applying the algorithm to the DES-SV data in the redshift range 0.2 < z < 0.8 , we identify 87 voids with comoving radii spanning the range 18-120 Mpc/h, and carry out a stacked weak lensing measurement. With a significance of 4.4σ, the lensing measurement confirms the voids are truly underdense in the matter field and hence not a product of Poisson noise, tracer density effects or systematics in the data. In conclusion, it also demonstrates, for the first time in real data, the viability of void lensing studies in photometric surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Blinded challenge for precision cosmology with large-scale structure: Results from effective field theory for the redshift-space galaxy power spectrum

An accurate theoretical template for the galaxy power spectrum is key for the success of ongoing and future spectroscopic surveys. We examine to what extent the effective field theory (EFT) of large-scale structure is able to provide such a template and correctly estimate cosmological parameters. To that end, we initiate a blinded challenge to infer cosmological parameters from the redshift-space power spectrum of high-resolution mock catalogs mimicking the BOSS galaxy sample but covering a 100 times larger cumulative volume. This gigantic simulation volume allows us to separate systematic bias due to theoretical modeling from the statistical error due to sample variance. The challenge is to measure three unknown input parameters used in the simulation: the Hubble constant, the matter density fraction, and the clustering amplitude. We present analyses done by two independent teams, who have fitted the mock simulation data generated by yet another independent group. This allows us to avoid any confirmation bias by analyzers and to pin down possible tuning of the specific EFT implementations. Both independent teams have recovered the true values of the input parameters within subpercent statistical errors corresponding to the total simulation volume.

79 ASTRONOMY AND ASTROPHYSICS↗

Tomographic constraints on gravity from angular redshift fluctuations in the late Universe

Fluctuations in sky maps of galaxy redshifts, dubbed angular redshift fluctuations (ARF), contain precise information about the growth rate of structures and the nature of gravity in the Universe. Specifically, ARFs constrain the combination of cosmological parameters H/H 0 fσ 8 (z), while being an intrinsically tomographic probe and largely insensitive to many observational systematic errors, all without requiring the assumption of any redshift-to-distance relation under a given fiducial cosmology. We present the first cosmological constraints derived from ARF by using Baryon Oscillation Spectroscopic Survey (BOSS) LOWZ(Low redshift(Z))+CMASS(Constant MASS) DR12(Data Release 12) galaxy samples, obtaining 7 per cent accurate constraints on H/H 0 fσ 8 (z) at more than 20 redshifts over the range z ϵ [0.26, 0.72]. Our best-fitting value is 10 per cent larger, but compatible at the 1.4σ level, than the ΛCDM expectation set by Planck observations of cosmic microwave background (CMB) radiation. In conclusion, our tomographic measurements, combined with these CMB data, provide one of the strongest constraints on the gravity index γ, $\gamma =0.44^{+0.09}_{-0.07}$, which lies within 2σ of the prediction of General Relativity (γGR ≃ 0.55).

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic equilibrium reconstructions of plasmas in the MAST database and preparation for reconstruction of the first plasmas in MAST upgrade

Reconstructions of plasma equilibria using magnetic sensors were routine during operation of the Mega Ampere Spherical Tokamak (MAST) device, but reconstructions using kinetic profiles were not. These are necessary for stability and disruption analysis of the MAST database, as well as for operation in the upgrade to the device, MAST-U. The three-dimensional (3D) code VALEN is used to determine eddy currents in the 3D vessel structures for vacuum coil test shots, which are then mapped to effective resistances in the two-dimensional vessel groupings in the EFIT equilibrium reconstruction code to be used in conjunction with nearby loop voltage measurements for estimated currents in the structures during reconstruction. Kinetic equilibrium reconstructions with EFIT, using all available magnetic sensors as well as Thomson scattering measurements of electron temperature and density, charge exchange recombination spectroscopy measurements of ion temperature, and internal magnetic field pitch angle measurements from a motional Stark effect (MSE) diagnostic are performed for a large database of MAST discharges. Excellent convergence errors are obtained for the portions of the discharges where the stored energy was not too low, and it is found that reconstructions performed with temperature and density measurements but without MSE data usually already match the pitch angle measurements well. A database of 275 kinetic equilibria is used to test the ideal MHD stability calculation capability for MAST. In conclusion, the necessary changes to conducting structure in VALEN, and diagnostic setup in EFIT have been completed for the upgrade from MAST to MAST-U, enabling kinetic reconstructions to commence from the first plasma discharges of the upgraded device.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Long-read, whole-genome shotgun sequence data for five model organisms

Single molecule, real-time (SMRT) sequencing from Pacific Biosciences is increasingly used in many areas of biological research including de novo genome assembly, structural-variant identification, haplotype phasing, mRNA isoform discovery, and base-modification analyses. High-quality, public datasets of SMRT sequences can spur development of analytic tools that can accommodate unique characteristics of SMRT data (long read lengths, lack of GC or amplification bias, and a random error profile leading to high consensus accuracy). In this paper, we describe eight high-coverage SMRT sequence datasets from five organisms (Escherichia coli, Saccharomyces cerevisiae, Neurospora crassa, Arabidopsis thaliana, and Drosophila melanogaster) that have been publicly released to the general scientific community (NCBI Sequence Read Archive ID SRP040522). Data were generated using two sequencing chemistries (P4C2 and P5C3) on the PacBio RS II instrument. The datasets reported here can be used without restriction by the research community to generate whole-genome assemblies, test new algorithms, investigate genome structure and evolution, and identify base modifications in some of the most widely-studied model systems in biological research.

59 BASIC BIOLOGICAL SCIENCES↗

Constraining the Multiscale Structure of Geophysical Fields in Machine Learning: The Case of Precipitation

The use of deep-learning algorithms for estimating the value of geophysical variables from remotely sensed information is rapidly expanding. The typical objective function minimized in such algorithms is the mean square error (MSE), which is known to lead to smooth estimates with compressed dynamical range as compared to the true distribution of the variable of interest. Here, we introduce and evaluate alternative objective functions, focusing on the retrieval of precipitation rates from satellite passive microwave radiometric measurements using a deep convolutional neural network. For this testbed application, the results show that explicitly imposing the preservation of the statistical distribution and spatial wavelet power spectrum of the target variable allows to accurately reproduce extreme values and sharp gradients across multiple scales in this study.

58 GEOSCIENCES↗

Molecular dynamics force-field refinement against quasi-elastic neutron scattering data

Quasi-elastic neutron scattering (QENS) is one of the experimental techniques of choice for probing the dynamics at length and time scales that are also in the realm of full-atom molecular dynamics (MD) simulations. This overlap enables extension of current fitting methods that use time-independent equilibrium measurements to new methods fitting against dynamics data. We present an algorithm that fits simulation-derived incoherent dynamical structure factors against QENS data probing the diffusive dynamics of the system. We showcase the difficulties inherent to this type of fitting problem, namely, the disparity between simulation and experiment environment, as well as limitations in the simulation due to incomplete sampling of phase space. We discuss a methodology to overcome these difficulties and apply it to a set of full-atom MD simulations for the purpose of refining the force-field parameter governing the activation energy of methyl rotation in the octa-methyl polyhedral oligomeric silsesquioxane molecule. Our optimal simulated activation energy agrees with the experimentally derived value up to a 5% difference, well within experimental error. We believe the method will find applicability to other types of diffusive motions and other representation of the systems such as coarse-grain models where empirical fitting is essential. In addition, the refinement method can be extended to the coherent dynamic structure factor with no additional effort.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗