Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Charge distribution shape”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
NASA NTRS1982

The sputtering of insulating materials by fast heavy ions

In this paper recent experimental results on sputtering of UF4 and H2O (ice) by fast heavy ions are reviewed. Measurements have been made of the dependence of the sputtering yield on the incident ion type, charge state, and energy. In the case of UF4, the energy spectra of neutral sputtered particles have been obtained as well. There is a clear dependence of the sputtering yield on the electronic part of the stopping power, and the yield is strongly affected by the charge state of the incident ion, which shows that in the near-surface region from which sputtered particles arise, the charge state of the incident ion has not reached equilibrium. The shape of the energy spectra observed for UF4 targets is that expected from a thermal distribution rather than that of the collision cascade form typical of ordinary sputtering. A model of the sputtering-track registration process that has arisen from these data is shown to provide a framework of understanding the close relationship of these observations to the 'plasma desorption' of marcromolecules.

Seiberling, L. E.↗
DOE OSTI2025

Charged-particle multiplicity distributions over a wide pseudorapidity range in p–Pb collisions at $\sqrt{s_{NN}}$ = 5.02 TeV

This paper presents the primary charged-particle multiplicity distributions in proton–lead collisions at a centre-of-mass energy per nucleon–nucleon collision of $\sqrt{s_{NN}}$ = 5.02 TeV. The distributions are reported for non-single diffractive collisions in different pseudorapidity ranges. The measurements are performed using the combined information from the Silicon Pixel Detector and the Forward Multiplicity Detector of ALICE. The multiplicity distributions are parametrised with a double negative binomial distribution function which provides satisfactory descriptions of the distributions for all the studied pseudorapidity intervals. The data are compared to models and analyzed quantitatively, evaluating the first four moments (mean, standard deviation, skewness, and kurtosis). The shape evolution of the measured multiplicity distributions is studied in terms of KNO variables and it is found that none of the considered models reproduces the measurements. This paper also reports on the average charged-particle multiplicity, normalised by the average number of participating nucleon pairs, as a function of the collision energy. The multiplicity results are then compared to measurements made in proton–proton and nucleus–nucleus collisions across a wide range of collision energies.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Particle production as a function of charged-particle flattenicity in 𝑝⁢𝑝 collisions at $\sqrt{s}$ = 13 TeV

This paper reports the first measurement of the transverse momentum (𝑝 T ) spectra of primary charged pions, kaons, (anti)protons, and unidentified particles as a function of the charged-particle flattenicity in pp collisions at $\sqrt{s}$ =13 TeV. Flattenicity is a novel event shape observable that is measured in the pseudorapidity intervals covered by the V0 detector, 2.8 < 𝜂 < 5.1 and −3.7 < 𝜂 < −1.7. According to QCD-inspired phenomenological models, it shows sensitivity to multiparton interactions and is less affected by biases toward larger 𝑝 T due to local multiplicity fluctuations in the V0 acceptance than multiplicity. The analysis is performed in minimum-bias (MB) as well as in high-multiplicity events up to 𝑝 T = 20 GeV/⁢𝑐. The event selection requires at least one charged particle produced in the pseudorapidity interval |𝜂| < 1. The measured 𝑝 T distributions, average 𝑝 T , kaon-to-pion and proton-to-pion particle ratios, presented in this paper, are compared to model calculations using pythia 8 based on color strings and EPOS LHC. The modification of the 𝑝 T -spectral shapes in low-flattenicity events that have large event activity with respect to those measured in MB events develops a pronounced peak at intermediate 𝑝 T (2 < 𝑝 T < 8 GeV/𝑐), and approaches the vicinity of unity at higher 𝑝 T . The results are qualitatively described by pythia, and they show different behavior than those measured as a function of charged-particle multiplicity based on the V0M estimator.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Shell effects in quasi-fission: Reactions forming actinide and superheavy compound nuclei

