Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Formal Methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,657 records · Page 92
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

SU-G-JeP2-08: Image-Guided Radiation Therapy Using Synthetic CTs in Brain Cancer

Purpose: Synthetic-CTs(synCTs) are essential for MR-only treatment planning. However, the performance of synCT for IGRT must be carefully assessed. This work evaluated the accuracy of synCT and synCT-generated DRRs and determined their performance for IGRT in brain cancer radiation therapy. Methods: MR-SIM and CT-SIM images were acquired of a novel anthropomorphic phantom and a cohort of 12 patients. SynCTs were generated by combining an ultra-short echo time (UTE) sequence with other MRI datasets using voxel-based weighted summation. For the phantom, DRRs from synCT and CT were compared via bounding box and landmark analysis. Planar (MV/KV) and volumetric (CBCT) IGRT performance were evaluated across several platforms. In patients, retrospective analysis was conducted to register CBCTs (n=34) to synCTs and CTs using automated rigid registration in the treatment planning system using whole brain and local registration techniques. A semi-automatic registration program was developed and validated to rigidly register planar MV/KV images (n=37) to synCT and CT DRRs. Registration reproducibility was assessed and margin differences were characterized using the van Herk formalism. Results: Bounding box and landmark analysis of phantom synCT DRRs were within 1mm of CT DRRs. Absolute 2D/2D registration shift differences ranged from 0.0–0.7mm for phantom DRRs on all treatment platforms and 0.0–0.4mm for volumetric registrations. For patient planar registrations, mean shift differences were 0.4±0.5mm (range: −0.6–1.6mm), 0.0±0.5mm, (range: −0.9–1.2mm), and 0.1±0.3mm (range: −0.7–0.6mm) for the superior-inferior(S-I), left-right(L–R), and anterior-posterior(A-P) axes, respectively. Mean shift differences in volumetric registrations were 0.6±0.4mm (range: −0.2–1.6mm), 0.2±0.4mm (range: −0.3–1.2mm), and 0.2±0.3mm (range: −0.2–1.2mm) for S-I, L–R, and A–P axes, respectively. CT-SIM and synCT derived margins were within 0.3mm. Conclusion: DRRs generated via synCT agreed well with CT-SIM. Planar and volumetric registrations to synCT-derived targets were comparable to CT. This validation is the next step toward clinical implementation of MR-only planning for the brain. The submitting institution has research agreements with Philips Healthcare. Research sponsored by a Henry Ford Health System Internal Mentored Grant.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2016

Consequences of Ca multisite occupation for the conducting properties of BaTiO{sub 3}

In combination with the dielectric modulus formalism and theoretical calculations, a newly developed defect incorporation mode, which is a combination of the standard A- and B-site doping mechanisms, is used to explain the conducting properties in 5 mol% Ca-doped BaTiO{sub 3}. Simulation results for Ca solution energies in the BaTiO{sub 3} lattice show that the new oxygen vacancy inducing mixed mode exhibits low defect energies. A reduction in dc conductivity compared with undoped BaTiO{sub 3} is witnessed for the incorporation of Ca. The conducting properties of 5 mol% Ca-doped BaTiO{sub 3} are analyzed using molecular dynamics and impedance spectroscopy. The ionic conductivity activation energies for each incorporation mode are calculated and good agreement with experimental data for oxygen migration is observed. The likely existence of the proposed defect configuration is also analyzed on the basis of these methods. - Graphical abstract: Oxygen vacancy formation as a result of self-compensation in Ca-doped BaTiO{sub 3}.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Constraining Polarized Foregrounds for EoR Experiments. II. Polarization Leakage Simulations in the Avoidance Scheme

