A method for the efficient construction and sampling of vector-valued translation random fields.
Abstract not provided.
Grigoriu, Mircea↗
SEARCH · Search NASA
Results for “Randomized 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.
Towards developing a compact model for magnetization switching in straintronics magnetic random access memory devices
Strain-mediated magnetization switching in a magnetic tunneling junction (MTJ) by exploiting a combination of piezoelectricity and magnetostriction has been proposed as an energy efficient alternative to spin transfer torque (STT) and field induced magnetization switching methods in MTJ-based magnetic random access memories (MRAM). Theoretical studies have shown the inherent advantages of strain-assisted switching, and the dynamic response of the magnetization has been modeled using the Landau-Lifshitz-Gilbert (LLG) equation. However, an attempt to use LLG for simulating dynamics of individual elements in large-scale simulations of multi-megabyte straintronics MRAM leads to extremely time-consuming calculations. Hence, a compact analytical solution, predicting the flipping delay of the magnetization vector in the nanomagnet under stress, combined with a liberal approximation of the LLG dynamics in the straintronics MTJ, can lead to a simplified model of the device suited for fast large-scale simulations of multi-megabyte straintronics MRAMs. In this work, a tensor-based approach is developed to study the dynamic behavior of the stressed nanomagnet. First, using the developed method, the effect of stress on the switching behavior of the magnetization is investigated to realize the margins between the underdamped and overdamped regimes. The latter helps the designer realize the oscillatory behavior of the magnetization when settling along the minor axis, and the dependency of oscillations on the stress level and the damping factor. Next, a theoretical model to predict the flipping delay of the magnetization vector is developed and tested against LLG-based numerical simulations to confirm the accuracy of findings. Lastly, the obtained delay is incorporated into the approximate solutions of the LLG dynamics, in order to create a compact model to liberally and quickly simulate the magnetization dynamics of the MTJ under stress. Using the developed delay equation, the efficiency of the straintronics switching over the STT method is highlighted by analytically investigating the energy-delay trade-off of both methodologies.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Overview of Nuclear Integration Activities for Mitigation of Radiation Loads to ITER Superconducting Magnets
The Nuclear Integration Unit (NIU) is a trans-departmental team from the ITER Organization and Domestic Agencies, responsible for all nuclear integration aspects of the facility and for ensuring that systems and components are suitable for their radiation environment and satisfy all related performance requirements. This is a very challenging task namely due to the large number of systems, their still evolving designs, the intricate and far-reaching impacts of these in nuclear responses everywhere, and the complex world-wide distribution of procurement responsibilities. A sound systems engineering approach that involves all relevant experts and stakeholders, and trades-off between different performance requirements in an optimum way, needs to be taken and sustained during design evolution and construction, always with the objective of final overall project success. One of the most important current occupations of the NIU is the evaluation of, and protection from, radiation conditions to superconducting magnets during plasma operations, and particularly nuclear heating. This requires neutronic analyses in massive computer resources, with state-of-the-art methods and complex, extensive computer models. The main ITER superconducting magnet systems subject to significant nuclear loads are the toroidal field coils (TFCs) and the poloidal field coils (PFCs). These are illustrated in Fig. 1. Nuclear responses in these depend mostly on the design of all in-vessel (blankets, manifolds, divertor, in-vessel coils, port plugs), vessel and thermal shield components. Sources to consider are the 500 MW DT burning plasma (14.1 MeV neutrons) and the activated water in the in-vessel cooling system (photons up to 7 MeV from N-16 and neutrons up to 3.5 MeV from N-17). This paper discusses a snapshot of the present understanding of nuclear heat conditions in both the TFC and PFC magnets. This includes integral values and distributions, as well as best estimates from computations and conservative upper estimates by consideration of systematic effects and random uncertainties. The computer models, methods and assumptions used for evaluation of these responses are described. Nuclear integration aspects and ongoing efforts to study and implement additional shielding and other mitigating measures are also discussed.
70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
Molecular ecological network analyses
Background: Understanding the interaction among different species within a community and their responses to environmental changes is a central goal in ecology. However, defining the network structure in a microbial community is very challenging due to their extremely high diversity and as-yet uncultivated status. Although recent advance of metagenomic technologies, such as high throughout sequencing and functional gene arrays, provide revolutionary tools for analyzing microbial community structure, it is still difficult to examine network interactions in a microbial community based on high-throughput metagenomics data. Results: Here, we describe a novel mathematical and bioinformatics framework to construct ecological association networks named molecular ecological networks (MENs) through Random Matrix Theory (RMT)-based methods. Compared to other network construction methods, this approach is remarkable in that the network is automatically defined and robust to noise, thus providing excellent solutions to several common issues associated with high throughput metagenomics data. We applied it to determine the network structure of microbial communities subjected to long-term experimental warming based on pyrosequencing data of 16 S rRNA genes. We showed that the constructed MENs under both warming and unwarming conditions exhibited topological features of scale free, small world and modularity, which were consistent with previously described molecular ecological networks. Eigengene analysis indicated that the eigengenes represented the module profiles relatively well. In consistency with many other studies, several major environmental traits including temperature and soil pH were found to be important in determining network interactions in the microbial communities examined. To facilitate its application by the scientific community, all these methods and statistical tools have been integrated into a comprehensive Molecular Ecological Network Analysis Pipeline (MENAP), which is open-accessible now (http://ieg2.ou.edu/MENA). Conclusions: The RMT-based molecular ecological network analysis provides powerful tools to elucidate network interactions in microbial communities and their responses to environmental changes, which are fundamentally important for research in microbial ecology and environmental microbiology.
59 BASIC BIOLOGICAL SCIENCES↗
Computing rank‐revealing factorizations of matrices stored out‐of‐core
This paper describes efficient algorithms for computing rank-revealing factorizations of matrices that are too large to fit in main memory (RAM), and must instead be stored on slow external memory devices such as disks (out-of-core or out-of-memory). Traditional algorithms for computing rank-revealing factorizations (such as the column pivoted QR factorization and the singular value decomposition) are very communication intensive as they require many vector-vector and matrix-vector operations, which become prohibitively expensive when data is not in RAM. Randomization allows to reformulate new methods so that large contiguous blocks of the matrix are processed in bulk. The paper describes two distinct methods. The first is a blocked version of column pivoted Householder QR, organized as a “left-looking” method to minimize the number of the expensive write operations. The second method results employs a UTV factorization. It is organized as an algorithm-by-blocks to overlap computations and I/O operations. As it incorporates power iterations, it is much better at revealing the numerical rank. Numerical experiments on several computers demonstrate that the new algorithms are almost as fast when processing data stored on slow memory devices as traditional algorithms are for data stored in RAM.
97 MATHEMATICS AND COMPUTING↗
Wave Measurements Using Multi-Frame Processing of Marine Radar Data
Marine radars have proven to be useful for measuring ocean waves, but the accuracy of the measurements is limited by several factors including the look-angle dependence of the radar signals as well as noise in the radar data. The look-angle dependence introduces a systematic error or bias in the measurements, and noise causes a random error. This paper describes a method of combining data from multiple radar frames that is optimal in the sense of minimizing the error for a set of biased measurements with random additive noise. The results are shown experimentally to increase the correlation of the radar estimates with buoy measurements.
47 OTHER INSTRUMENTATION↗
Spillover as a cause of bias in baseline evaluation methods for demand response programs
Prior research has shown load reduction estimates from residential event-driven demand response programs (e.g., Critical Peak Pricing)using X of the highest Y days with a weather adjustment method are the best performing within the class of currently used baseline methods. However, they are still biased relative to estimates produced from randomized control trials (RCTs), the unbiased “gold standard” evaluation method. In this paper we identify underlying factors that cause some of the bias found in one commonly used baseline method. Rather than simply quantifying bias, our research provides a deeper understanding of what causes the bias and thus can be used to develop more accurate methods that are not subject to these underlying factors. Previous studies have compared various baseline methods relative to each other; however, because all baseline methods are biased, it is impossible to determine the true bias that exists in them. We have access to a unique RCT dataset: the Sacramento Municipal Utility District's study of critical peak pricing. Our analysis of 23 event days over two summers allows us to identify the true bias on load reduction estimates by using the RCT estimates as the unbiased gold standard against which we compare the estimates from the baseline methods. We found that spillover of energy reductions, from hours targeted by a program onto other hours, is one underlying factor that is a major cause of bias in baseline methods. We discuss alternative baseline methods that may not be subject to this same bias.
24 POWER TRANSMISSION AND DISTRIBUTION↗
Using the Monte Carlo Method to Evaluate the Reliability of Screening Multifamily Housing for Radon
When screening for radon in a multifamily housing complex using a fixed sample density (e.g., testing 1 in 10 (10%) or 1 in 4 units (25%)), the statistical confidence is dependent upon the assumed elevated radon frequency. For example, testing 25% of units in a complex estimated to have eight units with elevated radon levels will provide 90% confidence. However, if it is assumed that, in the same complex, there are only three units with elevated radon levels, the confidence drops to around 58% for the same sample density. Furthermore, in the previous example, if elevated radon levels are not found during the screening, all that can be stated is that the screening provides 90% confidence that there are no more than seven units in the complex with elevated radon levels. To more fully illustrate this uncertainty, ten separate multifamily housing radon data sets with 1 to 10 units with radon levels ≥4 pCi/L (Table 1) were selected for analysis using the Monte Carlo statistical method. Unlike other mathematically based statistical approaches, the Monte Carlo statistical method relies on repeated analysis of randomly selected data. from a 100% sampled complex at various sample densities. Success for each simulation is defined as finding at least one unit with elevated radon levels. In this statistical method, confidence in the overarching conclusion can be greatly enhanced by repeating the simulated screening hundreds or even thousands of times. For this study, each of the ten data sets was simulated 1,000 times. For each simulation, the data set was randomized three times before the fixed percentage of data was selected.
38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
Timing of Salvage Hormonal Therapy in Prostate Cancer Patients With Unfavorable Prognosis Treated With Radiotherapy: A Secondary Analysis of Radiation Therapy Oncology Group 85-31
Purpose: Radiation Therapy Oncology Group 85-31 was a randomized trial comparing radiotherapy (RT) alone vs. RT plus adjuvant androgen suppression for life in unfavorable-prognosis carcinoma of the prostate. We examined the impact of early initiation of salvage hormonal therapy (HT) in relapsing patients randomized to RT alone arm. Methods and Materials: Patients were divided into two groups: early salvage HT and late salvage HT. The early salvage group was defined as receiving HT with a prostate-specific antigen (PSA) level of less than 10 ng/mL, and the late salvage HT group had a PSA level of 10 ng/mL or greater. The outcomes were overall survival (OS), cause-specific mortality (CSM), and local failure (LF). The Kaplan-Meier estimation and log-rank test were used for OS, and the cumulative incidence estimation and Gray's test were used for CSM and LF. Proportional hazards regression models were used to compare the outcomes adjusted for other covariates. Results: The median follow-up times of surviving patients in the early and late salvage HT groups were about 11 and 13 years, respectively. The late salvage HT group had significantly more post-prostatectomy patients and patients with high Gleason scores. After adjustment for all covariates, OS was significantly longer in the early salvage HT group (hazard ratio, 1.5; p = 0.01). However, there were no statistically significant differences in LF or CSM between the groups. Conclusions: The early introduction of salvage HT resulted in improved OS but not improved CSM and LF. A randomized trial to define the optimal salvage hormonal timing is warranted in this group of patients with PSA recurrence after RT.
62 RADIOLOGY AND NUCLEAR MEDICINE↗
Embedded random phase approximation for magnetic systems: H 2 dissociative adsorption on Fe(110)
The random phase approximation (RPA), a method for treating electron correlation, has been shown to be superior to standard density functional theory (DFT) approximations in numerous cases. However, the RPA’s computational cost is substantially higher than that of DFT, particularly restricting its application to extended surfaces. The recently introduced embedded RPA (emb-RPA) approach [Wei et al., J. Chem. Phys. 159(19), 194108 (2023)] reduces this computational cost by approximately two orders of magnitude. While previous applications of emb-RPA focused on non-spin-polarized systems, here we extend the approach to ferromagnetic ones. Unlike other embedded correlated wavefunction methods, such as embedded complete active space self-consistent field theory, emb-RPA is advantageous for spin-polarized systems because the RPA is compatible with unrestricted DFT solutions, which are eigenfunctions of the spin angular momentum operator S z but not the total spin-squared operator S 2 . By applying emb-RPA with specific magnetization constraints, we achieved a speedup of two to three orders of magnitude (one order when accounting for the one-time embedding potential optimization cost) with only small errors (∼50 meV) compared to full periodic RPA. Moreover, emb-RPA significantly reduces the over-binding errors of DFT approximations. In conclusion, we anticipate that the acceleration enabled by the spin-polarized emb-RPA approach will broaden the applicability of RPA to magnetic materials.
Density functional theory↗
Dose and Fractionation in Radiation Therapy of Curative Intent for Non-Small Cell Lung Cancer: Meta-Analysis of Randomized Trials
Purpose: The optimum dose and fractionation in radiation therapy of curative intent for non-small cell lung cancer remains uncertain. We undertook a published data meta-analysis of randomized trials to examine whether radiation therapy regimens with higher time-corrected biologically equivalent doses resulted in longer survival, either when given alone or when given with chemotherapy. Methods and Materials: Eligible studies were randomized comparisons of 2 or more radiation therapy regimens, with other treatments identical. Median survival ratios were calculated for each comparison and pooled. Results: 3795 patients in 25 randomized comparisons of radiation therapy dose were studied. The median survival ratio, higher versus lower corrected dose, was 1.13 (95% confidence interval [CI] 1.04-1.22) when radiation therapy was given alone and 0.83 (95% CI 0.71-0.97) when it was given with concurrent chemotherapy (P for difference=.001). In comparisons of radiation therapy given alone, the survival benefit increased with increasing dose difference between randomized treatment arms (P for trend=.004). The benefit increased with increasing dose in the lower-dose arm (P for trend=.01) without reaching a level beyond which no further survival benefit was achieved. The survival benefit did not differ significantly between randomized comparisons where the higher-dose arm was hyperfractionated and those where it was not. There was heterogeneity in the median survival ratio by geographic region (P<.001), average age at randomization (P<.001), and year trial started (P for trend=.004), but not for proportion of patients with squamous cell carcinoma (P=.2). Conclusions: In trials with concurrent chemotherapy, higher radiation therapy doses resulted in poorer survival, possibly caused, at least in part, by high levels of toxicity. Where radiation therapy was given without chemotherapy, progressively higher radiation therapy doses resulted in progressively longer survival, and no upper dose level was found above which there was no further benefit. These findings support the consideration of further radiation therapy dose escalation trials, making use of modern treatment methods to reduce toxicity.
62 RADIOLOGY AND NUCLEAR MEDICINE↗
Scalable Extended Dynamic Mode Decomposition Using Random Kernel Approximation
The Koopman operator is a linear, infinite-dimensional operator that governs the dynamics of system observables; extended dynamic mode decomposition (EDMD) is a data-driven method for approximating the Koopman operator using functions (features) of the system state snapshots. This paper investigates an approach to EDMD in which the features used provide random approximations to a particular kernel function. The objective of this is computational economy for large data sets: EDMD is generally ill-suited for problems with large state dimension, and its dual kernel formulation (KDMD) is well-suited for such problems only if the number of data snapshots is relatively small. We discuss two specific methods for generating features: random Fourier features and the Nÿstrom method. The first method is a data-independent method for translation-invariant kernels only and involves random sampling in feature space; the second method is a data-dependent empirical method that may be used for any kernel and involves random sampling of data. We first discuss how these ideas may be applied in an EDMD context, as well as a means for adaptively adding random Fourier features. As a result, we demonstrate these methods on two example problems and conclude with an analysis of the relative benefits and drawbacks of each method.
59 BASIC BIOLOGICAL SCIENCES↗
Simple and statistically sound recommendations for analysing physical theories
Physical theories that depend on many parameters or are tested against data from many different experiments pose unique challenges to statistical inference. Many models in particle physics, astrophysics and cosmology fall into one or both of these categories. These issues are often sidestepped with statistically unsound ad hoc methods, involving intersection of parameter intervals estimated by multiple experiments, and random or grid sampling of model parameters. Whilst these methods are easy to apply, they exhibit pathologies even in low-dimensional parameter spaces, and quickly become problematic to use and interpret in higher dimensions. In this article we give clear guidance for going beyond these procedures, suggesting where possible simple methods for performing statistically sound inference, and recommendations of readily-available software tools and standards that can assist in doing so. Our aim is to provide any physicists lacking comprehensive statistical training with recommendations for reaching correct scientific conclusions, with only a modest increase in analysis burden. Furthermore, our examples can be reproduced with the code publicly available at Zenodo.
72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
Semi-device-independent random-number expansion without entanglement
By testing the classical correlation violation between two systems, true random numbers can be generated and certified without applying classical statistical method. In this work, we propose a true random-number expansion protocol without entanglement, where the randomness can be guaranteed only by the two-dimensional quantum witness violation. Furthermore, we only assume that the dimensionality of the system used in the protocol has a tight bound, and the whole protocol can be regarded as a semi-device-independent black-box scenario. Compared with the device-independent random-number expansion protocol based on entanglement, our protocol is much easier to implement and test.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Efficient Ab initio Modeling of Random Multicomponent Alloys
Here, we present in this Letter a novel small set of ordered structures (SSOS) method that allows extremely efficient ab initio modeling of random multi-component alloys. Using inverse II-III spinel oxides and equiatomic quinary bcc (so-called high entropy) alloys as examples, we also demonstrate that a SSOS can achieve the same accuracy as a large supercell or a well-converged cluster expansion, but with significantly reduced computational cost. In particular, because of this efficiency, a large number of quinary alloy compositions can be quickly screened, leading to the identification of several new possible high entropy alloy chemistries. Furthermore, the SSOS method developed here can be broadly useful for the rapid computational design of multi-component materials, especially those with a large number of alloying elements, a challenging problem for other approaches.
36 MATERIALS SCIENCE↗
Application of stochastic Galerkin FEM to the complete electrode model of electrical impedance tomography
The aim of electrical impedance tomography is to determine the internal conductivity distribution of some physical body from boundary measurements of current and voltage. The most accurate forward model for impedance tomography is the complete electrode model, which consists of the conductivity equation coupled with boundary conditions that take into account the electrode shapes and the contact resistances at the corresponding interfaces. If the reconstruction task of impedance tomography is recast as a Bayesian inference problem, it is essential to be able to solve the complete electrode model forward problem with the conductivity and the contact resistances treated as a random field and random variables, respectively. In this work, we apply a stochastic Galerkin finite element method to the ensuing elliptic stochastic boundary value problem and compare the results with Monte Carlo simulations.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Safety and effectiveness of a non-electric infant warmer for hypothermia in Rwanda: A cluster-randomized stepped-wedge trial
Background: Neonatal hypothermia is a common source of morbidity and mortality in low resource settings. We developed the Dream Warmer, a low cost, re-usable non-electric infant warmer to prevent and treat hypothermia. Methods: We conducted a cluster-randomized stepped-wedge trial. The primary aim was to assess the effect on overall euthermia rates of introducing the warmer compared to standard of care in rural Rwandan hospitals. The secondary aims were to assess effects of warmer introduction on mortality, as well as the safety and feasibility of the warmer. Ten district hospitals participated in the study from November 19th 2019 to July 15th 2020. Patients were eligible to use the warmer if they were 1) hypothermic (temp < 36·5 °C) or 2) or at risk of hypothermia (weight < 2·5 kg or estimated post menstrual age < 35 weeks) when Kangaroo Mother Care was not available. An encounter was defined as the data from an individual infant on a single day. Trial of a Non Electric Infant Warmer for Prevention and Treatment of Hypothermia in Rwanda [NCT03890211]. Findings: Over the study period, 3179 patients were enrolled across the ten neonatal wards, yielding 12,748 encounters; 464 unique infants used the warmer 892 times, 79% eligible due to hypothermia. Because of limited study nurse resources, the warmer was used in only 18% of eligible encounters. Despite this low rate of warmer use, the rate of euthermia rose from 51% (95% CI 50–52%) of encounters pre-intervention to 67% (66–68%) post-intervention; p < 0·0001. Among the encounters in which the warmer was used, only 11% (9–13%) remained hypothermic. While mortality rates pre- and post-intervention did not change, mortality rate among those who used the warmer was significantly lower than among those who did not (0·9% vs 2·8%, p = 0·01). Use of the warmer did not affect hyperthermia rates. There were no safety concerns or instances of incorrect warmer use. Interpretation: Introduction of the warmer increased rates of euthermia with no associated safety concerns.
60 APPLIED LIFE SCIENCES↗
TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations
TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy–cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe–Salpeter methods, second-order Møller–Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE’s functionality, including excited-state methods, RPA and Green’s function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE’s current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE’s development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