Quasi-fission occurs in fully damped heavy-ion collisions following a significant mass transfer between the fragments, without formation of a compound nucleus. As in fission, quasi-fission is expected to be affected by quantum effects leading to asymmetric mass splits. Quantum shells stabilising fission fragments with octupole shapes have been invoked as a factor determining the distribution of nucleons between the fragments at scission, explaining the fact that the centroid of the heavy fragment charge distribution is found around Z = 54 protons in fission and quasi-fission of actinides. These shell effects are studied in microscopic studies of quasi-fission. In particular, time-dependent Hartree-Fock (TDHF) calculations have been performed for reactions forming actinides.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1998

Interpretation of Solar Wind Ion Composition Measurements from Ulysses

The ion compositions measured in situ in the solar wind are important since the ion fractions carry information on the plasma conditions in the inner corona. The conditions in the inner corona define the properties of the solar wind plasma flow. Thus, if the ion fraction measurements can be used to unravel some of the plasma parameters in the inner corona, they will provide a valuable contribution to solving the heating and acceleration problem of the solar wind. The ion charge states in the solar wind carry information on electron temperature, electron density and ion flow speed. They are also sensitive to the shape of the electron distribution function. Through carefully modeling the solar wind and calculating the ion fractions predicted for different solar wind conditions, constraints on the electron temperature and ion flow speeds can be placed if the electron density is measured using polarization brightness measurements.

Esser, Ruth↗
DOE OSTI2023

STRESS AND FATIGUE ANALYSIS OF AIR-TO-REFRIGERANT HEAT EXCHANGERS WITH NON-ROUND TUBE SHAPES

Air-to-refrigerant heat exchangers are fundamental components in Heating, Ventilation, Air Conditioning, and Refrigeration (HVAC&R) systems. Recent research has shown that heat exchangers utilizing small hydraulic diameter, non-round, shape-optimized tubes to be promising next design update. Such tubes can improve performance, reduce material and lower refrigerant charge. However, such tube designs may be more susceptible to fatigue issues. Therefore, stress and fatigue analyses of such heat exchangers are critical to ensure structural integrity, reliability, and manufacturability. In the present work, a mechanical model is developed to obtain the stress distributions in heat exchangers with round, ellipse, and airfoil-shaped tubes. The results indicate that under 3.45 MPa internal pressure, the maximum stresses don't surpass 190 MPa.

Zhang, Mingkan↗
NASA NTRS1972

Mass transport contamination study

A theoretical analysis was performed to determine the effects of outgassing and waste dumping on the contamination field around an orbiting spacecraft. The spacecraft was assumed to be spherical in shape with the mass flow emitting uniformly from the spherical surface at a constant rate and in a D'Lambertian spatial distribution. The outflow of gases were assumed to be neutrally charged and of a single species with a molecular weight characteristic of a composite of the actual species involved in the mass flow. The theoretical analysis showed that, for outgassing only, less than 1.5 percent of the outgas products will return to the Skylab spacecraft as a result of intermolecular collisions. When the total mass flow from the spacecraft, including waste dumps and reaction control motor firings, was considered, it was estimated that about 30 percent will return to the spacecraft.

Robertson, S. J.↗
NASA NTRS2008

Damage Tolerance Testing of a NASA TransHab Derivative Woven Inflatable Module

