Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Density matrix 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 973 records · Page 54

R-matrix calculations for opacities: I. Methodology and computations

Abstract An extended version of the R -matrix methodology is presented for calculation of radiative parameters for improved plasma opacities. Contrast and comparisons with existing methods primarily relying on the distorted wave approximation are discussed to verify accuracy and resolve outstanding issues, particularly with reference to the opacity project (OP). Among the improvements incorporated are: (i) large-scale Breit–Pauli R -matrix calculations for complex atomic systems including fine structure, (ii) convergent close coupling wave function expansions for the ( e + ion) system to compute oscillator strengths and photoionization cross sections, (iii) open and closed shell iron ions of interest in astrophysics and experiments, (iv) a treatment for plasma broadening of autoionizing resonances as function of energy-temperature-density dependent cross sections, (v) a ‘top-up’ procedure to compare convergence with R -matrix calculations for highly excited levels, and (vi) spectroscopic identification of resonances and bound ( e + ion) levels. The present R -matrix monochromatic opacity spectra are fundamentally different from OP and lead to enhanced Rosseland and Planck mean opacities. An outline of the work reported in other papers in this series and those in progress is presented. Based on the present re-examination of the OP work, opacities of heavy elements might require revisions in high temperature-density plasma sources.

Pradhan, A. K. (ORCID:0000000187753643)↗

Holographic View of Mixed-State Symmetry-Protected Topological Phases in Open Quantum Systems

We establish a holographic duality between d -dimensional mixed-state symmetry-protected topological (mSPT) phases and ( d + 1 ) -dimensional subsystem symmetry-protected topological (SSPT) states. Specifically, we show that the reduced density matrix of the boundary layer of a ( d + 1 ) -dimensional SSPT state with subsystem symmetry S and global symmetry G corresponds to a d -dimensional mSPT phase with strong S and weak G symmetries. Conversely, we demonstrate that the wave function of an SSPT state can be constructed by replicating the density matrix of the corresponding lower-dimensional mSPT phase. This mapping links the density matrix in lower dimensions to the entanglement properties of higher-dimensional wave functions, providing an approach for analyzing nonlinear quantities and quantum information metrics in mixed-state systems. Our duality offers a new perspective for studying intrinsic mSPT phases that are unique to open quantum systems, without pure-state analogues. We show that strange correlators and twisted Rényi- N correlators can diagnose these nontrivial phases and we explore their connection to strange correlators in pure-state SSPT phases. Furthermore, we discuss several implications of this holographic duality, including a method for preparing intrinsic mSPT phases through the duality. Published by the American Physical Society 2025

Sun, Shijun (ORCID:0000000168880030)↗

Simulating Electron Dynamics of Complex Molecules with Time-Dependent Complete Active Space Configuration Interaction

In this paper, time-dependent electronic structure methods are growing in popularity as tools for modeling ultrafast and/or nonlinear processes, for computing spectra, and as the electronic structure component of mean-field molecular dynamics simulations. Time-dependent configuration interaction (TD-CI) offers several advantages over the widely used real-time time-dependent density functional theory: namely, that it correctly models Rabi oscillations; it offers a spin-pure description of open-shell systems; and a hierarchy of TD-CI methods can be defined that systematically approach the exact solution of the time-dependent Schrodinger equation (TDSE). In this work, we present a novel TD-CI approach that extends TD-CI to large complete active-space configuration expansions. Such extension is enabled by use of a direct configuration interaction approach that eliminates the need to explicitly build, store, or diagonalize the Hamiltonian matrix. Graphics processing unit (GPU) acceleration enables fast solution of the TDSE even for large active spaces—up to 12 electrons in 12 orbitals (853776 determinants) in this work. A symplectic split operator propagator yields long-time norm conservation. We demonstrate the applicability of our approach by computing the response of a large molecule with a strongly correlated ground state, decacene (C 42 H 24 ), to various pulses (δ-function, transform limited, chirped). Our simulations predict that chirped pulses can be used to induce dipole-forbidden transitions. Simulations of decacene using the 6-31G(d) basis set and a 12 electrons/12 orbitals active space took 20.1 h to propagate for 100 fs with a 1 attosecond time step on a single NVIDIA K40 GPU. Convergence with respect to time step is found to depend on the property being computed and the chosen active space.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Thermally induced damage in composite space structure: Predictive methodology and experimental correlation

