Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “zigzag”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

The Crystal Structure of BaZn 2 Se 2 (OH) 2 Featuring Brownmillerite‐Type Layers

Abstract A bimetallic hydroxychalcogenide, BaZn 2 Se 2 (OH) 2 , was synthesized through hydrothermal pouch methods. The single crystal X‐ray diffraction and electron diffraction indicates that the phase crystallizes in the orthorhombic space group Pnma and is composed of anionic layers [ZnSe 3/3 (OH) 1/1 ] − that are separated and charged balanced by Ba 2+ cations. The [ZnSe 3/3 (OH) 1/1 ] – layer comprises two unique Zn sites, which form interpenetrating zigzag chains with an in‐plane dipole moment and adopts a brownmillerite‐type structural motif. The adjacent layers contain tetrahedrally coordinated Zn chains of opposite handedness related by an inversion center, which cancel the microscopic dipoles to minimize the macroscopic electric polarization. The adoption of a brownmillerite structural motif in BaZn 2 Se 2 (OH) 2 can be rationalized by the distinct charge difference between Se 2− and OH − anions, which creates a sufficient dipole moment in the ZnSe 3 (OH) tetrahedra to allow the occurrence of twisted chains. FTIR spectroscopy confirms the existence of OH − anions and DFT calculations indicate that BaZn 2 Se 2 (OH) 2 is a semiconductor with a direct band gap. This work expands the chemistry of the brownmillerite family from traditional homoanionic oxides to multianion hydroxychalcogenides, offering a new opportunity to explore tunable structural complexity for better design of functional materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Strain-modulated intercalated phases of Pb monolayer with dual periodicity in SiC(0001)-graphene interface

Intercalation of metal atoms at the SiC(0001)-graphene (Gr) interface can provide confined 2D metal layers with interesting electronic properties. The intercalated Pb monolayer (ML) has shown the coexistence of the Gr(10 x 10)-moiré and a stripe phase, which still lacks understanding. Using density functional theory calculation and thermal annealing with ab initio molecular dynamics as motivated by experiment, we have studied the formation energy of Gr/Pb/SiC(0001) for different Pb coverages. Near the coverage of a Pb(111)-like ML mimicking the (10 x 10)-moiré, we find a slightly more stable stripe structure, where one half of the structure has compressive strain with Pb occupying the Si-top sites and the other half has tensile strain with Pb off the Si-top sites. This stripe structure along the Gr zigzag direction has a periodicity of 2.3 nm across the [1$\overline{2}$10] direction agreeing with the previous observations using scanning tunneling microscopy. Analysis with electron density difference and density of states show the tensile region has a more metallic character than the compressive region, while both are dominated by charge transfer from Pb ML to SiC(0001). As a result, the small energy difference between the stripe and Pb(111)-like structures means the two phases are almost degenerate and can coexist, which explains the experimental observations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Mechanistic study of direct coupling of CO 2 and C 2 H 4 over atomically dispersed metal at graphene edges

Direct coupling of CO 2 and ethylene (hereinafter DCCE) to acrylic acid is valuable for valorizing CO 2 to manufacture acrylate-derived products. However, previous studies in DCCE have been limited on molecular catalysts with challenges in improving catalytic performance. In this work, we employed density functional theory calculations and ab initio molecular dynamics simulations to investigate the heterogeneous catalysis of DCCE over atomically dispersed metal centers at nitrogen-doped zigzag edge of graphene. Based on competitive adsorption and structural stability, Mo, Cr, V, Ru, and Ni active sites are chosen to explore the reaction kinetics. Here, we find that the activation barriers are determined by the charge redistribution at transition states, which explains the trend of activity for the C-C coupling and the hydrogen transfer, two key steps in DCCE. Furthermore, we show that the intramolecular hydrogen transfer (rate-limiting step) is hindered due to the lack of local coordinate at the active sites. We thus propose to use co-adsorbed water as a “proton-exchanger” following a water-assisted route, and show that the activation barriers are reduced over all metal centers. Particularly, water promotes the hydrogen transfer over metals with strong CO 2 -ethylene co-activation and facile C-C coupling kinetics, which could be considered promising for DCCE. In both mechanisms, the stability of metallactone intermediate can be used to predict the catalytic activity. It is anticipated that the insights from this work can provide guidelines for mimicking well-defined multifunctional active sites in molecular catalysts to design heterogeneous catalysts for such C-C coupling, which advances catalytic utilization of CO 2 .

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Angle-resolved polarized Raman study of layered b-As x P x-1 alloys: Identification of As-P vibrational modes

