Search NASA⌕ Search

SEARCH · Search NASA

Results for “First Principles Calculations”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Effects of aluminum content on thermoelectric performance of Al x CoCrFeNi high-entropy alloys

Introducing a non-regular distribution in the mass and bonding by including distinctly different elements can reduce the phonon transport even within structurally well-ordered materials. These distributions are a quality of all high-entropy alloys (HEAs), however the inclusion of aluminum in Al x CoCrFeNi is particularly impactful due to the large mismatch in atomic mass with other components. The resultant low phonon conductivity is a requirement for high thermoelectric performance, motivating the investigation of the effects of Al content on phonon transport as well as other thermoelectric properties. This work examines the phonon and electron transport and thermoelectric conversion properties with various Al contents (0 ≤ x Al ≤ 2) in this Cantor alloy using first-principles calculations, molecular dynamics, and semi-classical Boltzmann transport theory. The calculated phonon density of states and thermoelectric properties present reasonable agreements with experiments, including neutron scattering. A large reduction of phonon conductivity (k L ) is observed even with low x Al s, which we attribute to effective phonon scatterings by the large mass mismatch. However, its temperature dependence is not significant, demonstrating a minor contribution of interphonon scattering. In contrast, electrical conductivity (σ) and Seebeck coefficient (S) increase with temperature at higher x Al s with body-centered cubic structures. Therefore, the thermoelectric figure of merit (ZT) of Al x CoCrFeNi HEAs is enhanced by increasing the Al content mainly due to the increase of the thermoelectric power factor (σS 2 ) at high temperatures, while at low temperatures the phonon-scattering enhancement by mass mismatch is also important.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Structural, electronic, and magnetic properties of CoFeVGe-based compounds: Experiment and theory

We have carried out a combined theoretical and experimental investigation of both stoichiometric and nonstoichiometric CoFeVGe alloys. In particular, we have investigated CoFeVGe, Co 1.25 Fe 0.75 VGe, Co 0.75 Fe 1.25 VGe, and CoFe 0.75 VGe bulk alloys. Our first principles calculations suggest that all four alloys show ferromagnetic order, where CoFeVGe, Co 1.25 Fe 0.75 VGe, and Co 0.75 Fe 1.25 VGe are highly spin polarized with spin polarization values of over 80%. However, the spin polarization value of CoFe 0.75 VGe is only about 60%. We have synthesized all four samples using arc melting and high-vacuum annealing at 600 °C for 48 hours. The room temperature x-ray diffraction of these samples exhibits a cubic crystal structure with disorder. All the samples show single magnetic transitions at their Curie temperatures, where the Curie temperature and high field (3T) magnetization are 288 K and 42 emu/g; 305 K and 1.5 emu/g; 238 K and 39 emu/g; and 306 K and 35 emu/g for CoFeVGe, Co 1.25 Fe 0.75 VGe, Co 0.75 Fe 1.25 VGe, and CoFe 0.75 VGe, respectively.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Elastic Properties and Deformation Mechanisms in the van der Waals Single-Crystalline Indium Selenide

The bulk van der Waals (vdW) single-crystalline indium selenide (β-InSe) displays exceptional plasticity in a layered crystalline form at both micro- and macroscale. However, the nanoscale origin of plasticity remains unclear. Herein, an atomic-level study on the deformation mechanisms of InSe by using first-principles calculations is reported. Remarkable anisotropic elasticity is observed in the vdW InSe layered crystal, and the stiffness is dramatically softened because of the vdW gap in the layered structure. The simulations capture the distinct fracture modes in the uniaxial tensile deformation, depending on the loading directions—brittle fracture in the [100] and [110] directions while ductile failure in the [001] direction. The InSe layered crystal structure exhibits superplastic deformability under uniaxial compression. Different transition pathways, including interlayer tangling, amorphization, and cross-linking, are tracked along respective deformation directions. The unprecedented plasticity of InSe layered crystals can be attributed to the phase transition coupled with interlayer gliding and cross-layer dislocation slipping. Here, this study deepens our understanding of the deformation mechanisms of layered materials at the atomistic level and provides insights into tailoring material properties for low-dimensional material design based on its deformation mechanisms.

