Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “theoretical models”

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 325 records · Page 18

Measuring differential particle correlations in relativistic nuclear collisions

This study explores the transverse momentum ( p T ) dependencies of symmetric and asymmetric correlations (SC and ASC) with one and two particles of interest in Au + Au collisions at 200 GeV. Leveraging a multiphase transport (AMPT) model, the investigation delves into the sensitivity of these correlations to the final state effects, providing valuable insights into their potential for constraining the final state effects' p T dependencies. The HIJING model is employed as a benchmark for nonflow correlations, shedding light on their impact on interpreting SC and ASC data. Moreover, the study points out that differential SC and ASC with one and two particles of interest (POIs) typically incorporate contributions from event-plane angle fluctuations. Consequently, this work highlights the significance of SC and ASC with one and two POIs as valuable tools for investigating the p T nature of the final state effects and advocates for comprehensive experimental measurements across various beam energies and system sizes to enhance understanding and provide additional constraints for theoretical models. Published by the American Physical Society 2024

Physics↗

Guidelines to Achieving High Selectivity for the Hydrogenation of α,β-Unsaturated Aldehydes with Bimetallic and Dilute Alloy Catalysts: A Review

Selective hydrogenation of unsaturated aldehydes to unsaturated alcohols is a challenging class of reactions, yielding valuable intermediates for the production of pharmaceuticals, perfumes, and flavorings. On monometallic heterogeneous catalysts, the formation of the unsaturated alcohols is thermodynamically disfavored over the saturated aldehydes. Hence, new catalysts are required to achieve the desired selectivity. In this work, the literature of three major research areas in catalysis is integrated as a step toward establishing the guidelines for enhancing the selectivity: reactor studies of complex catalyst materials at operating temperature and pressure, surface science studies of crystalline surfaces in ultrahigh vacuum, and first-principles modeling using density functional theory calculations. Aggregate analysis shows that bimetallic and dilute alloy catalysts significantly enhance the selectivity to the unsaturated alcohols compared to monometallic catalysts. This comprehensive review focuses primarily on the role of different metal surfaces as well as the factors that promote the adsorption of the unsaturated aldehyde via its C=O bond, most notably by electronic modification of the surface and formation of the electrophilic sites. Furthermore, challenges, gaps, and opportunities are identified to advance the rational design of efficient catalysts for this class of reactions, including the need for systematic studies of catalytic processes, theoretical modeling of complex materials, and model studies under ambient pressure and temperature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GLUE Code: A framework handling communication and interfaces between scales

Many scientific applications are inherently multiscale in nature. Such complex physical phenomena often require simultaneous execution and coordination of simulations spanning multiple time and length scales. This is possible by combining expensive small-scale simulations (such as molecular dynamics simulations) with larger scale simulations (such continuum limit/hydro solvers) to allow for considerably larger systems using task and data parallelism. However, the granularity of the tasks can be very large and often leads to load imbalance. Traditionally, we use approximations to streamline the computation of the more costly interactions and this introduces trade-offs between simulation cost and accuracy. In recent years, the available computational power and the advances in machine learning have made computing these scale-bridging interactions and multiscale simulations more feasible. One driving application has been plasma modeling in inertial confinement fusion (ICF), which is fundamentally multiscale in nature. This requires deep understanding of how to extrapolate microscopic information into macroscopically relevant scales. For example, in ICF one needs an accurate understanding of the connection between experimental observables and the underlying microphysics. The properties of the larger scales are often affected by the microscale behavior incorporated usually into the equations of state and ionic and electronic transport coefficients (Liboff, 1959; Rinderknecht et al., 2014; Rosenberg et al., 2015; Ross et al., 2017). Instead of incorporating this information using reliable molecular dynamics (MD) simulations, one often needs to use theoretical models, due to the inability of MD to reach engineering scales (Glosli et al., 2007; Marinak et al., 1998). One approach to resolve this issue is by coupling two MD simulations of different scales via force interpolation, e.g., the AdResS method (Krekeler et al., 2018; Nagarajan et al., 2013). Another approach, which we will pursue in the scope of this work, is by enabling scale bridging between MD simulations and meso/macro-scale models through the development and support of application programming interfaces that these different applications can interact through.

54 ENVIRONMENTAL SCIENCES↗