A critical challenge in the observation of the redshifted 21 cm line is its separation from bright Galactic and extragalactic foregrounds. In particular, the instrumental leakage of polarized foregrounds, which undergo significant Faraday rotation as they propagate through the interstellar medium, may harmfully contaminate the 21 cm power spectrum. We develop a formalism to describe the leakage due to instrumental widefield effects in visibility-based power spectra measured with redundant arrays, extending the delay-spectrum approach presented in Parsons et al. We construct polarized sky models and propagate them through the instrument model to simulate realistic full-sky observations with the Precision Array to Probe the Epoch of Reionization. We find that the leakage due to a population of polarized point sources is expected to be higher than diffuse Galactic polarization at any k mode for a 30 m reference baseline. For the same reference baseline, a foreground-free window at k > 0.3 h Mpc{sup −1} can be defined in terms of leakage from diffuse Galactic polarization even under the most pessimistic assumptions. If measurements of polarized foreground power spectra or a model of polarized foregrounds are given, our method is able to predict the polarization leakage in actual 21 cm observations, potentially enabling its statistical subtraction from the measured 21 cm power spectrum.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Improving the User Interface of the DeepLynx Data Warehouse

DeepLynx is an open-source ontology-based data warehouse created by INL to support the creation and life cycle of digital engineering projects, with a particular emphasis on digital twins [1]. Digital twins are systems that represent physical assets and process in a real-time digital environment [1]. Most well-known commercial data warehouses use Graphical User Interfaces (GUIs) for users to interact with their systems [3]. Limited publications have addressed the design of these interfaces and understanding of their target users. The current users and development team acknowledge the need to improve the current UI, not just for aesthetics but to improve functionality and workflow of DeepLynx. Traditional data warehouse users are developers, data scientists and business analysts [2]. DeepLynx users have a vast range of experience using data warehouses, and diverse roles, including engineers, scientists and management positions. Because there is a broader audience of target users for DeepLynx than a typical data warehouse, it is essential that DeepLynx has a useable and intuitive user interface. To achieve this the team performed human-computer interaction methods, including a Heuristic Evaluation of current UI using Neilsen’s Usability Heuristic, create personas based on current users by designing a user survey, data analysis and develop of personas. Followed by a redesign of the UI following using Neilsen’s Usability Heuristic and Norman’s Principles of Interactive Design in industry standard software Figma. Lastly a Heuristic Evaluation of new UI design, using Neilsen’s Usability Heuristic and User testing of redesign UI and have a group of users complete a Thinking Aloud Test of the new UI. Preliminary results of the Heuristic Evaluation of current UI arise issue with Consistency and Standards, Visibility of System Status, Match System and Real World and Recognition Rather than Recall. These issues were addressed in the proposed redesign by applying Neilsen’s Usability Heuristic and Norman’s Principles of Interactive Design. Next steps include formalized list of lessons learned and design implications for future publications.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Statistical Methods for Effective Spatial and Temporal Scaling in Support of Probabilistic Performance Assessments - 17340

Performance assessments (PAs) for the disposal of low-level radioactive waste have traditionally been performed in a deterministic manner, with inputs specified as single values as if they carry no uncertainty. Modern PAs, however, explicitly incorporate uncertainty through use of probability distributions for inputs instead of single values. While this is clearly an improvement, challenges remain in specifying appropriate distributions, many of which relate to capturing the appropriate spatial and/or temporal scale given the PA model, and matching the scale(s) of data and information to that specified by the model. More careful attention to scaling will lead to more defensible characterizations of physical variability and uncertainty in model inputs, thus leading to more defensible information regarding outcomes of interest to be used in decision making, and hence more defensible decisions. This work introduces general strategies and considerations for explicitly incorporating considerations of scale into distribution development. The strategies are explained as a series of steps and then presented with an example that starts to illustrate the new statistical concepts. This work frames the problem in terms of hierarchical statistical modeling, providing important connections to current statistical approaches to formally and rigorously account for temporal and spatial scaling. The larger goal is to build a more formal and unified approach to distribution development to help teams of modelers, subject matter experts, and statisticians work together to develop distributions and justify the choices made in the process, ultimately leading to more informed decision-making. Without this rigor, statistical distributions are not specified correctly and can have unexpected consequences for modeling and sensitivity analysis. (authors)

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

Image Guided Radiation Therapy Using Synthetic Computed Tomography Images in Brain Cancer