2D materials↗
DOE OSTI2024

On the origin of magnetic anisotropy in 4f-free ferromagnets based on the CaCu 5 structure *

Using first-principles calculations, we investigate the origin of magnetocrystalline anisotropy in a series of 4f-electron-free intermetallics with CaCu 5 -based structures: YCo 5 , YCo 4 B, and Y 3 Co 13 B 2 . The electronic structure of these compounds is characterized by a set of narrow 3d bands near the Fermi level. In YCo 5 the easy-axis anisotropy originates primarily in the spin–orbit coupling-induced mixing of the electronic states with Co d x 2 -y 2 and d xy character. The analysis of k-resolved anisotropy shows that positive contributions accumulate from the entire Brillouin zone but are particularly large near the k z = 0 plane. The analysis of the single-site and two-site terms reveals a large positive single-site contribution to the magnetocrystalline anisotropy from the Co atoms on the honeycomb sublattice, along with two-site contributions from both honeycomb and kagome sublattices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

High-throughput studies of novel magnetic materials in borides

Borides are a versatile material family with various properties for valuable applications. Conventional magnetism, such as ferromagnetism and antiferromagnetism in borides, have been extensively studied. However, research on unconventional magnetism in borides where quantum effects are dominant is scarce. Here, we implement a high-throughput workflow combining first-principles calculations, materials prediction, and magnetic properties calculations to discover novel magnetism and magnetic materials in borides. Successfully applying the workflow, we report three families of novel magnetic borides, including two families of borides exhibiting quantum magnetism. One is a family of dimerized quantum magnets among YCrB 4 -type borides, which provides a rare platform for studying the spin-gap quantum critical point. The other is a family of altermagnets among FeMo2B 2 -type borides, extending the magnetic orderings exhibited by borides beyond conventional ferromagnetism and antiferromagnetism. We also predict a family of magnetic laminate transition metal borides, known as the MAB phases, in the AlFe 2 B 2 -type family, which provide pure-phase or alloying candidates for studying magnetocaloric materials and the associated magnetic transitions. The workflow is expected to be used in further studies of novel magnetism and magnetic materials.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Evidence for Ga clusters in β-Ga 2 O 3 from Raman spectroscopy and density functional theory

Monoclinic gallium oxide (β-Ga 2 O 3 ) single crystals have a Raman mode at ∼250 cm −1 that is strongly correlated with free-electron density. Prior work attributed this peak to an electronic excitation of a shallow donor impurity band. However, heavily n-type thin films grown by metalorganic chemical vapor deposition or molecular beam epitaxy do not have the peak. In the present work, an alternate model is proposed: the 250 cm −1 Raman peak arises from Ga clusters, defined as two or more Ga atoms that form Ga–Ga bonds. Raman mapping reveals variations in the frequency that are consistent with a distribution of cluster sizes. The intensity of the peak decreases as the temperature is raised, attributed to melting of the Ga clusters. First-principles calculations indicate that the 250 cm −1 mode is due to Ga–Ga bond-stretching vibrations. As the Fermi energy is raised, the formation of Ga–Ga dimers becomes energetically favorable, explaining the correlation between n-type conductivity and the appearance of the Raman peak.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Monolayer 1T-LaN 2 : Dirac spin-gapless semiconductor of p -state and Chern insulator with a high Chern number

Two-dimensional transition-metal dinitrides have attracted considerable attention in recent years due to their rich magnetic properties. Here, we focus on rare-earth-metal elements and propose a monolayer of lanthanum dinitride with a 1T structural phase, 1T-LaN 2 . Using first-principles calculations, we systematically investigated the structure, stability, magnetism, and band structure of this material. It is a flexible and stable monolayer exhibiting a low lattice thermal conductivity, which is promising for future thermoelectric devices. Here, the monolayer shows the ferromagnetic ground state with a spin-polarized band structure. Two linear spin-polarized bands cross at the Fermi level forming a Dirac point, which is formed by the p atomic orbitals of the N atoms, indicating that monolayer 1T-LaN 2 is a Dirac spin-gapless semiconductor of p-state. When the spin–orbit coupling is taken into account, a large nontrivial indirect bandgap (86/354 meV) can be opened at the Dirac point, and three chiral edge states are obtained, corresponding to a high Chern number of C = 3, implying that monolayer 1T-LaN 2 is a Chern insulator. Importantly, this kind of band structure is expected to occur in more monolayers of rare-earth-metal dinitride with a 1T structural phase.