A general analysis method is presented to predict matrix cracks in all plies of a composite laminate, and resulting degraded laminate properties, as functions of temperature or thermal cycles. A shear lag solution of the stresses in the vicinity of cracks and a fracture mechanics crack formation criteria are used to predict cracks. Damage is modeled incrementally, which allows the inclusion of the effects of temperature dependent material properties and softening of the laminate due to previous cracking. The analysis is incorporated into an easy-to-use computer program. The analysis is correlated with experimentally measured crack densities in a variety of laminates exposed to monotonically decreasing temperatures. Crack densities are measured at the edges of specimens by microscopic inspection, and throughout the specimen volumes by x-ray and sanding down of the edges. Correlation between the analytical results and the crack densities in the interiors of the specimens was quite good. Crack densities measured at specimen edges did not agree with internal crack densities (or analyses) in some cases. A free-edge stress analysis clarified the reasons for these discrepancies.

Park, Cecelia↗

Residual entropy and waterlike anomalies in the repulsive one dimensional lattice gas

The thermodynamics and kinetics of the one dimensional lattice gas with repulsive interaction are investigated using transfer matrix technique and Monte Carlo simulations. This simple model is shown to exhibit waterlike anomalies in density, thermal expansion coefficient, and self-diffusion. An unified description for the thermodynamic anomalies in this model is achieved based on the ground state residual entropy which appears in the model due to mixing entropy in a ground state phase transition.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Influence of alloy disorder scattering on the hole mobility of SiGe nanowires

In this work, we analyze the influence of the alloy disorder (AD) scattering on the low-field hole mobility of Si{sub 1-x}Ge{sub x} nanowires (NWs). To do it, the electrostatic description is achieved through a self-consistent solution of the Poisson equation and the six-band k⋅p method in the cross section of the NW. The momentum relaxation time approximation is used to calculate the hole mobility, including alloy disorder and phonon scattering mechanisms, and the use of approximations to calculate the overlap integrals for the scattering matrix elements is discussed. We study the influence of the alloy disorder scattering on the total mobility compared to the phonon contribution, for different values of the AD scattering parameter proposed in the literature, and analyze the performance of SiGe NWs as a function of the Ge molar fraction for both low and high inversion charge densities.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Towards Highly Scalable Ab Initio Molecular Dynamics (AIMD) Simulations on the Intel Knights Landing Manycore Processor

The Ab Initio Molecular Dynamics (AIMD) method allows scientists to treat the dynamics of molecular and condensed phase systems while retaining a first-principles-based description of their interactions. This extremely important method has tremendous computational requirements, because the electronic Schr¨odinger equation, approximated using Kohn-Sham Density Functional Theory (DFT), is solved at every time step. With the advent of manycore architectures, application developers have a significant amount of processing power within each compute node that can only be exploited through massive parallelism. A compute intensive application such as AIMD forms a good candidate to leverage this processing power. In this paper, we focus on adding thread level parallelism to the plane wave DFT methodology implemented in NWChem. Through a careful optimization of tall-skinny matrix products, which are at the heart of the Lagrange multiplier and nonlocal pseudopotential kernels, as well as 3D FFTs, our OpenMP implementation delivers excellent strong scaling on the latest Intel Knights Landing (KNL) processor. We assess the efficiency of our Lagrange multiplier kernels by building a Roofline model of the platform, and verify that our implementation is close to the roofline for various problem sizes. Finally, we present strong scaling results on the complete AIMD simulation for a 64 water molecules test case, that scales up to all 68 cores of the Knights Landing processor.

Environmental Molecular Sciences Laboratory↗

Improved Thermoplastic/Iron-Particle Transformer Cores