Purpose: The development of synthetic computed tomography (CT) (synCT) derived from magnetic resonance (MR) images supports MR-only treatment planning. We evaluated the accuracy of synCT and synCT-generated digitally reconstructed radiographs (DRRs) relative to CT and determined their performance for image guided radiation therapy (IGRT). Methods and Materials: Magnetic resonance simulation (MR-SIM) and CT simulation (CT-SIM) images were acquired of an anthropomorphic skull phantom and 12 patient brain cancer cases. SynCTs were generated using fluid attenuation inversion recovery, ultrashort echo time, and Dixon data sets through a voxel-based weighted summation of 5 tissue classifications. The DRRs were generated from the phantom synCT, and geometric fidelity was assessed relative to CT-generated DRRs through bounding box and landmark analysis. An offline retrospective analysis was conducted to register cone beam CTs (n=34) to synCTs and CTs using automated rigid registration in the treatment planning system. Planar MV and KV images (n=37) were rigidly registered to synCT and CT DRRs using an in-house script. Planar and volumetric registration reproducibility was assessed and margin differences were characterized by the van Herk formalism. Results: Bounding box and landmark analysis of phantom synCT DRRs were within 1 mm of CT DRRs. Absolute planar registration shift differences ranged from 0.0 to 0.7 mm for phantom DRRs on all treatment platforms and from 0.0 to 0.4 mm for volumetric registrations. For patient planar registrations, the mean shift differences were 0.4 ± 0.5 mm (range, −0.6 to 1.6 mm), 0.0 ± 0.5 mm (range, −0.9 to 1.2 mm), and 0.1 ± 0.3 mm (range, −0.7 to 0.6 mm) for the superior-inferior (S-I), left-right (L-R), and anterior-posterior (A-P) axes, respectively. The mean shift differences in volumetric registrations were 0.6 ± 0.4 mm (range, −0.2 to 1.6 mm), 0.2 ± 0.4 mm (range, −0.3 to 1.2 mm), and 0.2 ± 0.3 mm (range, −0.2 to 1.2 mm) for the S-I, L-R, and A-P axes, respectively. The CT-SIM and synCT derived margins were <0.3 mm different. Conclusion: DRRs generated by synCT were in close agreement with CT-SIM. Planar and volumetric image registrations to synCT-derived targets were comparable with CT for phantom and patients. This validation is the next step toward MR-only planning for the brain.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2024

Passivated Contacts for Direct Wafer Product (Final Technical Report)

This TCF project developed a thin-oxide (SiO 2 )/polycrystalline silicon (poly-Si) passivated contact solar cell on CubicPV's (formally 1366 Technologies, Inc.) Direct Wafer® Product (DWP) kerfless wafers. The project used two NREL-developed technologies described in U.S. Patent No. 9,911,873, Hydrogenation of Passivated Contacts and U.S. Patent Application Serial No. 15/890,172, Doped Passivated Contacts . The project was motivated by a potential higher efficiency cell (compared to a PERC cell) using passivated contacts on the ultra-low cost kerfless wafers grown using the Direct Wafer process. The hope was to accelerate market adoption of the cell and wafer by delivering the lowest LCOE in the PV industry. The project tested both n-type and p-type SiO 2 /poly-Si passivated contacts grown by thermal oxidation and plasma enhanced chemical vapor deposition (PECVD) of the poly-Si layer on DWP with varying wafer resistivities. Both deposition techniques are industry standards and thus economically viable methods for commercializing the contacts. The results indicated that both n-type and p-type poly-Si passivated contacts can be formed on polycrystalline DWP wafers, but implied open-circuit voltages (i Voc ) were limited to below 0.65 mV (compared with ~ 730 mV on n-Cz wafers). Diffusion of H to the Si/SiO 2 /poly-Si interface was key to obtaining high i Voc values. In this study, H was diffused from a high-temperature SiN x layer deposited over the poly-Si layer during a high-temperature firing step, similar to one used for screen printed metals. The study concluded that poly-Si passivated contacts on DWP wafers passivated the surface of the wafers as well as PERC passivated surfaces, which use a less expensive dielectric layer stack. The project showed that Direct Wafer Product wafers grown by CubicPV could produce high i Voc values (~0.647 mV), which could produce a cell over 20% efficient with proper processing and metallization. These cells, though not economically viable in 2024 as a stand-alone cell, could be integrated with a wide-bandgap top solar cell to form a two-junction tandem cell that could be viable under certain circumstances. This is because the bottom cell of a 30%, two-terminal tandem only needs to be a 20% cell under one-sun conditions. Thus, the DWP could be an ideal low-cost wafer for tandems. The project also revealed that a TOPCon type cell could be formed on a p-type DWP wafer using a P-diffused emitter and a p-type poly-Si contact. In fact, the p-type version of the poly-Si contact out-performed the n-type version for a variety of wafer resistivities, from highly doped to lowly doped. This curiosity requires more work to understand because on Cz wafers, the n-type poly-Si contact is of much higher quality than the p-type version.