In this study, the polarization-resolved Raman spectra of b-As x P x-1 flakes with varying arsenic concentration (x = 0, 0.4, and 0.8) are systematically investigated. A clear polarization dependence is observed in the Raman intensity of all the P-P, As-P, and As-As modes and confirmed by the first-principles calculations. The observed angle-dependence can be related to the structure of the Raman tensors of each peak. In particular, the Raman intensities corresponding to As-P vibration modes are studied in detail. The out of plane As-P mode ($A$$^{1}_{g}$) is clearly identified, but the butterfly-shaped behavior observed in the polar plot around 350 cm -1 is interpreted as the combination of the two in-plane vibrations, one along the armchair and the other along the zigzag direction. This work allows unambiguous assignment of the $A$$^{1}_{g}$, B 2g , and $A$$^{2}_{g}$ modes of the As-P vibrational modes and the importance of the consideration of the superposition of the in-plane B 2g and $A$$^{2}_{g}$ modes. Such study can be utilized to further characterize the anisotropic structural properties and analyze the local distribution of the As and P in b-As x P x-1 alloys which will play key roles in electronic and optical properties. Furthermore, this methodology can be further extended to alloyed two-dimensional material systems with in-plane anisotropy in general.

2D materials↗
DOE OSTI2021

Giant anisotropic in-plane thermal conduction induced by Anomalous phonons in pentagonal PdSe 2

In two-dimensional materials, different atomic stacking induces anisotropic atomic interactions and phonon dispersions, leading to the anisotropy of in-plane thermal transport. In this work, we report an exceptional case in layered pentagonal PdSe 2 , where the bonding, force constants, and lattice constants are nearly-equal along the in-plane crystallographic axis directions, while the thermal conductivity is surprisingly much greater along b-axis than along a-axis with a ratio up to 1.8. Such strong anisotropy is not only unexpected in in-plane uniform structured materials, but also 2 comparable to the record high in-plane anisotropic thermal conductivity in the nonuniform structured material reported to date (the ratio is ~2.0 in TiS 3 ). By combining the inelastic X-ray scattering measurement and the first-principles calculations, we attribute such high anisotropy to the low-energy phonons along a-axis, particularly their lower group velocities and “avoided-crossing” behavior. The different buckling structures between a- (zigzag-type) and b-axis (flat-type) are mainly responsible for such unique phonon dynamics properties of PdSe 2 . This finding helps to discover materials with high anisotropic in-plane thermal conductivity in uniform structures and reveals new physics of anisotropy of in-plane thermal conduction. Due to the unique features in structure and thermal transport properties, PdSe 2 may serve as a new platform for designing novel devices to route heat flow precisely at the nanoscale.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Unveiling field-transverse spin anisotropy in the spin liquid candidate α-RuCl3 via spin Hall magnetoresistance

α-RuCl3 has emerged as a possible candidate for a quantum spin liquid (QSL) that promises exotic quasiparticles relevant for fault-tolerant quantum computation. Here, we report spin-sensitive transport measurements using a proximal spin Hall metal, platinum (Pt), to probe magnetic moments in the insulator α-RuCl3. We observe spin Hall magnetoresistance (SMR), where both the transverse and longitudinal resistivities exhibit angular oscillations between the in-plane magnetic field and the current, driven by the interplay between the spin Hall effect in Pt and local magnetic moments in α-RuCl3. These oscillations occur from 1.5 to 18 T, covering the zigzag antiferromagnetic, putative QSL, and supposedly partially field-polarized phases. The phase of the SMR oscillations suggests that the local moments, whether static or fluctuating, develop spin anisotropy with a quantization axis in-plane and largely transverse to the magnetic field across all fields from 1.5 to 18 T. Temperature dependence indicates that the spin anisotropy operates at an energy scale similar to that of the reported QSL signatures in α-RuCl3.