Visualizing electronic structure of twisted bilayer MoTe 2 in devices

The pursuit of emergent quantum phenomena lies at the forefront of modern condensed matter physics. A recent breakthrough in this arena is the discovery of the fractional quantum anomalous Hall effect (FQAHE) in twisted bilayer MoTe₂ (tbMoTe₂), marking a paradigm shift and establishing a versatile platform for exploring the intricate interplay among topology, magnetism, and electron correlations. While significant progress has been made through both optical and electrical transport measurements, direct experimental insights into the electronic structure – crucial for understanding and modeling this system – have remained elusive. Here, using spatially and angle-resolved photoemission spectroscopy (μ-ARPES), we directly map the electronic band structure of tbMoTe₂. We identify the valence band maximum, whose partial filling underlies the FQAHE, at the K points, situated approximately 150 meV above the Γ valley. By fine-tuning the doping level via in-situ alkali metal deposition, we also resolve the conduction band minimum at the K point, providing direct evidence that tbMoTe₂ exhibits a direct band gap – distinct from all previously known moiré bilayer transition metal dichalcogenide systems. These results offer critical insights for theoretical modeling and advance our understanding of fractionalized excitations and correlated topological phases in this emergent quantum material.

Chen, Cheng [Univ. of Oxford (United Kingdom)]↗

Electromagnetic and two-photon transition form factors of the pseudoscalar mesons: An algebraic model computation

We compute electromagnetic and two-photon transition form factors of ground-state pseudoscalar mesons: π , K , η c , η b . To this end, we employ an algebraic model based upon the coupled formalism of Schwinger-Dyson and Bethe-Salpeter equations. Within this approach, the dressed quark propagator and the relevant Bethe-Salpeter amplitude encode the internal structure of the corresponding meson. Electromagnetic properties of the meson are probed via the quark-photon interaction. The algebraic model employed by us unifies the treatment of all ground-state pseudoscalar mesons. Its parameters are carefully fitted performing a global analysis of existing experimental data including the knowledge of the charge radii of the mesons studied. We then compute and predict electromagnetic and two-photon transition form factors for a wide range of probing photon momentum-squared which is of direct relevance to the experimental observations carried out thus far or planned at different hadron physics facilities such as the Thomas Jefferson National Accelerator Facility (JLab) and the forthcoming Electron-Ion Collider. We also present comparisons with other theoretical models and approaches and lattice quantum chromodynamics. Published by the American Physical Society 2024

Higuera-Angulo, I. M. (ORCID:0000000256008875)↗

SDSS-IV MaNGA: Understanding Ionized Gas Turbulence Using Integral Field Spectroscopy of 4500 Star-forming Disk Galaxies

The Sloan Digital Sky Survey MaNGA program has now obtained integral field spectroscopy for over 10,000 galaxies in the nearby universe. We use the final MaNGA data release DR17 to study the correlation between ionized gas velocity dispersion and galactic star formation rate, finding a tight correlation in which σ Hα from galactic H ιι regions increases significantly from ~18–30 km s –1 , broadly in keeping with previous studies. In contrast, σ Hα from diffuse ionized gas increases more rapidly from 20–60 km s –1 . Using the statistical power of MaNGA, we investigate these correlations in greater detail using multiple emission lines and determine that the observed correlation of σ Hα with local star formation rate surface density is driven primarily by the global relation of increasing velocity dispersion at higher total star formation rate, as are apparent correlations with stellar mass. Assuming H ιι region models consistent with our finding that σ [O III] < σ Hα < σ [O I] , we estimate the velocity dispersion of the molecular gas in which the individual H ιι regions are embedded, finding values σ Mol = 5–30 km s –1 consistent with ALMA observations in a similar mass range. Finally, we use variations in the relation with inclination and disk azimuthal angle to constrain the velocity dispersion ellipsoid of the ionized gas σ z /σ r = 0.84 ± 0.03 and σ Φ /σ r = 0.91 ± 0.03, similar to that of young stars in the Galactic disk. Our results are most consistent with the theoretical models in which turbulence in modern galactic disks is driven primarily by star formation feedback.

79 ASTRONOMY AND ASTROPHYSICS↗

Toward a Machine Learning Approach to Interpreting X-ray Spectra of Trace Impurities by Converting XANES to EXAFS