2D materials↗
DOE OSTI2022

Role of defects in ultra-high gain in fast planar tin gallium oxide UV-C photodetector by MBE

Here we report ultra-high responsivity of epitaxial (Sn x Ga 1–x ) 2 O 3 (TGO) Schottky UV-C photodetectors and experimentally identified the source of gain as deep-level defects, supported by first principles calculations. Epitaxial TGO films were grown by plasma-assisted molecular beam epitaxy on (–201) oriented n-type β-Ga 2 O 3 substrates. Fabricated vertical Schottky devices exhibited peak responsivities as high as 3.5 ×10 4 A/W at –5 V applied bias under 250 nm illumination with sharp cutoff shorter than 280 nm and fast rise/fall time in milliseconds order. Hyperspectral imaging cathodoluminescence (CL) spectra were examined to find the mid-bandgap defects, the source of this high gain. Irrespective of different tin mole fractions, the TGO epilayer exhibited extra CL peaks at the green band (~2.20 eV) not seen in β-Ga 2 O 3 along with enhancement of the blue emission-band (~2.64 eV) and suppression of the UV emission-band. Based on hybrid functional calculations of the optical emission expected for defects involving Sn in β-Ga 2 O 3 , V Ga –Sn complexes are proposed as potential defect origins of the observed green and blue emission-bands. Such complexes behave as acceptors that can efficiently trap photogenerated holes and are predicted to be predominantly responsible for the ultra-high photoconductive gain in the Sn-alloyed Ga 2 O 3 devices by means of thermionic emission and electron tunneling. Regenerating the V Ga –Sn defect complexes by optimizing the growth techniques, we have demonstrated a planar Schottky UV-C photodetector of the highest peak responsivity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Optical properties of vacancies in aluminum oxide (𝛼−Al 2 ⁢O 3 ) from first principles

We employ first-principles calculations based on hybrid density functional theory to investigate the structural and optical properties of the oxygen vacancy (𝑉 O ) and aluminum vacancy (𝑉 Al ) in 𝛼−Al 2 ⁢O 3 , the most stable (corundum) phase of alumina. Our calculations facilitate the identification of experimental excitation and luminescence spectra with specific electronic transitions at the vacancy sites. The absorption line shape for excitation of an electron at 𝑉$^0_O$ to the conduction-band minimum (CBM) is in excellent agreement with the 6.1 eV band detected by optical absorption spectroscopy, and we find that the 5.9 eV absorption band is generated by an internal electron transition at 𝑉$^0_O$. We confirm that the slowly decaying 3.0 eV emission band of the 𝐹 center is due to a triplet-singlet transition at 𝑉$^0_O$. Our calculations also reveal that the 4.8 eV/5.4 eV absorption and 3.8 eV emission bands assigned to the 𝐹 + center are generated by internal transitions at 𝑉$^{+1}_O$. The line shape for the excitation of an electron from 𝑉$^{+1}_O$ to the CBM agrees well with an absorption band at 6.4 eV, while the recombination of an electron with 𝑉$^{+2}_O$ also produces luminescence around 3.8 eV, but with a considerably broader line shape. For Al vacancies, we confirm that the most prevalent configuration in 𝛼−Al 2 ⁢O 3 is a split-vacancy configuration 𝑉 Al,s , and we calculate migration barriers for different directions in the corundum crystal. We predict absorption and emission spectra for the excitation and recombination of an electron localized at 𝑉$^{−3}_{Al}$ and 𝑉$^{−3}_{Al,s}$ sites with the CBM. The line shapes of the two 𝑉 Al configurations overlap and are considerably broader than the spectra corresponding to 𝑉 O .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Competing magnetic phases in Li-Fe-Ge kagome systems