Idzuchi, Hiroshi [Tohoku University, Japan]↗
DOE OSTI2026

Notes on remanent magnetization measurements in superconductors and hard ferromagnets

Data on zero applied field measurements of remanent magnetization and magnetic relaxation in a BCS superconductor LuNi 2 B 2 C and several hard ferromagnets are presented and compared. Apparent similarities and differences, in particular in Thermoremanent Magnetization (TRM) - like, Isothermal Remanent Magnetization (IRM) - like, and remanent magnetization measurements with zigzag temperature sweep measurements are outlined. It is discussed how these results could be relevant for the magnetization measurements in diamond anvil cells.

Ferromagnets↗
DOE OSTI2021

In-situ TEM observation of shear induced microstructure evolution in Cu-Nb alloy

Phase boundaries in multiphase alloys govern defect interaction and chemical intermixing across different phases during plastic deformation. Dynamic interaction of defects with phase boundaries in multiphase alloys, especially for immiscible alloys, has been topic of significant research interest in recent years. In this work, we describe a novel approach for carrying out in-situ TEM shear deformation to directly observe interfacial microstructural evolution of a Cu-Nb alloy. A unique double shear specimen geometry is microfabricated by a focused ion beam technique to apply shear deformation upon push loading inside the TEM. From the real-time observation, we discover that the phase boundary with a zigzag morphology effectively blocks stacking faults nucleated in a Cu grain from slipping into a Nb grain. Meanwhile, the Cu phase bears the most plastic deformation through slip or twinning mechanisms. This work sheds light on understanding the shear deformation and the behavior of phase boundaries in multiphase alloys during shear deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Adsorption, charge transfer and a coverage-driven transition of alkali metals on rutile TiO 2 (110)

The interaction of alkali metals with metal oxide surfaces is central to tuning surface reactivity in heterogeneous catalysis and photocatalysis. Here we present a comprehensive DFT+U study of the adsorption of alkali metals (Li, Na, K, Rb, Cs) on the (110) surface of rutile TiO 2 . At low coverage (θ = 1/8), all alkali metals bind preferentially to bridging oxygen sites with adsorption energies in the range −4.06 to −3.33 eV, transferring nearly one full electron (0.92–0.99 |e|) to the substrate and inducing Ti 4+ → Ti 3+ reduction. The excess charge localizes preferentially at subsurface Ti sites in the form of small polarons. Diffusion barriers indicate facile motion along bridging-oxygen rows, whereas inter-row hopping is strongly hindered. Coverage effects were examined systematically for potassium: adsorption energy and charge transfer per K atom decrease monotonically with increasing θ. Strikingly, a sharp energy discontinuity occurs between θ = 4/8 and θ = 5/8 (ΔE ≈ 1 eV per atom), which we identify as a coverage-driven structural transition arising from steric packing constraints and enhanced K–K electrostatic repulsion once every (1×1) surface cell is occupied. This structural transition perfectly correlates with a dramatic drop in the work function down to an ultra-low minimum of 0.84 eV at θ=5/8, followed by a metallization- driven recovery at higher coverages. Ab initio molecular dynamics simulations confirm zigzag K arrangements at moderate coverage (θ = 1/3), while at high coverage (θ = 2/3) short-range K–K correlations emerge without long-range order. These results provide atomistic insight into the structure–activity relationships underlying alkali promotion effects on oxide-supported catalysts.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Ba 4 RuMn 2 O 10 : A Noncentrosymmetric Polar Crystal Structure with Disordered Trimers