The fact that the photoabsorption spectrum of a material contains information about the atomic structure, commonly understood in terms of multiple scattering theory, is the basis of the popular extended X-ray absorption spectroscopy (EXAFS) technique. How much of the same structural information is present in other complementary spectroscopic signals is not obvious. Here we use a machine learning approach to demonstrate that within theoretical models that accurately predict the EXAFS signal, the extended near-edge region does indeed contain the EXAFS-accessible structural information. We do this by exhibiting deep operator neural networks (DeepONets) that have learned the relationship between the extended and near edge portions of the X-ray absorption spectrum to predict the former from the latter. We find that we can accurately predict the EXAFS spectrum between 6 and 14 Å –1 from the first 6 Å –1 (≈100 eV) of the absorption spectrum of Cu 2 + substitutional defects in the Fe 3+ mineral hematite (α-Fe 2 O 3 ). This surprising finding implies that theoretical analyses of X-ray absorption spectra could be implemented that extract the same conclusions as high-quality EXAFS studies from spectra collected over a much smaller range of photon energies. This relaxes a host of experimental limitations related to the X-ray source and measurement sample, including collection time, minimum dopant concentration, source brilliance, and energy range. We describe the theoretical data sets and DeepONet construction and show that the resulting DeepONets produce EXAFS that recovers linear combination fits to experimental data with accuracy approaching the original ab initio calculations. We discuss the implications of our findings for minor constituent characterization and for understanding the information content of spectroscopic data more broadly, including how this approach might be applied to measured experimental spectra. In conclusion, to encourage similar efforts, the simulated X-ray spectra, machine learning, and fitting code are publicly available.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Exploring a New Regime of Molecular Orbital Physics in 4⁢d Cluster Magnets with Resonant Inelastic X-Ray Scattering

Molecular orbital systems with clusters of heavy transition metal (TM) ions are one of the most important classes of model materials for studying the interplay between local physics and effects of itinerancy. Despite a large number of candidates identified in the family of 4⁢d TM materials, an understanding of their physics from competing microscopic energy scales is still missing. Here, we bridge this gap by reporting the first resonant inelastic x-ray scattering (RIXS) measurement on a well-known series of Ru cluster magnets with a 6⁢H-perovskite structure Ba 3 ⁢M⁢Ru 2 ⁢O 9 (M 3+ = In 3+ , Y 3+ , La 3+ ) comprised of Ru dimers. In addition to providing a microscopic explanation for their anomalous magnetic properties, our RIXS measurements combined with theoretical modeling uncover a new regime of molecular orbital physics where the combined effect of large hopping and small spin-orbit coupling results in highly fragile electronic states in the Ru-dimer compounds directly manifested as an abrupt change in the RIXS spectrum accompanying a tiny change in the local structure tuned by the M-site ion. This unique combination of energy scales found only in the 4⁢d but not the 5⁢d cluster magnets highlights the value of these materials as a new platform for studying quantum phase transition involving molecular orbitals.

36 MATERIALS SCIENCE↗

Weak lensing in DES Y3: redshift distributions, shape catalogue, and mass mapping

The size of galaxy surveys have rapidly escalated with time, growing from thousands of galaxies of the first pioneering surveys to current and future surveys (DES, KiDS, HSC, LSST, Euclid, WFIRST), which map (or plan to map) a few hundreds millions of galaxies. These ambitious programs have the potential to answer many open questions in Cosmology (or, depending on their results, to open new exciting scenarios, involving new physics beyond the standard cosmological model and the standard model of particle physics). In particular, this thesis focuses on the first three years of data from the Dark Energy Survey (DES Y3), which comprises 100 million galaxies observed in 5 different optical bands over 5000 square degrees of the southern hemisphere. DES and the data set used in this thesis are discussed in detail in Chapter 3, which together with Chapter 1 and 2 constitutes Part I of this thesis. The scientifically positive outcome of DES and other photometric galaxy surveys goes h and in h and with a supreme understanding of the data, of the observational systematics affecting the measurements, and of the limits of the theoretical modelling of the observables. One of the main challenge of galaxy photometric surveys like DES is to correctly estimate the distance to the galaxies. Photometric surveys observe the galaxies in a limited number of broad optical bands, and the estimate of their \textit{redshifts} is somewhat uncertain, which constitutes one of the main limiting factor of their cosmological analyses. Part II and its Chapter 4 are devoted to a new auxiliary technique recently developed that allows to infer (or calibrate) redshift distributions exploiting cross-correlation with small reference samples with secure redshifts. The so called ``clustering-redshift'' method is discussed in depth, as well as its role in the overall DES Y3 redshift calibration strategy. Part III and Chapter 5 are devoted to another key aspect of photometric surveys, that enables the stu dy of weak gravitational lensing: the shape catalogue. The official DES Y3 shape catalogue is presented, together with all the systematic tests performed to validate it such that it can be considered science-ready. Last, part IV focuses on how to use weak gravitational lensing and the DES Y3 shape catalogue to create weak lensing mass maps, i.e., maps of the mass distribution of the Universe. Chapter 6 discusses different map making techniques, whereas Chapter 7 presents a simulated cosmological analysis using the high order moments of the mass maps.