Competing interlayer magnetic interactions in kagome magnets can lead to diverse magnetic phases, which enable promising topological or quantum material properties. In this work, the electronic structure and magnetic properties are studied using first-principles calculations for the Li⁢Fe6⁢G⁡e6, Li⁢Fe6⁢G⁡e4, and Li⁢Fe6⁢G⁡e5 compounds sharing the kagome F⁡e3⁢Ge layer motif but with different interlayer arrangements. For Li⁢Fe6⁢G⁡e4 and Li⁢Fe6⁢G⁡e5, the predicted magnetic ground states are collinear antiferromagnetic (AFM) states involving a mix of ferromagnetic (FM) and AFM interlayer orientations. In contrast, for Li⁢Fe6⁢G⁡e6, an incommensurate cycloidal spin spiral ground state is stabilized, although its energy is close to a collinear A-type AFM state. The analysis of magnetic RKKY exchange coupling confirms the results of calculations. The values of atomic magnetic moments are in good agreement with existing experimental estimations. Our experiments on Li⁢Fe6⁢G⁡e6 single crystals reveals AFM ordering at ∼540 K and transition to another magnetic phase with a small FM component, possibly with spin canting, below ∼270 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Enhanced mobility of cations and anions in the redox state: The polaronium mechanism

Hugely enhanced slow-ion diffusivity has been widely observed under extreme redox conditions and for unclear reasons. Aided by first-principles calculations on model systems of ZrO 2 , CeO 2 , BaTiO 3 and Li 4/3 Mn 2/3 O 2 , here we successfully explained the intriguing phenomenon by a polaronium mechanism. We found a polaronium, defined as a transitory complex of a polaronic electron or hole and a migrating counterion, becomes highly mobile when the counterion comes from a d 0 or f 0 cation (e.g., Zr 4+ , Ti 4+ and Ce 4 ) or a p 6 anion (e.g., O 2- ) in the host compound. Upon a redox reaction, the complex attains a d 1 /f 1 or p 5 configuration, which spontaneously forms because it is favored by an electron-phonon interaction (manifest as the Jahn-Teller effect in high symmetry systems) that enables local relaxation and lowers the system energy. Here our calculations found such interaction reaching its peak at the saddle point where the local environment is softest, so soft that it allows a reorientation of the anisotropic d/f/p orbital to minimize the electron repulsion locally. Since the complex may dissolve after a successful ion-migration event, the redox electron/hole can be recycled to form another free-radical-like polaronium elsewhere, thereby enhancing ion migration repeatedly. The proposed polaronium mechanism, which also operates in ceramics doped with mixed-valence cations, is most relevant under dynamic and extreme thermal/field/irradiation conditions where extra electrons/holes are abundantly generated by non-equilibrium redox reactions. For such operations, some with emerging applications, our diffusion-enhancing mechanism may provide new theoretical insight to help understand their material/microstructure stability and performance.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Highly Anisotropic Quasi‐Direct Organic Metal Halide Hybrids: A Platform for Polarization‐Sensitive Optoelectronics

Low-dimensional organic–inorganic metal halide hybrids (OMHHs) exhibit remarkable optical properties and enhanced environmental stability. We investigate a 1D OMHH with formula C 4 N 2 H 14 PbBr 4 , consisting of Pb–Br chains separated by organic cations, which shows a large Stokes shift (0.83 eV) and broadband emission. Through first-principles calculations and polarized Raman spectroscopy, we characterize the material's vibrational properties and identify the specific phonon modes that drive exciton self-trapping. Our novel GW/Bethe-Salpeter equation force formalism reveals that low-frequency phonons (∼ 100 cm −1 , primarily involving Pb–Br motions) couple strongly with excitons, with a remarkably high Huang-Rhys factor of 137 ± 4, and gives a pathway for ultrafast structural analysis during the absorption process. This phonon-exciton coupling mechanism explains the material's broadband emission and provides a pathway for controlling optical properties through vibrations and for tuning vibrations through optical excitations. The material also exhibits highly anisotropic optical properties and electronic transport, with bands that are dispersive along the Pb–Br chains but nearly flat in perpendicular directions, resulting in direction-dependent electrical conductivity that is calculated to be an order of magnitude higher along the chain direction and consistent with measurements. These combined properties make this system an excellent platform for polarization-sensitive optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