Current options for Lunar habitat architecture include inflatable habitats and airlocks. Inflatable structures can have mass and volume advantages over conventional structures. Inflatable structures are perceived to carry additional risk because they are at a lower Technical Readiness Level (TRL) than conventional metallic structures. One of the risks associated with inflatable structures is understanding the tolerance to component damage and the resulting behavior of the system after the damage is introduced. The Damage Tolerance Test (DTT) is designed to study the structural integrity of an expandable structure during and subsequent to induced damage. The TransHab Project developed an experimental inflatable module developed at Johnson Space Center in the 1990's. The TransHab design was originally envisioned for use in Mars Transits but was also studied as a potential habitat for the International Space Station (ISS). The design of the TransHab module was based on a woven design using an Aramid fabric. Testing of this design demonstrated a high level of predictability and repeatability and good correlation with analytical predictions of stresses and deflections. Based on JSC's experience with the design and analysis of woven inflatable structures, the Damage Tolerance Test article was designed and fabricated using a woven design. The Damage Tolerance Test Article consists of a load bearing restraint layer, a bladder or gas barrier, and a structural metallic core. The test article restraint layer is fabricated from one inch wide Kevlar webbing that is woven in a basket weave pattern. Underneath the structural restraint layer is the bladder or gas barrier. For this test the bladder was required to maintain pressure for testing only and was not representative of a flight design. The bladder and structural restraint layer attach to the structural core of the module at steel bulkheads at each end. The two bulkheads are separated by a 10 foot center tube which provides the structural support for the module when in a non-inflated state as well as resists a portion of the axial load when pressurized. The longitudinal members of the structural restraint layer are attached to the bulkheads using a series of clevises that are bolted to the bulkheads. Strain gages are placed on the clevises that can measure change in load when the structural restraint is inflated. The test module is 88 inches in diameter and 120 inches in height. The objectives of the DTT are to (1) verify the structural integrity of the assembled and pressurized structure when a section of the structural restraint layer is cut by a foreign object, and (2) verify the load distribution of the structural restraint layer during pressurization, before and after the structural restraint layer is severed. For this test, a longitudinal structural restraint strap will be severed using a linear shape charge. The linear shape charge was designed specifically for this application to cut only a single longitudinal strap, while not damaging the bladder. An array of strain gages were located at the bulkhead mounted clevises where the longitudinal restraint layer straps are attached. The DTT article was inflated to 45 psig, 25% of the ultimate design pressure, and one of the one-inch wide longitudinal structural members was severed. Strain gage measurements of loading in an array of longitudinal straps were taken throughout pressurization of the module to 45 psig, before firing of the linear shape charge, and after firing of the shape charge and separation of the strap. During testing not only were the original objectives met but better than expected results occurred. This paper will discuss space inflatable structures, damage tolerance analysis, test results, and applicability to the Lunar architecture.

Edgecombe, John↗
NASA NTRS2012

GEMS X-ray Polarimeter Performance Simulations

The Gravity and Extreme Magnetism Small explorer (GEMS) is an X-ray polarization telescope selected as a NASA small explorer satellite mission. The X-ray Polarimeter on GEMS uses a Time Projection Chamber gas proportional counter to measure the polarization of astrophysical X-rays in the 2-10 keV band by sensing the direction of the track of the primary photoelectron excited by the incident X-ray. We have simulated the expected sensitivity of the polarimeter to polarized X-rays. We use the simulation package Penelope to model the physics of the interaction of the initial photoelectron with the detector gas and to determine the distribution of charge deposited in the detector volume. We then model the charge diffusion in the detector,and produce simulated track images. Within the track reconstruction algorithm we apply cuts on the track shape and focus on the initial photoelectron direction in order to maximize the overall sensitivity of the instrument, using this technique we have predicted instrument modulation factors nu(sub 100) for 100% polarized X-rays ranging from 10% to over 60% across the 2-10 keV X-ray band. We also discuss the simulation program used to develop and model some of the algorithms used for triggering, and energy measurement of events in the polarimeter.

Baumgartner, Wayne H.↗
DOE OSTI2025

Measurement of the Neutron Electromagnetic Form Factors at High Momentum Transfer Using a Polarized 3He Target at Jefferson Lab