79 ASTRONOMY AND ASTROPHYSICS↗

Materials synthesis at terapascal static pressures

Theoretical modelling predicts very unusual structures and properties of materials at extreme pressure and temperature conditions. Hitherto, their synthesis and investigation above 200 gigapascals have been hindered both by the technical complexity of ultrahigh-pressure experiments and by the absence of relevant in situ methods of materials analysis. Here we report on a methodology developed to enable experiments at static compression in the terapascal regime with laser heating. We apply this method to realize pressures of about 600 and 900 gigapascals in a laser-heated double-stage diamond anvil cell, producing a rhenium–nitrogen alloy and achieving the synthesis of rhenium nitride Re 7 N 3 —which, as our theoretical analysis shows, is only stable under extreme compression. Full chemical and structural characterization of the materials, realized using synchrotron single-crystal X-ray diffraction on microcrystals in situ, demonstrates the capabilities of the methodology to extend high-pressure crystallography to the terapascal regime.

36 MATERIALS SCIENCE↗

β -delayed neutron emission of Mn 64 , Cr 62 , and Fe 65

The β-decay properties of nuclei near the second nuclear “island of inversion” around neutron rich nuclei with neutron number 40 are important tests of nuclear structure models and interactions. In particular, the β-delayed neutron emission branch (Ρ n ), is useful for investigating β-strength and neutron-γ competition above the neutron separation energies of the daughter nuclei. We report new constraints for Ρ n values for three nuclei in the region: 62 Cr (Ρ n <1%), 64 Mn (Ρ n = 1.5⁢(6)%), and 65 Fe (Ρ n < 1%), measured with the Neutron Emission Ratio Observer (NERO) neutron long counter system and the Beta Counting Station (BCS) at the National Superconducting Cyclotron Laboratory (NSCL). Our results resolve the large discrepancy between previous direct and indirect measurements for 64 Mn and confirm the predictions of global theoretical models when a statistical treatment of the γ and neutron decays of the daughter states is included. Here we also obtain improved half-lives for 62 Cr [206(5) ms] and the short-lived isomer in the 62 Fe daughter [112(7) ms] from β-delayed γ emission data obtained in the same experiment with the Summing NaI (SuN) total absorption spectrometer. Finally, we use γ emission data to obtain a new upper limit for the 62 Cr β-decay population of the long-lived isomeric state in 62 Mn.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

𝑁 = 8 Shell Breaking in 12 Be from a Single-Particle Perspective

Experimental observations of the low-lying states in 12 Be and their accurate modeling play an essential role in understanding the disappearance of the 𝑁 = 8 magic number. Long-standing experimental ambiguities have been clarified using an one-neutron adding (𝑑, 𝑝) reaction on 11 Be using the ISOLDE Solenoidal Spectrometer at CERN’s HIE-ISOLDE facility. The single-particle energies of 1⁢𝑠 1/2 , 0⁢𝑑 5/2 , and 0⁢𝑝 1/2 orbitals in 12 Be have been determined from the extracted spectroscopic factors. A significant reduction between the separation of 1⁢𝑠 1/2 and 0⁢𝑝 1/2 orbitals is found in comparison with the carbon isotones, highlighting the breakdown of the 𝑁 = 8 shell. These observations serve as an important test of different effects incorporated in theoretical models. It is found that two synergistic mechanisms, core deformation and weak binding, are responsible for the 𝑁 = 8 shell breaking and the exotic near-threshold phenomena observed in 12 Be, including the narrow unnatural-parity resonance $0^{-}_1$ and the possible halolike nature of the $0^+_2$ isomer.