Phase-pure polycrystalline Ba 4 RuMn 2 O 10 was prepared and determined to adopt the noncentrosymmetric polar crystal structure (space group Cmc2 1 ) based on results of second harmonic generation, convergent beam electron diffraction, and Rietveld refinements using powder neutron diffraction data. The crystal structure features zigzag chains of corner-shared trimers, which contain three distorted face-sharing octahedra. The three metal sites in the trimers are occupied by disordered Ru/Mn with three different ratios: Ru1:Mn1 = 0.202(8):0.798(8), Ru2:Mn2 = 0.27(1):0.73(1), and Ru3:Mn3 = 0.40(1):0.60(1), successfully lowering the symmetry and inducing the polar crystal structure from the centrosymmetric parent compounds Ba 4 T 3 O 10 (T = Mn, Ru; space group Cmca). The valence state of Ru/Mn is confirmed to be +4 according to X-ray absorption near-edge spectroscopy. Ba 4 RuMn 2 O 10 is a narrow bandgap (~0.6 eV) semiconductor exhibiting spin-glass behavior with strong magnetic frustration and antiferromagnetic interactions.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Anomalous Ferromagnetic Behavior in Orthorhombic Li 3 Co 2 SbO 6

Monoclinic Li 3 Co 2 SbO 6 has been proposed as a Kitaev spin liquid candidate and investigated intensively, whereas the properties of its polymorph, the orthorhombic phase, are less known. In this work, we report the magnetic properties of orthorhombic Li 3 Co 2 SbO 6 as revealed by dc and ac magnetic susceptibility, muon spin relaxation (μSR), and neutron diffraction measurements. Successive magnetic transitions at 115, 89, and 71 K were observed in the low-field dc susceptibility measurements. The transitions below T N (115 K) are suppressed at higher applied fields. However, zero-field ac susceptibility measurements reveal distinct frequency-independent transitions at about 114, 107, 97, 79, and 71 K. A long-range magnetic ordered state was confirmed by specific heat, μSR, and neutron diffraction measurements, all indicating a single transition at about 115 K. The discrepancy between different measurements is attributed to possible stacking faults and/or local disorders of the ferromagnetic zigzag chains, resulting in ferromagnetic boundaries within the overall antiferromagnetic matrix.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

The Role of Methylaryl Radicals in the Growth of Polycyclic Aromatic Hydrocarbons: The Formation of Five-Membered Rings

The regions of the C 13 H 11 potential energy surface (PES) related to the unimolecular isomerization and decomposition of the 1-methylbiphenylyl radical and accessed by the 1-/2-methylnaphthyl + C 2 H 2 reactions have been explored by ab initio G3(MP2,CC)//B3LYP/6-311G(d,p) calculations. The kinetics of these reactions relevant to the growth of polycyclic aromatic hydrocarbons (PAH) under high-temperature conditions in circumstellar envelopes and in combustion flames has been studied employing the RRKM-Master Equation approach. The unimolecular reaction of 1-methylbiphenylyl proceeding via a five-membered ring closure followed by H elimination is predicted to be very fast, on a submicrosecond scale above 1000 K and to result in the formation of an embedded five-membered ring in the 9H-fluorene product. The 1-/2-methylnaphthyl + C 2 H 2 reaction mechanism involves acetylene addition to the radical on the methylene group followed by a six- or five-membered ring closure and aromatization via an H atom loss. Despite of the complexity of the C 13 H 11 PES, these straightforward pathways are dominant in the high-temperature regime (above ~1000 K), with the prevailing products being phenalene, with a significant contribution of 1H-cyclopenta(a)naphthalene, for 1-methylnaphthyl + C 2 H 2 , and 1H-cyclopenta(b)naphthalene and 3H-cyclopenta(a)naphthalene, for 2-methylnaphthyl + C 2 H 2 . The methylnaphthyl reactions with acetylene represent a clean source of the three-ring PAHs, but they are relatively slow owing to the high entrance barriers of ~10 kcal/mol, with the rate constants of about an order of magnitude lower as compared to those for naphthyl + allene and σ-aryl + C 2 H 2 . Here, the 1-methylnaphthyl + C 2 H 2 and 2-methylnaphthyl + C 2 H 2 reactions represent prototypes for PAH growth by an extra six- and five-membered ring on a zigzag edge or a corner of PAH and the generated modified Arrhenius expressions are recommended for kinetic modeling of PAH expansion by the mechanism of acetylene addition to methylaryl radicals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Helical Organic and Inorganic Polymers