Ever since the composite structure of nucleons was uncovered, many experiments have been commissioned with the purpose of extracting parameters that shape our understanding of this structure. Form factors are powerful tools that rest at the forefront of the nucleon structure study. In the non-relativistic limit, they describe the magnetic and charge distributions inside the nucleon. The GEn-II experiment, which ran in Hall A of Jefferson Lab as a part of the SuperBigbite Spectrometer (SBS) program, aimed to extract the ratio of the neutron’s electric to magnetic Sachs form factors (GnE/GnM. This was accomplished using the double polarization technique: scattering polarized electrons with initial energies up to 8.45 GeV off high-density polarized 3He targets. The experiment was a coincidence measurement, with the scattered electrons detected in the BigBite Spectrometer and the recoil nucleons in a hadron calorimeter, a part of the SuperBigbite Spectrometer. The convection-style target cells used in this experiment have a 60cm long target chamber, an increase from the 40 cm cell used in the most recent polarized 3He target experiment. In addition, the cells were optically pumped from two sides of the pumping chamber. These improvements in design and functionality introduced world-record-breaking luminosity. This thesis presents preliminary results of the GEn-II experiment for the following 4-momentum transfer (Q2) values: 2.93, 6.76, and 9.78 GeV2, in addition to target performance and preliminary polarimetry results.

Presley, Hunter [Univ. of Virginia, Charlottesvill↗
DOE OSTI2024

Effect of Molecular Conformation on Intermolecular Interactions and Photovoltaic Performances of Giant Molecule Acceptors

The molecular conformation of giant molecule acceptors (GMAs) plays a significant role in regulating the intermolecular interactions and their photovoltaic performances in organic solar cells (OSCs). For the linear GMA GT‐l, the stronger homo‐molecular interaction causes its aggregation being weakly affected by the donor, thus forming an ordered molecular stacking and proper phase separation in its blend film. The star‐shaped GMA GT‐s‐based blend film shows a dominant hetero‐molecular interaction that suppresses the aggregation of the donor and acceptor, resulting in smaller phase separation and more uniform vertical phase distribution. While for another star‐shaped GMA GTs, the weakest hetero‐molecular interaction causes its blend film to form larger phase separation. Therefore, the GT‐l based OSC with PM6 as donor shows the highest charge mobilities, the fastest charge transfer (CT) process, reduced energy loss and less charge recombination, contributing to a higher power conversion efficiency (PCE) of 19.03%. Comparatively, the PCEs of the OSCs based on GTs and GT‐s are 18.05% and 17.58% respectively. Notably, all the three GMAs based OSCs show excellent thermal stability and long‐term storage stability. This study provides a facile strategy by tuning the linking unit and its connecting mode for designing highly efficient and stable organic photovoltaic materials.

Chemistry↗
NASA NTRS2015

Study of Static Microchannel Plate Saturation Effects for the Fast Plasma Investigation Dual Electron Spectrometers on NASA's Magnetospheric MultiScale Mission

Imaging detecting systems based on microchannel plates (MCPs) are the most common for low energy plasma measurements for both space borne and ground applications. One of the key parameters of these detection systems is the dynamic range of the MCP's response to the input fluxes of charged particles. For most applications the dynamic range of the linear response should be as wide as possible. This is especially true for the Dual Electron Spectrometers (DESs) of the Fast Plasma Investigation (FPI) on NASA's Magnetospheric MultiScale (MMS) mission because a wide range of input fluxes are expected. To make use of the full available dynamic range, it is important to understand the MCP response behavior beyond the linear regime where the MCPs start to saturate. We have performed extensive studies of this during the characterization and calibration of the DES instruments and have identified several saturation effects of the detection system. The MCP itself exhibits saturation when the channels lack the ability to replenish charge sufficiently rapidly. It is found and will be shown that the ground system can significantly impact the correct measurement of this effect. As the MCP starts to saturate, the resulting pulse height distribution (PHD) changes shape and location (with less pulse height values), which leads to truncation of the PHD by the threshold set on the detection system discriminator. Finally, the detection system pulse amplifier exhibits saturation as the input flux drives pulse rates greater than its linear response speed. All of these effects effectively change the dead time of the overall detection system and as a result can affect the quality and interpretation of the flight data. We present results of detection system saturation effects and their interaction with special emphasis on the MCP related effects.

detecting system↗
DOE OSTI2022

Inverse Calculation of Burden Distribution Matrix Using B-spline Model Based PDF control in Blast Furnace Burden Charging Process

The inverse calculation of burden distribution matrix (BDM) is one of the most important challenges in the blast furnace operation in iron-making processes. In general, blast furnace consumes 65% of the total energy for the whole steel-making. Focusing on this practical challenge, this article proposes a new burden distribution spatial model in calculating burden charging process, and develops a B-spline approximation-based probability density function (PDF) control algorithm to assign the expected thickness distribution of burden layer and, thus, develops a new method for the required inverse calculation of BDM. First, a novel method for the thickness distribution of burden layer is given using B-spline model to produce an expected distribution shape subjected to a desired tracking within a specific spatial constraint. Then, according to the coexistence of continuous and bounded discrete variables in BDM, a novel hybrid optimization control method by combining integer programming and PDF tracking is further established for the effective inverse calculation of BDM. Finally, the proposed PDF-based iterative inverse calculation of BDM using B-spline models are tested using various data from industrial examples. Furthermore, the simulation results show that the proposed method is well suited to solve the BDM inverse calculation problem in practice.

42 ENGINEERING↗
DOE OSTI2024

Shape-persistent ladder molecules exhibit nanogap-independent conductance in single-molecule junctions

Molecular electronic devices require precise control over the flow of current in single molecules. However, the electron transport properties of single molecules critically depend on dynamic molecular conformations in nanoscale junctions. Here, in this work, we report a unique strategy for controlling molecular conductance using shape-persistent molecules. Chemically diverse, charged ladder molecules, synthesized via a one-pot multicomponent ladderization strategy, show a molecular conductance (d[log( G/G 0 )]/d x ≈ -0.1 nm -1 ) that is nearly independent of junction displacement, in stark contrast to the nanogap-dependent conductance (d[log( G/G 0 )]/d x ≈ -7 nm -1 ) observed for non-ladder analogues. Ladder molecules show an unusually narrow distribution of molecular conductance during dynamic junction displacement, which is attributed to the shape-persistent backbone and restricted rotation of terminal anchor groups. These principles are further extended to a butterfly-like molecule, thereby demonstrating the strategy's generality for achieving gap-independent conductance. Overall, our work provides important avenues for controlling molecular conductance using shape-persistent molecules. Achieving robust and controllable conductance in single-molecule junctions is challenging due to the dynamic nature of molecular conformations that fluctuate over operational timescales. A strategy using shape-persistent molecules has now been developed that demonstrates nearly junction-displacement-independent conductance, providing a stable solution for single-molecule electronic properties.

molecular electronics↗
DOE OSTI2022

Noncoordinating Secondary Sphere Ion Modulates Supramolecular Clustering of Lanthanides

The role of counterions in molecular recognition of lanthanides is underexplored, especially when they exhibit weak interactions with the metal cations. Here, we report a complementary and comprehensive investigation integrating theoretical calculations with X-ray absorption fine structure spectroscopy, dynamic light scattering, and small-angle X-ray scattering to reveal atomic-scale structural features beyond the immediate coordination sphere of a system used for rare-earth element separations. Here, our results indicate the formation of an unusual T-shaped outer-sphere lanthanide complex, containing two ligands and two nitrate ions in the first coordination sphere, whereas the third nitrate is weakly coordinated and resides in the second shell. This unique structural arrangement causes inhomogeneous charge distribution, leading to self-assembly of the complexes into larger nanoclusters through sterically directed electrostatic interactions in the nonpolar medium. Our findings point to the importance of "noncoordinating" anions in defining the degree of supramolecular aggregation and ion cluster assembly.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

System Assessment of a High Power 3-U CubeSat

The Advanced eLectrical Bus (ALBus) CubeSat project is a technology demonstration mission of a 3-UCubeSat with an advanced, digitally controlled electrical power system capability and novel use of Shape Memory Alloy (SMA) technology for reliable deployable solar array mechanisms. The objective of the project is to, through an on orbit demonstration, advance the state of power management and distribution (PMAD) capabilities to enable future missions requiring higher power, flexible and reliable power systems. The goals of the mission include demonstration of: 100 Watt distribution to a target electrical load, efficient battery charging in the orbital environment, flexible power system distribution interfaces, adaptation of power system control on orbit, and reliable deployment of solar arrays and antennas utilizing re-settable SMA mechanisms. The power distribution function of the ALBus PMAD system is unique in the total power to target load capability of 100 W, the flexibility to support centralized or point-to-load regulation and ability to respond to fast transient power requirements. Power will be distributed from batteries at 14.8 V, 6.5 A to provide 100 W of power directly to a load. The deployable solar arrays utilize NASA Glenn Research Center superelastic and activated Nitinol(Nickel-Titanium alloy) Shape Memory Alloy (SMA) technology for hinges and a retention and release mechanism. The deployable solar array hinge design features utilization of the SMA material properties for dual purpose. The hinge uses the shape memory properties of the SMA to provide the spring force to deploy the arrays. The electrical conductivity properties of the SMA also enables the design to provide clean conduits for power transfer from the deployable arrays to the power management system. This eliminates the need for electrical harnesses between the arrays and the PMAD system in the ALBus system design. The uniqueness of the SMA retention and release mechanism design is the ability to reset the mechanism, allowing functional tests of the mechanisms prior to flight with no degradation of performance. The project is currently in preparation at the NASA Glenn Research Center for a launch in late calendar year of 2017. The 100 Watt power distribution and dual purpose, re-settable SMA mechanisms introduced several system level challenges due to the physical constraints in volume, mass and surface area of 3-U CubeSats. Several trade studies and design cycles have been completed to develop a system which supports the project objectives. This paper is a report on the results of the system level trade studies and assessments. The results include assessment of options for thermal control of 100 Watts of power dissipation, data from system analyses and engineering development tests, limitations of the 3-U system and extensibility to larger scale CubeSat missions.

CubeSat↗
DOE OSTI2025

Scalability and Effectiveness of Smart Charge Management

The rise in electric vehicle (EV) adoption presents growing challenges for power grids, particularly from simultaneous residential charging, which can cause voltage fluctuations and increase feeder peak loads. Baltimore Gas and Electric (BGE), with support from the U.S. Department of Energy, initiated a pilot program to evaluate managed residential EV charging through Smart Charge Management (SCM). This study analyzes real-world charging behavior data from the pilot and feeder-level base loads from BGE to simulate residential charging scenarios through 2035 across the Washington, DC–Baltimore region. Grid impacts under unmanaged charging are compared to three SCM strategies: TOU-immediate, TOU-distributed, and Load Balancing. Results show that the magnitude of peak reduction is highly feeder-dependent. Some feeders achieve reductions of more than 40% at high enrollment levels, while others show improvements closer to 10–15%. This heterogeneity reflects differences in baseline feeder load shapes, EV penetration, and plug-in behavior across customers. Results also highlight trade-offs between shifting load away from peak periods and minimizing secondary demand peaks, offering practical insights for future utility program design.

Electric vehicle↗
DOE OSTI2023

Effective Landau-type model of a Hf x Zr 1 - x O 2 -graphene nanostructure

Here, to describe the charge-polarization coupling in the nanostructure formed by a thin Hf x Zr 1-x O 2 film with a single-layer graphene as a top electrode, we develop the “effective” Landau-Ginzburg-Devonshire model. This approach is based on the parametrization of the Landau expansion coefficients for the polar (FE) and antipolar (AFE) orderings in thin Hf 1-x Zr x O 2 films from a limited number of polarization-field curves and hysteresis loops. The Landau expansion coefficients are nonlinearly dependent on the film thickness h and Zr/[Hf+Zr] ratio x, in contrast to h-independent and linearly x-dependent expansion coefficients of a classical Landau energy. We explain the dependence of the Landau expansion coefficients by the strong nonmonotonic dependence of the Hf 1-x Zr x O 2 film polar properties on the film thickness, grain size and surface energy. The proposed Landau free energy with five “effective” expansion coefficients, which are interpolation functions of x and h, describes the continuous transformation of polarization dependences on applied electric field and hysteresis loop shapes induced by the changes of x and h in the range 0 < x < 1 and 5 nm < h < 35 nm. Using the effective free energy, we demonstrated that the polarization of Hf 1-x Zr x O 2 film influences strongly on the graphene conductivity, and the full correlation between the distribution of polarization and charge carriers in graphene is revealed. In accordance with our modeling, the polarization of the (5 – 25) nm thick Hf 1-x Zr x O 2 films, which are in the ferroelectric-like or antiferroelectric-like states for the chemical compositions 0.35 ≤ x ≤ 0.95, determine the concentration of carriers in graphene and can control its field dependence. The result can be promising for creation of next generation Si-compatible nonvolatile memories and graphene-ferroelectric FETs, because the working voltages applied to the Hf 1-x Zr x O 2 film (which acts as a gate) can be relatively low (less than 2 V). These low voltages are sufficient to induce the pronounced hysteresis of ferroelectric polarization in the Hf 1-x Zr x O 2 gate, which, due to the strong electric coupling, induces the hysteresis of the graphene charge.

36 MATERIALS SCIENCE↗