The critical role of synthesis conditions on small polaron carrier concentrations in hematite—A first-principles study

Achieving highly efficient energy conversion with transition metal oxides necessitates overcoming conductivity limitations due to the formation of small polarons. Detailed understanding of the interplay among intrinsic defects, dopants, and electron polarons can help devise strategies for achieving higher carrier concentrations, therefore improving carrier conductivity. This work employs first-principles calculations to reliably predict electron polaron concentrations in a prominent polaronic oxide, hematite (Fe 2 O 3 ), by resolving interactions between charged defects and electron polarons and keeping charge neutrality condition among all charged species. This work addresses that both V O and Fe i can be primary donors in undoped hematite depending on the synthesis conditions, such as synthesis temperature and oxygen partial pressure, despite the fact that V O owns an extremely high ionization energy compared to k B T. Furthermore, from calculations of a plethora of n-type dopants (group IV and V elements), we find that Ti, Ge, Sb, and Nb are able to raise electron polaron concentrations in hematite significantly without considering dopant clustering. However, the magnitude of electron polaron concentration increase would be smaller if the dopant has a high tendency of clustering, such as Ti. Here we reveal the critical role of synthesis conditions on tuning electron polaron concentrations of both undoped and doped hematite. Our theoretical analysis provides important insights and general design principles for engineering more conductive polaronic oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Computationally guided experimental validation of divacancy defect formation in 4H-SiC

Recent research into solid-state qubits for quantum information science has focused on optically addressable spin defects such as the negatively charged nitrogen-vacancy center in diamond and the neutrally charged divacancy (VV) in 4H-SiC as scalable quantum sensors and networking qubits. Within this context, direct investigations of the structural origin and defect formation dynamics of a sub-set of the VV center in 4H-SiC remain lacking. Here, we take a systematic experimental approach guided by predictions from first-principles simulations to gain a thorough mechanistic understanding of the VV defect formation and control in 4H-SiC. We study the effect of annealing time and temperature on VV formation in high-purity semi-insulating 4H-SiC samples following electron irradiation. Three different temperatures (1123, 1273, and 1473 K) and annealing duration (from 0.5 to 72 h) are chosen to explore VV formation in different regions. We find that samples annealed at 1273 K give the highest VV-related photoluminescence (PL) intensities, in agreement with the prediction from first-principles calculations. Furthermore, the logarithmic dependence of VV-related PL intensities on the annealing duration at 1273 K indicates that 1273 K provides sufficient thermal energy for silicon vacancy migration but not for VV migration. Together, these results suggest that efficient VV formation occurs above the V Si migration temperature and below the VV migration threshold.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Nanoscale modulation of flat bands via controllable charge density wave defects in 4⁢𝐻⁡𝑏−Ta⁢S 2

Electron correlation is a main driver of exotic quantum phases and their interplay. The 4H b-TaS 2 system, possessing an intrinsic heterostructure of 1T- and 1H -TaS 2 monolayers, offers a unique opportunity to control electron correlation by distorting the atomic lattice or tuning interlayer coupling. Here, we investigated intrinsically deformed charge-density waves (CDWs) in the 1T layer of 4H b-TaS 2 to elucidate and control their effects on flat bands using scanning tunneling microscopy and spectroscopy (STM/S) combined with first-principles calculations. We identified two types of CDW defects: Type 1 has structural distortion and locally suppressed flat bands, while Type 2 features an increased flat band filling factor of intact CDW structure. Density functional theory calculations indicate that a sulfur vacancy in the 1T layer distorts the CDW structure and gives rise to a Type 1, whereas a sulfur vacancy in the 1H layer reduces the interlayer charge transfer and leads to a Type 2. Furthermore, we demonstrated creating and erasing individual CDW defects via STM manipulation. Here, our findings provide a pathway to not only tune flat bands but also selectively manipulate the interaction between CDW, the atomic lattice, and interlayer coupling in strongly correlated systems with atomic precision.