A method of fabricating improved transformer cores from composites of thermoplastic matrices and iron-particles has been invented. Relative to commercially available laminated-iron-alloy transformer cores, the cores fabricated by this method weigh less and are less expensive. Relative to prior polymer-matrix/ iron-particle composite-material transformer cores, the cores fabricated by this method can be made mechanically stronger and more magnetically permeable. In addition, whereas some prior cores have exhibited significant eddy-current losses, the cores fabricated by this method exhibit very small eddy-current losses. The cores made by this method can be expected to be attractive for use in diverse applications, including high-signal-to-noise transformers, stepping motors, and high-frequency ignition coils. The present method is a product of an experimental study of the relationships among fabrication conditions, final densities of iron particles, and mechanical and electromagnetic properties of fabricated cores. Among the fabrication conditions investigated were molding pressures (83, 104, and 131 MPa), and molding temperatures (250, 300, and 350 C). Each block of core material was made by uniaxial-compression molding, at the applicable pressure/temperature combination, of a mixture of 2 weight percent of LaRC (or equivalent high-temperature soluble thermoplastic adhesive) with 98 weight percent of approximately spherical iron particles having diameters in the micron range. Each molded block was cut into square cross-section rods that were used as core specimens in mechanical and electromagnetic tests. Some of the core specimens were annealed at 900 C and cooled slowly before testing. For comparison, a low-carbon-steel core was also tested. The results of the tests showed that density, hardness, and rupture strength generally increased with molding pressure and temperature, though the correlation was rather weak. The weakness of the correlation was attributed to the pores in the specimens. The maximum relative permeabilities of cores made without annealing ranged from 30 to 110, while those of cores made with annealing ranged from 900 to 1,400. However, the greater permeabilities of the annealed specimens were not associated with noticeably greater densities. The major practical result of the investigation was the discovery of an optimum distribution of iron-particle sizes: It was found that eddy-current losses in the molded cores were minimized by using 100 mesh (corresponding to particles with diameters less than or equal to 100 m) iron particles. The effect of optimization of particle sizes on eddy-current losses is depicted in the figure.

Wincheski, Russell A.↗

Renormalization of quantum discord and Bell nonlocality in the XXZ model with Dzyaloshinskii–Moriya interaction

In this paper, we have investigated the dynamical behaviors of the two important quantum correlation witnesses, i.e. geometric quantum discord (GQD) and Bell–CHSH inequality in the XXZ model with DM interaction by employing the quantum renormalization group (QRG) method. The results have shown that the anisotropy suppresses the quantum correlations while the DM interaction can enhance them. Meanwhile, using the QRG method we have studied the quantum phase transition of GQD and obtained two saturated values, which are associated with two different phases: spin-fluid phase and the Néel phase. It is worth mentioning that the block–block correlation is not strong enough to violate the Bell–CHSH inequality in the whole iteration steps. Moreover, the nonanalytic phenomenon and scaling behavior of Bell inequality are discussed in detail. As a byproduct, the conjecture that the exact lower and upper bounds of Bell inequality versus GQD can always be established for this spin system although the given density matrix is a general X state.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Generalized Sampling and Preconditioning Scheme for Sparse Approximation of Polynomial Chaos Expansions

We propose an algorithm for recovering sparse orthogonal polynomial expansions via collocation. A standard sampling approach for recovering sparse polynomials uses Monte Carlo sampling, from the density of orthogonality, which results in poor function recovery when the polynomial degree is high. Our proposed approach aims to mitigate this limitation by sampling with respect to the weighted equilibrium measure of the parametric domain and subsequently solves a preconditioned $\ell^1$-minimization problem, where the weights of the diagonal preconditioning matrix are given by evaluations of the Christoffel function. Our algorithm can be applied to a wide class of orthogonal polynomial families on bounded and unbounded domains, including all classical families. We present theoretical analysis to motivate the algorithm and numerical results that show our method is superior to standard Monte Carlo methods in many situations of interest. In conclusion, numerical examples are also provided to demonstrate that our proposed algorithm leads to comparable or improved accuracy even when compared with Legendre- and Hermite-specific algorithms.

97 MATHEMATICS AND COMPUTING↗

Magnetic properties of γ-Fe 2 O 3 nanoparticles in a porous SiO 2 shell for drug delivery