14 SOLAR ENERGY↗
DOE OSTI2014

Properties of gravitationally equilibrated Yukawa systems—A molecular dynamics study

Using 2D Molecular Dynamics simulation, the equilibrium and dynamical properties of a gravitationally equilibrated Yukawa liquid are investigated. We observe that due to asymmetry introduced in one direction by gravity, several interesting features arise. For example, for a given value of coupling parameter Γ, screening parameter κ, and according to a chosen value of gravitational force g (say in y-direction), the system is seen to exhibit super-, sub- or normal diffusion. Interestingly, x-averaged density profiles, unlike a barotropic fluid, acquires sharp, free surface with scale free linear y-dependence. As can be expected for a system with macroscopic gradients, self-diffusion calculated from Green-Kubo’s formalism does not agree with that obtained from Einstein-Smoluchowski diffusion. A 2D angular-radial pair correlation function g(r, θ) clearly indicates asymmetric features induced by gravity. We observe that due to compression in y-direction, though in liquid state for all values of gravity considered, the transverse mode is found to predominant as compared to the longitudinal mode, leading to a novel Anisotropic Solid-like Yukawa liquid.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

TU-FG-BRB-02: The Impact of Using Dual-Energy CT for Determining Proton Stopping Powers: Comparison Between Theory and Experiments

Purpose: To evaluate the clinical performance of dual-energy CT (DECT) in determining proton stopping power ratios (SPR) and demonstrate advantages over conventional single-energy CT (SECT). Methods: SECT and DECT scans of tissue-equivalent plastics as well as animal meat samples are performed with a Siemens SOMATOM Definition Flash. The methods of Schneider et al. (1996) and Bourque et al. (2014) are used to determine proton SPR on SECT and DECT images, respectively. Waterequivalent path length (WEPL) measurements of plastics and tissue samples are performed with a 195 MeV proton beam. WEPL values are determined experimentally using the depth-dose shift and dose extinction methods. Results: Comparison between CT-based and experimental WEPL is performed for 12 tissue-equivalent plastic as well as 6 meat boxes containing animal liver, kidney, heart, stomach, muscle and bones. For plastic materials, results show a systematic improvement in determining SPR with DECT, with a mean absolute error of 0.4% compared to 1.7% for SECT. For the meat samples, preliminary results show the ability for DECT to determine WEPL with a mean absolute value of 1.1% over all meat boxes. Conclusion: This work demonstrates the potential in using DECT for determining proton SPR with plastic materials in a clinical context. Further work is required to show the benefits of DECT for tissue samples. While experimental uncertainties could be a limiting factor to show the benefits of DECT over SECT for the meat samples, further work is required to adapt the DECT formalism in the context of clinical use, where noise and artifacts play an important role.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2015

THE BOLOCAM GALACTIC PLANE SURVEY. XII. DISTANCE CATALOG EXPANSION USING KINEMATIC ISOLATION OF DENSE MOLECULAR CLOUD STRUCTURES WITH {sup 13}CO(1-0)