Charge density waves↗
DOE OSTI2021

Mexican-hat potential energy surface in two-dimensional III 2 -VI 3 materials and the importance of entropy barrier in ultrafast reversible ferroelectric phase change

First-principles calculations reveal a Mexican-hat potential energy surface (PES) for two-dimensional (2D) In 2 Se 3 . This unique PES leads to a pseudo-centrosymmetric paraelectric β phase that resolves the current controversy between theory and experiment. We further show that while the α-to-β (ferroelectric-to-paraelectric) phase transition is fast and coherent, assisted by an in-plane shear phonon mode, a random distribution of the atoms in the trough of the PES acts as an entropy barrier against the reverse β-to-α transition. This will be the origin of the speed limitation of current In 2 Se 3 ferroelectric devices. Furthermore, if one orders the β phase (due to the formation of in-plane ferroelectric domains), the reverse transition can take place within tens of picoseconds in the presence of a perpendicular electric field. Finally, the Mexican-hat PES is a general feature for the entire family of 2D III 2 -VI 3 materials. Our finding offers a critical physical picture in controlling the ultrafast reversible phase transition in 2D In 2 Se 3 and other III 2 -VI 3 materials, which will benefit their practical industrial development for advanced ferroelectric devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Tuning the Magnetic Properties of CrI 3 Using Ni Thin Film Deposition for Applications in Spintronic Devices

Chromium triiodide (CrI 3 ), a van der Waals magnet, has recently been shown to host Ising ferromagnetism down to the monolayer limit. It is a potentially important material in 2D magnet-based applications, such as magnetic sensors and spintronic devices. Prior studies have revealed the coexistence of two different types of interlayer magnetic coupling, the antiferromagnetic (AFM) coupling near the surface and ferromagnetic (FM) coupling in the deep bulk layers, in pristine CrI 3 crystals below the Curie temperature. In this study, we used Ni thin film deposition to tune the surface magnetic states in bulk CrI 3 . A nanometer thickness (4 nm thick) of Ni was deposited on the surface of CrI 3 using electron-beam evaporation to form a Ni/CrI 3 heterostructure. The deposited Ni thin layer forms nanoclusters that completely cover the CrI 3 surface. Magnetic states of CrI 3 before and after Ni deposition are probed by ultralow-frequency magneto-Raman spectroscopy. Instead of seeing three spin wave branches in Raman scattering below 2 T, as in pristine bulk CrI 3 , in Ni/CrI 3 , we observe only a single spin wave branch, which is softened compared to that in pristine CrI 3 and displays a Zeeman shift under an out-of-plane magnetic field up to 7 T. This observation reveals that the AFM layers on pristine CrI 3 crystal surfaces transit into FM layers, so the entire CrI 3 crystal is in the FM state after Ni deposition. First-principles calculations show that Ni atoms tend to diffuse into the CrI 3 lattice, and the FM interlayer coupling has a much lower energy than AFM coupling in the presence of Ni atoms. As a result, our studies show that the magnetic state of CrI 3 can be modified through deposition of a thin metal layer on the surface, offering a route for controlling the spin degree of freedom in van der Waals magnets.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Topological Phases in Graphene Nanoribbons Tuned by Electric Fields

Graphene nanoribbons (GNRs) possess distinct symmetry-protected topological phases. Here we show, through first-principles calculations, that by applying an experimentally accessible transverse electric field, certain boron and nitrogen periodically codoped GNRs have tunable topological phases. The tunability arises from a field-induced band inversion due to an opposite response of the conduction- and valence-band states to the electric field. With a spatially varying applied field, segments of GNRs of distinct topological phases are created, resulting in a field-programmable array of topological junction states, each may be occupied with charge or spin. Our findings not only show that electric field may be used as an easy tuning knob for topological phases in quasi-one-dimensional systems, but also provide new design principles for future GNR-based quantum electronic devices through their topological characters.

1-dimensional systems↗