Despite being a staple of synthetic plastics and biomolecules, helical polymers are scarcely studied with Gaussian-basis-set ab initio electron-correlated methods on an equal footing with molecules. This article introduces an ab initio second-order many-body Green’s-function [MBGF(2)] method with nondiagonal, frequency-dependent Dyson self-energy for infinite helical polymers using screw-axis-symmetry-adapted Gaussian-spherical-harmonics basis functions. Together with the Gaussian-basis-set density-functional theory for energies, analytical atomic forces, translational-period force, and helical-angle force, it can compute correlated energy, quasiparticle energy bands, structures, and vibrational frequencies of an infinite helical polymer, which smoothly converge at the corresponding oligomer results. These methods can handle incommensurable structures, which have an infinite translational period and are hard to characterize by any other method, just as efficiently as commensurable structures. We apply them to polyethylene (2/1 helix), polyacetylene (Peierls’ system), and polytetrafluoroethylene (13/6 helix) to establish the quantitative accuracy of MBGF(2)/cc-pVDZ in simulating their (angle-resolved) ultraviolet photoelectron spectra, and of B3LYP/cc-pVDZ or 6-31G** in reproducing their structures, infrared and Raman band positions, phonon dispersions, and (coherent and incoherent) inelastic neutron scattering spectra. We then predict the same properties for infinitely catenated chains of nitrogen or oxygen and discuss their possible metastable existence under ambient conditions. They include planar zigzag polyazene (N2)x (Peierls’ system), 11/3-helical isotactic polyazane (NH)x, 9/4-helical isotactic polyfluoroazane (NF)x, and 7/2-helical polyoxane (O)x as potential high-energy-density materials.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Nature of Oxygen Adsorption on Defective Carbonaceous Materials

Here, plane-wave density functional theory has been used to study oxygen adsorption on graphene, graphite, and (12,0) zigzag single-walled carbon nanotubes with and without Stone–Wales (SW) and single-vacancy (SV) defects to understand the role of defects on carbonaceous material reactivity. Atomic oxygen adsorption leads to the formation of an epoxide on defect-free graphene and graphite and an ether on the exterior wall of carbon nanotubes and SW-defected materials. O 2 chemisorption is endothermic on defect-free graphene and graphite and slightly exothermic on defect-free nanotubes. O 2 chemisorption energies are predicted to be –1.1 to –1.4 eV on an SW defect and –6.0 to –8.0 eV on an SV defect. An SW defect lowers the energy barriers by 0.90 and 0.50 eV for O 2 chemisorption on graphene and nanotubes, respectively. The formation of a C–O–O–C group is important for O 2 dissociation on defect-free and SW-defected materials. The energy barrier is less than 0.30 eV on an SV defect. The more reactive SW defect toward O adsorption on graphene is mostly due to the strained defective carbon atoms being able to donate more electrons to an O to form an ether. The larger 2s character in the hybrid orbitals in an ether than in an epoxide makes the ether C–O bond stronger. Stronger C–O binding on an SW-defective carbon nanotube than on a defect-free nanotube is in part due to more flexibility of the defect to release the epoxide ring strain to form an ether.

25 ENERGY STORAGE↗
DOE OSTI2022

Understanding Lithium Local Environments in LiMn 0.5 Ni 0.5 O 2 Cathodes: A DFT-Supported 6 Li Solid-State NMR Study