We present an expanded distance catalog for 1710 molecular cloud structures identified in the Bolocam Galactic Plane Survey (BGPS) version 2, representing a nearly threefold increase over the previous BGPS distance catalog. We additionally present a new method for incorporating extant data sets into our Bayesian distance probability density function (DPDF) methodology. To augment the dense-gas tracers (e.g., HCO{sup +}(3-2), NH{sub 3}(1,1)) used to derive line-of-sight velocities for kinematic distances, we utilize the Galactic Ring Survey (GRS) {sup 13}CO(1-0) data to morphologically extract velocities for BGPS sources. The outline of a BGPS source is used to select a region of the GRS {sup 13}CO data, along with a reference region to subtract enveloping diffuse emission, to produce a line profile of {sup 13}CO matched to the BGPS source. For objects with a HCO{sup +}(3-2) velocity, ≈95% of the new {sup 13}CO(1-0) velocities agree with that of the dense gas. A new prior DPDF for kinematic distance ambiguity (KDA) resolution, based on a validated formalism for associating molecular cloud structures with known objects from the literature, is presented. We demonstrate this prior using catalogs of masers with trigonometric parallaxes and H II regions with robust KDA resolutions. The distance catalog presented here contains well-constrained distance estimates for 20% of BGPS V2 sources, with typical distance uncertainties ≲ 0.5 kpc. Approximately 75% of the well-constrained sources lie within 6 kpc of the Sun, concentrated in the Scutum-Centaurus arm. Galactocentric positions of objects additionally trace out portions of the Sagittarius, Perseus, and Outer arms in the first and second Galactic quadrants, and we also find evidence for significant regions of interarm dense gas.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

FOURIER ANALYSIS OF BLAZAR VARIABILITY: KLEIN–NISHINA EFFECTS AND THE JET SCATTERING ENVIRONMENT

The strong variability of blazars can be characterized by power spectral densities (PSDs) and Fourier frequency-dependent time lags. In previous work, we created a new theoretical formalism for describing the PSDs and time lags produced via a combination of stochastic particle injection and emission via the synchrotron, synchrotron self-Compton, and external Compton (EC) processes. This formalism used the Thomson cross section and simple δ-function approximations to model the synchrotron and Compton emissivities. Here we expand upon this work, using the full Compton cross section and detailed and accurate emissivities. Our results indicate good agreement between the PSDs computed using the δ-function approximations and those computed using the accurate expressions, provided the observed photons are produced primarily by electrons with energies exceeding the lower limit of the injected particle population. Breaks are found in the PSDs at frequencies corresponding to the cooling timescales of the electrons primarily responsible for the observed emission, and the associated time lags are related to the difference in electron cooling timescales between the two energy channels, as expected. If the electron cooling timescales can be determined from the observed time lags and/or the observed EC PSDs, then one could in principle use the method developed here to determine the energy of the external seed photon source for EC, which is an important unsolved problem in blazar physics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

A performance model for the communication in fast multipole methods on high-performance computing platforms

Exascale systems are predicted to have approximately 1 billion cores, assuming gigahertz cores. Limitations on affordable network topologies for distributed memory systems of such massive scale bring new challenges to the currently dominant parallel programing model. Currently, there are many efforts to evaluate the hardware and software bottlenecks of exascale designs. It is therefore of interest to model application performance and to understand what changes need to be made to ensure extrapolated scalability. The fast multipole method (FMM) was originally developed for accelerating N-body problems in astrophysics and molecular dynamics but has recently been extended to a wider range of problems. Its high arithmetic intensity combined with its linear complexity and asynchronous communication patterns make it a promising algorithm for exascale systems. In this paper, we discuss the challenges for FMM on current parallel computers and future exascale architectures, with a focus on internode communication. We focus on the communication part only; the efficiency of the computational kernels are beyond the scope of the present study. We develop a performance model that considers the communication patterns of the FMM and observe a good match between our model and the actual communication time on four high-performance computing (HPC) systems, when latency, bandwidth, network topology, and multicore penalties are all taken into account. To our knowledge, this is the first formal characterization of internode communication in FMM that validates the model against actual measurements of communication time. The ultimate communication model is predictive in an absolute sense; however, on complex systems, this objective is often out of reach or of a difficulty out of proportion to its benefit when there exists a simpler model that is inexpensive and sufficient to guide coding decisions leading to improved scaling. The current model provides such guidance.

Computer Science↗
DOE OSTI2013