Chen, Jie [Southern University of Science and Tech↗

Beam Energy Dependence of the Linear and Mode-Coupled Flow Harmonics Using the a Multi-Phase Transport Model

In the framework of the A Multi-Phase Transport (AMPT) model, the multi-particle azimuthal cumulant method is used to calculate the linear and mode-coupled contributions to the quadrangular flow harmonic (v4) and the mode-coupled response coefficient as functions of centrality in Au+Au collisions at sNN = 200, 39, 27 and 19.6 GeV. This study indicates that the linear and mode-coupled contributions to v4 are sensitive to beam energy change. Nevertheless, the correlations between different-order flow symmetry planes and the mode-coupled response coefficients show weak beam energy dependence. In addition, the presented results suggest that the experimental measurements that span a broad range of beam energies can be an additional constraint for the theoretical model calculations.

79 ASTRONOMY AND ASTROPHYSICS↗

Characterizing the impact of water uptake on light absorption by aerosol particles

Absorbing aerosols (AA) have an important impact on the global radiation budget. The composition and properties of AA can vary substantially throughout the atmosphere, depending on the particle source and the influence of chemical aging. Uncertainties associated with the radiative effects of AA remain substantial. A key contributor to this uncertainty is understanding the extent to which water uptake alters absorption by AA particles and how this depends on particle composition. We have used new and existing experimental tools to systematically characterize the relationship(s) between AA chemical and physical characteristics, relative humidity, and mixing of AA with other aerosol components through a combination of laboratory and field measurements. We have developed fundamental understanding through laboratory experiments that will enable interpretation of field observations at DOE ARM sites and facilitate process-based improvements in the simulation of absorption by carbonaceous aerosols in climate models. In the laboratory experiments we have considered a variety of different AA types, including fullerene soot (a surrogate for refractory black carbon, or BC) and nigrosine (a surrogate for partially soluble brown carbon, or BrC). By examining a range of AA types we will develop insight into how the chemical and physical properties of AA particles (e.g. whether they are refractory or soluble) impact the influence of water uptake on absorption. These AA particles will also be mixed—both internally and externally—with compounds having a wide range of hygroscopic properties to establish how the impact of water uptake on absorption will vary with photochemical processing. These laboratory observations have been compared with commonly used theoretical models, enabling development of empirical models that will facilitate improved representation of absorption by AA particles in regional and global climate models. Overall, our study has contributed to the mission of the Atmospheric System Research program by quantifying how interactions between aerosols and radiation depend on relative humidity. Through this, our work will improve understanding and model representation of aerosol processes as they affect the Earth’s radiation budget.

54 ENVIRONMENTAL SCIENCES↗

SDSS-IV MaStar: Stellar parameter determination with continuum-supplemented full-spectrum fitting

We present a stellar parameter catalog built to accompany the MaStar Stellar Library, which is a comprehensive collection of empirical, medium-resolution stellar spectra. We constructed this parameter catalog by using a multicomponent χ 2 fitting approach to match MaStar spectra to models generated by interpolating the ATLAS9-based BOSZ model spectra. The total χ 2 for a given model is defined as the sum of components constructed to characterize narrow-band features of observed spectra (e.g., absorption lines) and the broadband continuum shape separately. Extinction and systematics due to flux calibration were taken into account in the fitting. The χ 2 distribution for a given region of model space was sampled using a Markov chain Monte Carlo (MCMC) algorithm, the data from which were then used to extract atmospheric parameter estimates (T eff , log g, [Fe/H], and [α/Fe]), their corresponding uncertainties, and direct extinction measurements. Two methods were used to extract parameters and uncertainties: one that accepts the MCMC’s prescribed minimum-χ 2 result, and one that uses Bayesian inference to compute a likelihood-weighted mean from the χ 2 distribution sampled by the MCMC. Results were evaluated for internal consistency using repeat observations where available and by comparing them with external data sets (e.g., APOGEE-2 and Gaia DR2). Our spectral-fitting exercise reveals possible deficiencies in current theoretical model spectra, illustrating the potential power of MaStar spectra for helping to improve the models. This paper represents an update to the parameters that were originally presented with SDSS-IV DR17. The MaStar parameter catalog containing our BestFit results is available on the SDSS-IV DR17 website as part of version 2 of the MaStar stellar parameter value-added catalog.

79 ASTRONOMY AND ASTROPHYSICS↗

Bayesian Conavigation: Dynamic Designing of the Material Digital Twins via Active Learning

Scientific advancement is universally based on the dynamic interplay between theoretical insights, modeling, and experimental discoveries. However, this feedback loop is often slow, including delayed community interactions and the gradual integration of experimental data into theoretical frameworks. This challenge is particularly exacerbated in domains dealing with high-dimensional object spaces, such as molecules and complex microstructures. Hence, the integration of theory within automated and autonomous experimental setups, or theory in the loop-automated experiment, is emerging as a crucial objective for accelerating scientific research. The critical aspect is to use not only theory but also on-the-fly theory updates during the experiment. Furthermore, we introduce a method for integrating theory into the loop through Bayesian conavigation of theoretical model space and experimentation. Our approach leverages the concurrent development of surrogate models for both simulation and experimental domains at the rates determined by latencies and costs of experiments and computation, alongside the adjustment of control parameters within theoretical models to minimize epistemic uncertainty over the experimental object spaces. This methodology facilitates the creation of digital twins of material structures, encompassing both the surrogate model of behavior that includes the correlative part and the theoretical model itself. While being demonstrated here within the context of functional responses in ferroelectric materials, our approach holds promise for broader applications, such as the exploration of optical properties in nanoclusters, microstructure-dependent properties in complex materials, and properties of molecular systems.

Microscopy↗

Three-Wave Mixing Kinetic Inductance Traveling-Wave Amplifier with Near-Quantum-Limited Noise Performance

We present a theoretical model and experimental characterization of a microwave kinetic inductance traveling-wave (KIT) amplifier, whose noise performance, measured by a shot-noise tunnel junction (SNTJ), approaches the quantum limit. Biased with a dc current, this amplifier operates in a three-wave mixing fashion, thereby reducing by several orders of magnitude the power of the microwave pump tone and associated parasitic heating compared to conventional four-wave mixing KIT amplifier devices. It consists of a 50 Ω artificial transmission line whose dispersion allows for a controlled amplification bandwidth. We measure $16.5^{+1}_{–1.3}$ dB of gain across a 2 GHz bandwidth with an input 1 dB compression power of –63 dBm, in qualitative agreement with theory. Using a theoretical framework that accounts for the SNTJ-generated noise entering both the signal and idler ports of the KIT amplifier, we measure the system-added noise of an amplification chain that integrates the KIT amplifier as the first amplifier. This system-added noise, 3.1 ± 0.6 quanta (equivalent to 0.66 ± 0.15 K) between 3.5 and 5.5 GHz, is the one that a device replacing the SNTJ in that chain would see. This KIT amplifier is therefore suitable to read large arrays of microwave kinetic inductance detectors and promising for multiplexed superconducting qubit readout.

42 ENGINEERING↗

Inference for Nonparanormal Partial Correlation via Regularized Rank-Based Nodewise Regression

Abstract Partial correlation is a common tool in studying conditional dependence for Gaussian distributed data. However, partial correlation being zero may not be equivalent to conditional independence under non-Gaussian distributions. In this paper, we propose a statistical inference procedure for partial correlations under the high-dimensional nonparanormal (NPN) model where the observed data are normally distributed after certain monotone transformations. The NPN partial correlation is the partial correlation of the normal transformed data under the NPN model, which is a more general measure of conditional dependence. We estimate the NPN partial correlations by regularized nodewise regression based on the empirical ranks of the original data. A multiple testing procedure is proposed to identify the nonzero NPN partial correlations. The proposed method can be carried out by a simple coordinate descent algorithm for lasso optimization. It is easy-to-implement and computationally more efficient compared to the existing methods for estimating NPN graphical models. Theoretical results are developed to show the asymptotic normality of the proposed estimator and to justify the proposed multiple testing procedure. Numerical simulations and a case study on brain imaging data demonstrate the utility of the proposed procedure and evaluate its performance compared to the existing methods. Data used in preparation of this article were obtained from the Alzheimer's Disease Neuroimaging Initiative (ADNI) database.

97 MATHEMATICS AND COMPUTING↗