A method is presented for synthesizing core–shell nanoparticles with a magnetic core and a porous shell suitable for drug delivery and other medical applications. The core contains multiple γ -Fe 2 O 3 nanoparticles (~15 nm) enclosed in a SiO 2 (~100–200 nm) matrix using either methyl (denoted TMOS- γ -Fe 2 O 3 ) or ethyl (TEOS- γ -Fe 2 O 3 ) template groups. Low-temperature Mössbauer spectroscopy showed that the magnetic nanoparticles have the maghemite structure, γ -Fe 2 O 3 , with all the vacancies in the octahedral sites. Saturation magnetization measurements revealed that the density of γ -Fe 2 O 3 was greater in the TMOS- γ -Fe 2 O 3 nanoparticles than TEOS- γ -Fe 2 O 3 nanoparticles, presumably because of the smaller methyl group. Magnetization measurements showed that the blocking temperature is around room temperature for the TMOS- γ -Fe 2 O 3 and around 250 K for the TEOS- γ -Fe 2 O 3 . Three dimensional topography analysis shows clearly that the magnetic nanoparticles are not only at the surface but have penetrated deep in the silica to form the core–shell structure.

Physics↗

Analytical derivatives of the individual state energies in ensemble density functional theory method. I. General formalism

The state-averaged (SA) spin restricted ensemble referenced Kohn-Sham (REKS) method and its state interaction (SI) extension, SI-SA-REKS, enable one to describe correctly the shape of the ground and excited potential energy surfaces of molecules undergoing bond breaking/bond formation reactions including features such as conical intersections crucial for theoretical modeling of non-adiabatic reactions. Until recently, application of the SA-REKS and SI-SA-REKS methods to modeling the dynamics of such reactions was obstructed due to the lack of the analytical energy derivatives. Here, the analytical derivatives of the individual SA-REKS and SI-SA-REKS energies are derived. The final analytic gradient expressions are formulated entirely in terms of traces of matrix products and are presented in the form convenient for implementation in the traditional quantum chemical codes employing basis set expansions of the molecular orbitals. Finally, we will describe the implementation and benchmarking of the derived formalism in a subsequent article of this series.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Metal Foam Analysis: Improving Sandwich Structure Technology for Engine Fan and Propeller Blades

The Life Prediction Branch of the NASA Glenn Research Center is searching for ways to construct aircraft and rotorcraft engine fan and propeller blades that are lighter and less costly. One possible design is to create a sandwich structure composed of two metal faces sheets and a metal foam core. The face sheets would carry the bending loads and the foam core would have to resist the transverse shear loads. Metal foam is ideal because of its low density and energy absorption capabilities, making the structure lighter, yet still stiff. The material chosen for the face sheets and core was 17-4PH stainless steel, which is easy to make and has appealing mechanical properties. This material can be made inexpensively compared to titanium and polymer matrix composites, the two current fan blade alternatives. Initial tests were performed on design models, including vibration and stress analysis. These tests revealed that the design is competitive with existing designs; however, some problems were apparent that must be addressed before it can be implemented in new technology. The foam did not hold up as well as expected under stress. This could be due to a number of issues, but was most likely a result of a large number of pores within the steel that weakened the structure. The brazing between the face sheets and the foam was also identified as a concern. The braze did not hold up well under shear stress causing the foam to break away from the face sheets. My role in this project was to analyze different options for improving the design. I primarily spent my time examining various foam samples, created with different sintering conditions, to see which exhibited the most favorable characteristics for our purpose. Methods of analysis that I employed included examining strut integrity under a microscope, counting the number of cells per inch, measuring the density, testing the microhardness, and testing the strength under compression. Shear testing will also be done to examine the strengths of different types of brazes.

Fedor, Jessica L.↗

Liquid-Crystal Thermosets, a New Generation of High-Performance Liquid-Crystal Polymers