Data Sharing Report for the Quantification of Removable Activity in Various Surveillance and Maintenance Facilities at the Oak Ridge National Laboratory Oak Ridge TN

The U.S. Department of Energy (DOE) Oak Ridge Office of Environmental Management (OR-EM) requested that Oak Ridge Associated Universities (ORAU), working under the Oak Ridge Institute for Science and Education (ORISE) contract, provide technical and independent waste management planning support using American Recovery and Reinvestment Act (ARRA) funds. Specifically, DOE OR-EM requested that ORAU plan and implement a sampling and analysis campaign targeting potential removable radiological contamination that may be transferrable to future personal protective equipment (PPE) and contamination control materials—collectively referred to as PPE throughout the remainder of this report—used in certain URS|CH2M Oak Ridge, LLC (UCOR) Surveillance and Maintenance (S&M) Project facilities at the Oak Ridge National Laboratory (ORNL). Routine surveys in Bldgs. 3001, 3005, 3010, 3028, 3029, 3038, 3042, 3517, 4507, and 7500 continuously generate PPE. The waste is comprised of Tyvek coveralls, gloves, booties, Herculite, and other materials used to prevent worker exposure or the spread of contamination during routine maintenance and monitoring activities. This report describes the effort to collect and quantify removable activity that may be used by the ORNL S&M Project team to develop radiation instrumentation “screening criteria.” Material potentially containing removable activity was collected on smears, including both masselin large-area wipes (LAWs) and standard paper smears, and analyzed for site-related constituents (SRCs) in an analytical laboratory. The screening criteria, if approved, may be used to expedite waste disposition of relatively clean PPE. The ultimate objectives of this effort were to: 1) determine whether screening criteria can be developed for these facilities, and 2) provide process knowledge information for future site planners. The screening criteria, if calculated, must be formally approved by Federal Facility Agreement parties prior to use for ORNL S&M Project PPE disposal at the Environmental Management Waste Management Facility (EMWMF). ORAU executed the approved sampling and analysis plan (SAP) (DOE 2013) while closely coordinating with ORNL S&M Project personnel and using guidelines outlined in the Waste Handling Plan for Surveillance and Maintenance Activities at the Oak Ridge National Laboratory, DOE/OR/01-2565&D2 (WHP) (DOE 2012). WHP guidelines were followed because the PPE waste targeted by this SAP is consistent with that addressed under the approved Waste Lot (WL) 108.1 profile for disposal at EMWMF—this PPE is a “future waste stream” as defined in the WHP. The SAP presents sampling strategy and methodology, sample selection guidelines, and analytical guidelines and requirements necessary for characterizing future ORNL S&M Project PPE waste. This report presents a review of the sample and analysis methods including data quality objectives (DQOs), required deviations from the original design, summary of field activities, radiation measurement data, analytical laboratory results, a brief presentation of results, and process knowledge summaries.

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

Calculation of room temperature conductivity and mobility in tin-based topological insulator nanoribbons

Monolayers of tin (stannanane) functionalized with halogens have been shown to be topological insulators. Using density functional theory (DFT), we study the electronic properties and room-temperature transport of nanoribbons of iodine-functionalized stannanane showing that the overlap integral between the wavefunctions associated to edge-states at opposite ends of the ribbons decreases with increasing width of the ribbons. Obtaining the phonon spectra and the deformation potentials also from DFT, we calculate the conductivity of the ribbons using the Kubo-Greenwood formalism and show that their mobility is limited by inter-edge phonon backscattering. We show that wide stannanane ribbons have a mobility exceeding 10{sup 6} cm{sup 2}/Vs. Contrary to ordinary semiconductors, two-dimensional topological insulators exhibit a high conductivity at low charge density, decreasing with increasing carrier density. Furthermore, the conductivity of iodine-functionalized stannanane ribbons can be modulated over a range of three orders of magnitude, thus rendering this material extremely interesting for classical computing applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Permutation-invariant distance between atomic configurations