LiMn 0.5 Ni 0.5 O 2 has been plagued by the structural failure driven by the irreversible phase change upon subsequent cycling and the structural degradation caused by the seemingly unavoidable Li/Ni exchange. This begs for a better understanding of the correlation between lithium local structures, synthesis conditions, and overall electrochemical properties of materials. In this study, we use density functional theory (DFT) to assist experimental 6 Li solid-state magic angle spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy to understand the nature and change in lithium local structures with extended annealing times (168 vs 15 h) synthesized at 900 °C. To correlate the NMR peak changes after different synthesis conditions, we calculated the NMR spectra of representative low-energy model configurations of LiMn 0.5 Ni 0.5 O 2 composition–zigzag, row, and flower structures–as well as the influence of Ni/Mn ordering and Li/Ni exchange on each spectrum using a supercell that is large enough to allow us to consider the effect of different relative configurations of adjacent transition metal (TM) layers. The analysis based on a combination of the calculated NMR and the experimental spectra suggests that most configurations contributing to the experimental spectra are in fact composed of a blend of configurations higher in energy than the classical well-ordered ground-state structures. Further, the extended annealing of LiMn 0.5 Ni 0.5 O 2 slightly enhances Ni/Mn and Li/Mn orderings while reducing the Li/Ni mixing ratio. This study shows how DFT calculations are crucial in providing a better understanding of the experimental characterization data and therefore local environments and domain structures in oxide materials, such as Li-ion battery cathodes.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Strain-Induced Shifts in Defective Graphite Phonon Modes Predicted by Density Functional Theory

Carbon fiber composites have gained attention as a structural material because of their high strength-to-weight ratio, and understanding the effect of defects on reactivity and mechanical properties is important for the longevity and safety of the composite. Although it is known that strain causes the underlying graphitic vibrational modes to redshift, it is not clear how strain may alter reactivity and defect-induced vibrational changes. To investigate the strain-induced phonon changes of defective carbon fiber composites, density functional theory calculations of graphite are used, including intercalated hydrogen and fluorine defects. By comparing changes in the bond lengths, formation energies, and phonon density of states for uniaxially and biaxially strained graphite, strain was found to generally make defect formation more favorable and the specific behavior changes are dependent on the strain direction and defect identity. Specifically, intercalated fluorine phonons are more sensitive to strain than hydrogen intercalation phonons, and strain applied along the zigzag direction alters the calculated properties more than strain along the armchair direction. Finally, these results highlight the importance of understanding the microstructural effect of deviations from the ideal material because small changes in strain or defect type can significantly alter the behavior of the carbon fiber composite core.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Theoretical Study of High Harmonic Generation in Monolayer NbSe 2

High harmonic generation (HHG) is a powerful probe of electron dynamics on attosecond to femtosecond time scales and has been successfully used to detect electronic and structural changes in solid-state quantum materials, including transition-metal dichalcogenides (TMDs). Among TMDs, bulk NbSe 2 exhibits charge density wave (CDW) order below 33 K and becomes superconducting below 7.3 K. Monolayer NbSe 2 also has superconducting and CDW behavior and is therefore interesting as a material whose different structural and electronic properties could be probed via HHG. Here, we predict the HHG response of the pristine 2H and CDW phases of monolayer NbSe 2 using real-time time-dependent density functional theory under the application of a simulated laser pulse excitation. We find that due to the lack of inversion symmetry in both monolayer phases, it is possible to excite even harmonics and that the even harmonics appear as the transverse components of the current response under excitations polarized along the zigzag direction of the monolayer, while odd harmonics arise from the longitudinal current response in all excitation directions. This suggests that the even and odd harmonic responses can be controlled via the polarization of the probing field, opening an avenue for potentially useful applications in optoelectronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Chemical Vapor Transport Route toward Black Phosphorus Nanobelts and Nanoribbons

Chemical vapor transport (CVT) method is widely used for bulk black phosphorus (BP) fabrication. In this work, we demonstrate that CVT provides a route for the fabrication of BP nanoribbons and nanobelts. This method consists of a two-step procedure, including initial BP column growth using the CVT technique, followed by ultrasonic treatment and centrifugation. The obtained nanostructures preserve BP column dimensions, forming ultralong ribbon-like structures with the length to the width aspect ratio of up to 500. Computational modeling of the growth mechanism of a BP flake is also presented in support of the observed columnar growth. Calculation of the average energy of the molecule in the asymmetric flakes shows that the growth of the structure in the zigzag direction is more energetically favorable than in the armchair direction.

36 MATERIALS SCIENCE↗