One of the major challenges for NASA's next-generation reusable-launch-vehicle (RLV) program is the design of a cryogenic lightweight composite fuel tank. Potential matrix resin systems need to exhibit a low coefficient of thermal expansion (CTE), good mechanical strength, and excellent barrier properties at cryogenic temperatures under load. In addition, the resin system needs to be processable by a variety of non-autoclavable techniques, such as vacuum-bag curing, resin-transfer molding (RTM), vacuum-assisted resin-transfer molding (VaRTM), resin-film infusion (RFI), pultrusion, and advanced tow placement (ATP). To meet these requirements, the Advanced Materials and Processing Branch (AMPB) at NASA Langley Research Center developed a new family of wholly aromatic liquid-crystal oligomers that can be processed and thermally cross-linked while maintaining their liquid-crystal order. All the monomers were polymerized in the presence of a cross-linkable unit by use of an environmentally benign melt-condensation technique. This method does not require hazardous solvents, and the only side product is acetic acid. The final product can be obtained as a powder or granulate and has an infinite shelf life. The obtained oligomers melt into a nematic phase and do not exhibit isotropization temperatures greater than the temperatures of decomposition (Ti > T(sub dec)). Three aromatic formulations were designed and tested and included esters, ester-amides, and ester-imides. One of the major advantages of this invention, named LaRC-LCR or Langley Research Center-Liquid Crystal Resin, is the ability to control a variety of resin characteristics, such as melting temperature, viscosity, and the cross-link density of the final part. Depending on the formulation, oligomers can be prepared with melt viscosities in the range of 10-10,000 poise (100 rad/s), which can easily be melt-processed using a variety of composite-processing techniques. This capability provides NASA with custom-made matrix resins that meet the required processing conditions for the fabrication of textile composites. Once the resin is in place, the temperature is raised to 375 C and the oligomers are cross-linked into a high-glass-transition-temperature (Tg) nematic network without releasing volatiles. The mechanical properties of the fully crosslinked, composite articles are comparable to typical composites based on commercially available epoxy resins.

Dingemans, Theo↗

Spontaneous Raman–LIF–CO–OH measurements of species concentration in turbulent spray flames

This paper presents new measurements of species concentrations, temperature and mixture fraction in selected regions of a turbulent ethanol spray flame. The line-Raman–LIF–COsingle bondOH setup developed at the Sandia's Combustion Research Facility is utilised to probe regions of a spray flame where laser breakdown of liquid droplets is avoided and the remaining interferences can be corrected. The spray flame is stabilised on the piloted Sydney needle spray burner, where axial translation of the liquid injecting needle in the air-blast stream can transition the spray from dilute to dense. The solution to obtaining successful measurements is found to be multifaceted and includes: the appropriate selection of flame conditions; high sensitivity of the Raman detection system permitting reduced laser energies; development of a pre-processing algorithm to reject strong droplet interferences; and application of the hybrid matrix inversion method combined with wavelet denoising to account for interference corrections and noise at the very low signal levels obtained. Unique and necessary for the successful measurements reported in this paper, a pre-processing algorithm is outlined that removes data points corrupted with strong interferences from droplets. These interferences arise from a range of sources, but the most intense are due to the laser interaction with surrounding mist or liquid fragments, such that measurements near the jet centreline are corrupted and hence discarded. Reliable measurements of mixture fraction, temperature obtained from the sum of the species number densities, and species mole fractions are reported for regions in the flames sufficiently far from the centreline. The paper demonstrates the feasibility of the judicious use of Raman scattering in turbulent spray flames, the results of which will be extremely useful for validating numerical simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The propagation characteristics of electromagnetic waves through plasma in the near-field region of low-frequency loop antenna

A high-speed vehicle flying through the atmosphere between 100 and 20 km may suffer from a “communication blackout.” In this paper, a low frequency system with an on-board loop antenna to receive signals is presented as a potential blackout mitigation method. Because the plasma sheath is in the near-field region of the loop antenna, the traditional scattering matrix method that is developed for the far-field region may overestimate the electromagnetic (EM) wave's attenuation. To estimate the EM wave's attenuation in the near-field region, EM interference (EMI) shielding theory is introduced. Experiments are conducted, and the results verify the EMI shielding theory's effectiveness. Simulations are also conducted with different plasma parameters, and the results obtained show that the EM wave's attenuation in the near-field region is far below than that in the far-field region. The EM wave's attenuation increases with the increase in electron density and decreases with the increase in collision frequency. The higher the frequency, the larger is the EM wave's attenuation. During the entire re-entry phase of a RAM-C module, the EM wave's attenuations are below 10 dB for EM waves with a frequency of 1 MHz and below 1 dB for EM waves with a frequency of 100 kHz. Therefore, the low frequency systems (e.g., Loran-C) may provide a way to transmit some key information to high-speed vehicles even during the communication “blackout” period.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Structure and Sulfur: Tuning the Viscoelastic and Surface Properties of Natural Keratin Fibers