We present a permutation-invariant distance between atomic configurations, defined through a functional representation of atomic positions. This distance enables us to directly compare different atomic environments with an arbitrary number of particles, without going through a space of reduced dimensionality (i.e., fingerprints) as an intermediate step. Moreover, this distance is naturally invariant through permutations of atoms, avoiding the time consuming associated minimization required by other common criteria (like the root mean square distance). Finally, the invariance through global rotations is accounted for by a minimization procedure in the space of rotations solved by Monte Carlo simulated annealing. A formal framework is also introduced, showing that the distance we propose verifies the property of a metric on the space of atomic configurations. Two examples of applications are proposed. The first one consists in evaluating faithfulness of some fingerprints (or descriptors), i.e., their capacity to represent the structural information of a configuration. The second application concerns structural analysis, where our distance proves to be efficient in discriminating different local structures and even classifying their degree of similarity.

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

Stochastic resonance during a polymer translocation process

We have studied the occurrence of stochastic resonance when a flexible polymer chain undergoes a single-file translocation through a nano-pore separating two spherical cavities, under a time-periodic external driving force. The translocation of the chain is controlled by a free energy barrier determined by chain length, pore length, pore-polymer interaction, and confinement inside the donor and receiver cavities. The external driving force is characterized by a frequency and amplitude. By combining the Fokker-Planck formalism for polymer translocation and a two-state model for stochastic resonance, we have derived analytical formulas for criteria for emergence of stochastic resonance during polymer translocation. We show that no stochastic resonance is possible if the free energy barrier for polymer translocation is purely entropic in nature. The polymer chain exhibits stochastic resonance only in the presence of an energy threshold in terms of polymer-pore interactions. Once stochastic resonance is feasible, the chain entropy controls the optimal synchronization conditions significantly.

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

Unveiling acoustic physics of the CMB using nonparametric estimation of the temperature angular power spectrum for Planck

Estimation of the angular power spectrum is one of the important steps in Cosmic Microwave Background (CMB) data analysis. Here, we present a nonparametric estimate of the temperature angular power spectrum for the Planck 2013 CMB data. The method implemented in this work is model-independent, and allows the data, rather than the model, to dictate the fit. Since one of the main targets of our analysis is to test the consistency of the ΛCDM model with Planck 2013 data, we use the nuisance parameters associated with the best-fit ΛCDM angular power spectrum to remove foreground contributions from the data at multipoles ℓ ≥50. We thus obtain a combined angular power spectrum data set together with the full covariance matrix, appropriately weighted over frequency channels. Our subsequent nonparametric analysis resolves six peaks (and five dips) up to ℓ ∼1850 in the temperature angular power spectrum. We present uncertainties in the peak/dip locations and heights at the 95% confidence level. We further show how these reflect the harmonicity of acoustic peaks, and can be used for acoustic scale estimation. Based on this nonparametric formalism, we found the best-fit ΛCDM model to be at 36% confidence distance from the center of the nonparametric confidence set—this is considerably larger than the confidence distance (9%) derived earlier from a similar analysis of the WMAP 7-year data. Another interesting result of our analysis is that at low multipoles, the Planck data do not suggest any upturn, contrary to the expectation based on the integrated Sachs-Wolfe contribution in the best-fit ΛCDM cosmology.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Nucleon resonance structure in the finite volume of lattice QCD

An approach for relating the nucleon resonances extracted from πN reaction data to lattice QCD calculations has been developed by using the finite-volume Hamiltonian method. Within models of πN reactions, bare states are introduced to parametrize the intrinsic excitations of the nucleon. We show that the resonance can be related to the probability P N* (E) of finding the bare state, N*, in the πN scattering states in infinite volume. We further demonstrate that the probability P V N* (E) of finding the same bare states in the eigenfunctions of the underlying Hamiltonian in finite volume approaches P N* (E) as the volume increases. Our findings suggest that the comparison of P N* (E) and P V N* (E) can be used to examine whether the nucleon resonances extracted from the πN reaction data within the dynamical models are consistent with lattice QCD calculation. We also discuss the measurement of P V N* (E) directly from lattice QCD. Furthermore, the practical differences between our approach and the approach using the Lüscher formalism to relate LQCD calculations to the nucleon resonance poles embedded in the data are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