Natural keratin fibers, such as wool, possess a complex hierarchical structure that governs their mechanical properties and surface energy. However, the extent to which these characteristics are influenced by combined contributions of structural variations (e.g., fiber diameter, intermediate filament (IF) packing) and chemical composition (e.g., disulfide bond density) remains poorly understood. In this study, we investigate wool fibers from five sheep breeds (Merino, Polwarth, Cheviot, Eider, and Devon) to elucidate how these factors influence viscoelasticity and surface interactions. Using a multimodal approach integrating interfacial and bulk characterization methods, including inverse gas chromatography (IGC), atomic force microscopy-infrared spectroscopy (AFM-IR), X-ray photoelectron spectroscopy (XPS), uniaxial tensile testing, and synchrotron small-angle X-ray scattering (SAXS), we show that the nanometer-thick 18-methyleicosanoic acid (18-MEA) layer is consistently present across all wool types and plays a key role in governing hydrophobicity and surface heterogeneity. A controlled isothermal treatment at 200 °C, designed to cleave disulfide bonds, results in a nearly 40% reduction in specific surface area across all fiber types, accompanied by a significant decrease in tensile strength and 80% reduction in elongation at break for Merino and Devon wool, but limited influence on the mechanical properties of Eider fibers. Furthermore, rate-dependent tensile testing within the elastic regime reveals distinct viscoelastic responses among the fiber types, suggesting that the sulfur-rich protein matrix surrounding IFs and its structure contribute actively to stress partitioning. Altogether, when combined with conclusions from SAXS measurements of IF spacing, our work offers compelling insights into the role of the keratin-associated protein (KAP) matrix in shaping wool fiber mechanics. Differences in mechanical behavior among wool types, despite similar IF spacing or sulfur content, highlight the importance of matrix composition and cross-linking density, suggesting that the molecular architecture of the KAP network may be a dominant factor in determining fiber performance.

X-ray scattering↗

Ultrasensitive Label-free Electronic Chip for DNA Analysis Using Carbon Nanotube Nanoelectrode Arrays

There is a strong need for faster, cheaper, and simpler methods for nucleic acid analysis in today s clinical tests. Nanotechnologies can potentially provide solutions to these requirements by integrating nanomaterials with biofunctionalities. Dramatic improvement in the sensitivity and multiplexing can be achieved through the high-degree miniaturization. Here, we present our study in the development of an ultrasensitive label-free electronic chip for DNA/RNA analysis based on carbon nanotube nanoelectrode arrays. A reliable nanoelectrode array based on vertically aligned multi-walled carbon nanotubes (MWNTs) embedded in a SiO2 matrix is fabricated using a bottom-up approach. Characteristic nanoelectrode behavior is observed with a low-density MWNT nanoelectrode array in measuring both the bulk and surface immobilized redox species. The open-end of MWNTs are found to present similar properties as graphite edge-plane electrodes, with a wide potential window, flexible chemical functionalities, and good biocompatibility. A BRCA1 related oligonucleotide probe with 18 bases is covalently functionalized at the open ends of the MWNTs and specifically hybridized with an oligonucleotide target as well as a PCR amplicon. The guanine bases in the target molecules are employed as the signal moieties for the electrochemical measurements. Ru(bpy)3(2+) mediator is used to further amplify the guanine oxidation signal. This technique has been employed for direct electrochemical detection of label-free PCR amplicon through specific hybridization with the BRCAl probe. The detection limit is estimated to be less than approximately 1000 DNA molecules, approaching the limit of the sensitivity by laser-based fluorescence techniques in DNA microarray. This system provides a general electronic platform for rapid molecular diagnostics in applications requiring ultrahigh sensitivity, high-degree of miniaturization, simple sample preparation, and low- cost operation.

Li, Jun↗