SEARCH · Engineering Papers

Results for “Cross points”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,099 records · Page 61

All sources NASA NTRS DOE OSTI

NASA NTRS1991

MHD simulations of magnetic reconnection in a skewed three-dimensional tail configuration

The dynamic evolution of a nonsymmetric magnetotail configuration initiated by the sudden occurrence of (anomalous) resistivity are examined using a three-dimensional resistive MHD simulation developed by Birn (1990) that includes a net cross-tail field and thus breaks the mirror symmetry around the neutral sheet. Results show that the field evolution is similar to that of a symmetric configuration studied by Birn and Hones (1981), pointing to the formation and ejection of a plasmoid. On the other hand, the topological structure of the magnetic field, defined by the field line connections, was remarkably different from the symmetric case. The plasmoid in this case became 'open', connected initially with the earth but gradually becoming connected with the interplanetary field. The openness of the plasmoid and its magnetic connection with interplanetary field lines suggest the possibility of a heat flux out of the plasmoid on interconnected flux tubes.

Birn, J.↗

NASA NTRS1992

Negative-ion formation in the explosives RDX, PETN, and TNT using the Reversal Electron Attachment Detection (READ) technique

In the search for high sensitivity and direct atmospheric sampling of trace species, techniques have been developed such as atmospheric-sampling, glow-discharge ionization (ASGDI), corona discharge, atmospheric pressure ionization (API), electron-capture detection (ECD), and negative-ion chemical ionization (NICI) that are capable of detecting parts-per-billion to parts-per-trillion concentrations of trace species. These techniques are based on positive- or negative-ion formation via charge-transfer to the target, or electron capture under multiple-collision conditions in a Maxwellian distribution of electron energies at the source temperature. One drawback of the high-pressure, corona- or glow-discharge devices is that they are susceptible to interferences either through indistinguishable product masses, or through undesired ion-molecule reactions. The ASGDI technique is relatively immune from such interferences, since at target concentrations of less than 1 ppm the majority of negative ions arises via electron capture rather than through ion-molecule chemistry. A drawback of the conventional ECD, and possibly of the ASGDI, is that they exhibit vanishingly small densities of electrons with energies in the range 0-10 millielectron volts (meV), as can be seen from a typical Maxwellian electron energy distribution function at T = 300 K. Slowing the electrons to these subthermal (less than 10 meV) energies is crucial, since the cross section for attachment of several large classes of molecules is known to increase to values larger than 10(exp -12) sq cm at near-zero electron energies. In the limit of zero energy these cross sections are predicted to diverge as epsilon(exp -1/2), where epsilon is the electron energy. In order to provide a better 'match' between the electron energy distribution function and attachment cross section, a new concept of attachment in an electrostatic mirror was developed. In this scheme, electrons are brought to a momentary halt by reversing their direction with electrostatic fields. At this turning point the electrons have zero or near-zero energy. A beam of target molecules is introduced, and the resultant negative ions extracted. This basic idea has been recently improved to allow for better reversal geometry, higher electron currents, lower backgrounds, and increased negative-ion extraction efficiency. We present herein application of the so-called reversal electron attachment detector (READ) to the study of negative-ion formation in the explosives molecules RDX, PETN, and TNT under single-collision conditions.

Chutijian, Ara↗

NASA NTRS1995

Numerical Study of Boundary Layer Interaction with Shocks: Method Improvement and Test Computation

The objective is the development of a high-order and high-resolution method for the direct numerical simulation of shock turbulent-boundary-layer interaction. Details concerning the spatial discretization of the convective terms can be found in Adams and Shariff (1995). The computer code based on this method as introduced in Adams (1994) was formulated in Cartesian coordinates and thus has been limited to simple rectangular domains. For more general two-dimensional geometries, as a compression corner, an extension to generalized coordinates is necessary. To keep the requirements or limitations for grid generation low, the extended formulation should allow for non-orthogonal grids. Still, for simplicity and cost efficiency, periodicity can be assumed in one cross-flow direction. For easy vectorization, the compact-ENO coupling algorithm as used in Adams (1994) treated whole planes normal to the derivative direction with the ENO scheme whenever at least one point of this plane satisfied the detection criterion. This is apparently too restrictive for more general geometries and more complex shock patterns. Here we introduce a localized compact-ENO coupling algorithm, which is efficient as long as the overall number of grid points treated by the ENO scheme is small compared to the total number of grid points. Validation and test computations with the final code are performed to assess the efficiency and suitability of the computer code for the problems of interest. We define a set of parameters where a direct numerical simulation of a turbulent boundary layer along a compression corner with reasonably fine resolution is affordable.

Adams, N. A.↗

NASA NTRS2004

Topography of the Flattest Surface on Earth: using ICESAT, GPS, and MISR to Measure Salt Surface Topography on Salar de Uyuni, Bolivia

Salt flats are aptly named: they are composed largely of salt, and are maintained as nearly equipotential surfaces via frequent flooding. The salar de Uyuni, on the Altiplano in southwestern Bolivia, is the largest salt flat on Earth, with an area of 9,800 sq km. Except for a few bedrock islands, it has less than 40 cm of relief. The upper-most salt unit averages 5 m thick and contains 50 cu km of nearly pure halite. It includes most of the salt that was in solution in paleolake Minchin, which attained a maximum area of 60,000 sq km and a maximum depth of 150 m, roughly 15 kyr ago. Despite approx. 10 m of differential isostatic rebound since deposition, the salar surface has been actively maintained as an extraordinarily flat and smooth surface by annual flooding during the rainy season. We have used the strong optical absorption properties of water in the visible band to map spatial variations in water depth during a time when the salar was flooded. As water depth increases, the initially pure white surface appears both darker and bluer. We utilized MISR images taken during the interval from April to November 2001. The red and infra-red bands (672 and 867 nm wavelength) were most useful since the water depth is small and the absorption at those wavelengths is quite strong. Nadir pointed MISR images have 275 m spatial resolution. To aid in our evaluation of water depth variations over the saiar surface, we utilized two sources of direct topographic measurements: several ICESAT altimetry tracks cross the area, and a 40x50 km GPS grid was surveyed to calibrate ICESAT. A difficulty in using these data types is that both give salt surface elevations relative to the ellipsoid, whereas the water surface will, in the absence of wind or tidal disturbances, follow an equipotential surface. Geoid height is not known to the required accuracy of a few cm in the central Andes. As a result, before comparing optical absorption from MISR to salt surface topography from GPS or ICESAT, we removed the longest wavelengths from both.

Comstock, Robert L.↗

NASA NTRS2005

Analysis of Bare-Tether Systems as a Thruster for MXER Studies

The concept of the electrodynamic propulsion has a number of attractive features and has been widely discussed for different applications. A number of system designs have been proposed and compared during the last ten years. In spite of this, the choice of the proper design, for a specific mission, is far from evident. Such characteristics of tether performance as system acceleration, efficiency, etc. should be calculated and compared. The code that calculates the current for bare and partly insulated tethers with circular (wire) and rectangle (tape) cross-sections is presented. It takes into account the corrections to the OML current due to the tether cross-section geometry and the magnetic field produced by the tether current. There are two options in this code: for current calculation with the prescribed energy supply and with the prescribed end-point potential. This permits to calculate the parameters characterizing tether performance. Results for the current calculated for tethers with different designs for the currently proposed Momentum exchange Electrodynamic Reboost (MXER) Tether System are presented.

Khazanov, G. V.↗

NASA NTRS2005

Analysis of Bare-Tether Systems as a Thruster for MXER Studies

The concept of the electrodynamic propulsion has a number of attractive features and has been widely discussed for different applications. A number of system designs have been proposed and compared during the last ten years. In spite of this, the choice of the proper design, for a specific mission, is far from evident. Such characteristics of tether performance as system acceleration, efficiency, etc. should be calculated and compared. The code that calculates the current for bare and partly insulated tethers with circular (wire) and rectangle (tape) cross-sections is presented. It takes into account the corrections to the OML current due to the tether cross-section geometry and the magnetic field produced by the tether current. There are two options in this code: for current calculation with the prescribed energy supply and with the prescribed end-point potential. This permits us to calculate the parameters characterizing tether performance. Results for the current calculated for tethers with different designs for the currently proposed Momentum exchange Electrodynamic Reboost (MXER) Tether System are presented.

Khazanov, G. V.↗

NASA NTRS2009

Further Development of Ko Displacement Theory for Deformed Shape Predictions of Nonuniform Aerospace Structures

The Ko displacement theory previously formulated for deformed shape predictions of nonuniform beam structures is further developed mathematically. The further-developed displacement equations are expressed explicitly in terms of geometrical parameters of the beam and bending strains at equally spaced strain-sensing stations along the multiplexed fiber-optic sensor line installed on the bottom surface of the beam. The bending strain data can then be input into the displacement equations for calculations of local slopes, deflections, and cross-sectional twist angles for generating the overall deformed shapes of the nonuniform beam. The further-developed displacement theory can also be applied to the deformed shape predictions of nonuniform two-point supported beams, nonuniform panels, nonuniform aircraft wings and fuselages, and so forth. The high degree of accuracy of the further-developed displacement theory for nonuniform beams is validated by finite-element analysis of various nonuniform beam structures. Such structures include tapered tubular beams, depth-tapered unswept and swept wing boxes, width-tapered wing boxes, and double-tapered wing boxes, all under combined bending and torsional loads. The Ko displacement theory, combined with the fiber-optic strain-sensing system, provide a powerful tool for in-flight deformed shape monitoring of unmanned aerospace vehicles by ground-based pilots to maintain safe flights.

Ko, William L.↗

NASA NTRS2010

X-Ray Diffraction Wafer Mapping Method for Rhombohedral Super-Hetero-Epitaxy

A new X-ray diffraction (XRD) method is provided to acquire XY mapping of the distribution of single crystals, poly-crystals, and twin defects across an entire wafer of rhombohedral super-hetero-epitaxial semiconductor material. In one embodiment, the method is performed with a point or line X-ray source with an X-ray incidence angle approximating a normal angle close to 90 deg, and in which the beam mask is preferably replaced with a crossed slit. While the wafer moves in the X and Y direction, a narrowly defined X-ray source illuminates the sample and the diffracted X-ray beam is monitored by the detector at a predefined angle. Preferably, the untilted, asymmetric scans are of {440} peaks, for twin defect characterization.

Park, Yoonjoon↗

NASA NTRS2000

Deformation and Flexibility Equations for Idealized ARIS Umbilicals, Under Planar End-Loading Conditions

The International Space Station (ISS) relies on the Active Rack Isolation System (ARIS) as the central component of an integrated, station-wide strategy to isolate microgravity space-science experiments. ARIS uses electromechanical actuators to isolate an International Standard Payload Rack (ISPR) from disturbances due to the motion of the ISS. Disturbances to microgravity experiments on ARIS-isolated racks are primarily transmitted via the ARIS power and vacuum umbilicals. Recent experimental tests indicate that these umbilicals resonate at frequencies outside the ARIS controller's bandwidth. at levels of potential concern for certain microgravity experiments. Reduction in the umbilical resonant frequencies could help to address this issue. This paper develops equations for the in-plane deflections and flexibilities of an idealized umbilical (thin, flexible, cantilever beam) under end-point, in-plane loading (inclined-force and moment). The effect of gravity is neglected due to the on:orbit application. The analysis assumes an initially straight. cantilevered umbilical with uniform cross-section. which undergoes large deflections with no plastic deformation, such that the umbilical terminus remains in a single quadrant and the umbilical slope changes monotonically. The analysis is applicable to the ARIS power and vacuum umbilicals. under the indicated assumptions.

Hampton, R. David↗

DOE OSTI2026

Dark Energy Survey Year 6 Results: Cosmological Constraints from Galaxy Clustering and Weak Lensing

We present cosmology results combining galaxy clustering and weak gravitational lensing measured in the full six years (Y6) of observations by the Dark Energy Survey (DES) covering $\sim$5000 deg$^2$. We perform a large-scale structure analysis using three two-point correlation functions (3$\times$2pt): (i) cosmic shear from 140 million source galaxy shapes, (ii) galaxy clustering of 9 million lens galaxy positions, and (iii) galaxy-galaxy lensing from their cross-correlation. We model the data in flat $Λ$CDM and $w$CDM cosmologies. The combined analysis yields $S_8\equiv σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.789^{+0.012}_{-0.012}$ and matter density $Ω_{\rm m} = 0.333^{+0.023}_{-0.028}$ in $Λ$CDM (68% CL), where $σ_8$ is the clustering amplitude. These constraints show a (full-space) parameter difference of 1.8$σ$ from a combination of cosmic microwave background (CMB) primary anisotropy datasets from Planck 2018, ACT-DR6, and SPT-3G DR1. Projected only into $S_8$ the difference is $2.6σ$. In $w$CDM the Y6 3$\times$2pt results yield $S_8 = 0.782^{+0.021}_{-0.020}$, $Ω_{\rm m} = 0.325^{+0.032}_{-0.035}$, and dark energy equation-of-state parameter $w = -1.12^{+0.26}_{-0.20}$. For the first time, we combine all DES dark-energy probes: 3$\times$2pt, SNe Ia, BAO and Clusters. In $Λ$CDM this combination yields a $2.8σ$ parameter difference from the CMB. When combining DES 3$\times$2pt with other low-redshift datasets (DESI DR2 BAO, DES SNe Ia, SPT clusters), we find a 2.3$σ$ parameter difference with CMB. A joint fit of Y6 3$\times$2pt, CMB, and those low-redshift datasets produces the tightest $Λ$CDM constraints to date: $S_8 = 0.806^{+0.006}_{-0.007}$, $Ω_{\rm m} = 0.302^{+0.003}_{-0.003}$, $h = 0.683^{+0.003}_{-0.002}$, and $\sum m_ν< 0.14$ eV (95% CL). In $w$CDM, this combination yields $w = -0.981^{+0.021}_{-0.022}$.

Abbott, T. M.C. [Cerro-Tololo InterAmerican Obs.]↗

DOE OSTI2025

Deployment of neural-network-based neutron microscopic cross sections in the Griffin reactor physics application

The capability to utilize neural networks to predict macroscopic and microscopic cross section parametric spaces has been developed for the Griffin reactor physics application. The LibTorch interface enables Griffin's MOOSE-based materials to interact with LibTorch-trained models, allowing for the evaluation of complex macroscopic or microscopic cross section spaces, which are then used to evaluate the neutronic properties of the Griffin finite element model. This study benchmarks traditional ISOXML-formatted tabulation libraries against neural network-based models for 279 nuclides on 20,160 grid points for zero-dimensional and two-dimensional reactor models. Benchmark metrics include the fundamental mode eigenvalue, fission and absorption rates, and various temperature coefficients of reactivity (isothermal, fuel, and moderator). From the perspective of storage space, the complete set of LibTorch models uses 11 MB on disk, compared to the 10 GB for the ISOXML multigroup library that covers the same grid space. For the two-dimensional performance case considered in Griffin, the Torch model uses 97% less RAM than the reference ISOXML dataset while runtime increases by a factor of 3 when using the LibTorch model compared to the ISOXML dataset with multi-linear interpolation. The LibTorch model consistently yields errors within 0.01% for most analyzed quantities except for the temperature coefficients of reactivity where the maximum discrepancies are up to 0.3 $\frac{pcm}{K}$. Due to the neural network attempting to best predict quantities with no regard for a positive or negative bias for any given quantity, predictions may experience random fluctuations, resulting in both positive and negative errors. Future work will entail both depletion and coupled transient analysis to determine the predictive capabilities of Griffin with neural network-based cross sections.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Simulating the Multi-epoch Direct Detection Technique to Isolate the Thermal Emission of the Non-transiting Hot Jupiter HD187123b

We report the 6.5σ detection of water from the hot Jupiter HD187123b with a Keplerian orbital velocity K {sub p} of 53 ± 13 km s{sup −1}. This high-confidence detection is made using a multi-epoch, high-resolution, cross-correlation technique, and corresponds to a planetary mass of 1.4{sub −0.3}{sup +0.5} M {sub J} and an orbital inclination of 21° ± 5°. The technique works by treating the planet/star system as a spectroscopic binary and obtaining high signal-to-noise, high-resolution observations at multiple points across the planet’s orbit to constrain the system’s binary dynamical motion. All together, seven epochs of Keck/NIRSPEC L-band observations were obtained, with five before the instrument upgrade and two after. Using high-resolution SCARLET planetary and PHOENIX stellar spectral models, we were able to drastically increase the confidence of the detection by running simulations that could reproduce, and thus remove, the nonrandom structured noise in the final likelihood space well. The ability to predict multi-epoch results will be extremely useful for furthering the technique. Here, we use these simulations to compare three different approaches to combining the cross correlations of high-resolution spectra and find that the Zucker log(L) approach is least affected by unwanted planet/star correlation for our HD187123 data set. Furthermore, we find that the same total signal-to-noise ratio (S/N) spread across an orbit in many, lower S/N epochs rather than fewer, higher S/N epochs could provide a more efficient detection. This work provides a necessary validation of multi-epoch simulations, which can be used to guide future observations and will be key to studying the atmospheres of farther separated, non-transiting exoplanets.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

The Importance of Avoided Crossings in Understanding High Valley Degeneracy in Half‐Heusler Thermoelectric Semiconductors

Abstract Half‐Heusler (hH) compounds are promising candidates for inexpensive, low‐toxicity thermoelectric materials. It is well known that engineering electronic bands with high valley degeneracy is an effective approach for enhancing the performance of thermoelectric materials, and there are several routes for achieving high valley degeneracy in hH systems. For instance, there are multiple locations in the first Brillouin zone where the valence band maximum can be found (at the Γ‐, L‐, or W‐point), and there are two competing low‐lying conduction bands at the X‐point, where the conduction band minimum is located. By converging the multiple valence band and conduction band extrema, the valley degeneracy, and hence, performance of these materials can be improved. Here, group theoretical and tight‐binding approaches, in addition to first‐principles density functional theory calculations, are used to study the chemical origins of various band extrema in both the n‐type and p‐type compounds, with particular focus on ZrNiSn and NbFeSb. Specifically, the importance of avoided crossings is explained. The results of this work can be used to better understand and develop design strategies for engineering better performing hH thermoelectrics.

Brod, Madison K.↗

DOE OSTI2023

Four-dimensional microstructurally small fatigue crack growth in cyclically loaded nickel superalloy specimen

Here, we nondestructively image a low-solvus high refractory (LSHR) nickel base specimen microstructure during cyclic loading to characterize crack propagation during the microstructurally small crack (MSC) growth regime. Micro-computed tomography (uCT) is used to characterize the material density and identify crack interfaces, while high-energy diffraction microscopy (HEDM) provides snapshots of the microstructural orientation field changes. From these measurements, data are computationally reconstructed and coregistered to describe the specimen’s configuration before load, before crack nucleation, and throughout the crack propagation/loading process. An MSC forms at a FIB-cut notch, which was previously placed to concentrate stress at the center of the specimen gauge, and propagates along a facet parallel to the local {1 1 1} crystallographic plane. Analysis of available crystallographic slip directions and measured point-to-point disorientation axes indicate probable single and coplanar slip activity along FCC octahedral planes within this region. Similar analyses of other crack propagation facets reveal generally slower growth along directions favoring zig-zag cross and codirectional cross slip, which are often disfavored at room temperature.

36 MATERIALS SCIENCE↗

DOE OSTI2025

A generalizable machine learning-assisted fast Fourier transform algorithm to simulate the large strain phenomena in polycrystalline materials

Machine learning methods have shown initial promise in constitutive modeling for single crystals or homogenized polycrystals, delivering notable computational efficiency. However, existing machine learning-based constitutive models often lack generalizability, limiting their application across diverse boundary value problems. This study introduces a thermodynamics-informed artificial neural network model to accelerate rate-tangent crystal plasticity fast Fourier transform simulations for cross-scale deformation behaviors of polycrystals under complex loading. Our model integrates microstructural variability and local interactions effectively. To address local effects in each grain, we employ K-means clustering to group Gauss points within the microstructure into clusters assumed to be in similar mechanical states. This approach, based on self-clustering analysis, extends model scope from macroscopic stress response to the granular level, capturing mechanical responses and orientation evolution across grains. This reduces the number of nonlinear problems to solve, with cluster responses propagated throughout each group. The thermodynamics-based artificial neural network-extracted features are further processed using local material state clusters to account for history-dependent deformation and evolving microstructures. Additionally, representative volume element simulations with rate-tangent crystal plasticity fast Fourier transform provide reliable datasets for model training. The proposed model demonstrates high efficiency, accuracy, self-consistency, and enhanced generalizability in predicting strain–stress responses and orientation evolution at both individual grain and aggregate scales under complex loading conditions, such as biaxial tension and arbitrary loading scenarios.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Biogeochemistry of upland to wetland soils, sediments, and surface waters across Mid-Atlantic and Great Lakes coastal interfaces

Transferable and mechanistic understanding of cross-scale interactions is necessary to predict how coastal systems respond to global change. Cohesive datasets across geographically distributed sites can be used to examine how transferable a mechanistic understanding of coastal ecosystem control points is. To address the above research objectives, data were collected by the EXploration of Coastal Hydrobiogeochemistry Across a Network of Gradients and Experiments (EXCHANGE) Consortium – a regionally distributed network of researchers that collaborated on experimental design, methodology, collection, analysis, and publication. The EXCHANGE Consortium collected samples from 52 coastal terrestrial-aquatic interfaces (TAIs) during Fall of 2021. At each TAI, samples collected include soils from across a transverse elevation gradient (i.e., coastal upland forest, transitional forest, and wetland soils), surface waters, and nearshore sediments across research sites in the Great Lakes and Mid-Atlantic regions (Chesapeake and Delaware Bays) of the continental USA. The first campaign measures surface water quality parameters, bulk geochemical parameters on water, soil, and sediment samples, and physicochemical parameters of sediment and soil.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

First-order and gradual phase transitions of ethane confined in MCM-41

The first-order phase transition of ethane confined in MCM-41, i.e., capillary condensation, has been measured using an isochoric cooling procedure by differential scanning calorimetry (DSC) under conditions ranging from 206 K and 1.1 bar up to the pore critical point (PCP). The PCP has also been determined using the three-line method developed earlier based on the vanishing heat of phase transition. As in the bulk phase, no first-order phase transition can occur above the critical point, which also implies that vapor can transform into liquid gradually by following a path around the critical point through the supercritical region. For the first time, the gradual phase transition is demonstrated with ethane in MCM-41, which is achieved through a multistep process with paths proceeding around the PCP without crossing the capillary-condensation curve. Finally, the occurrence of the gradual phase transition in nanopores, thus the confined supercriticality, is confirmed while our consistent DSC measurements are also well demonstrated.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2021

Chiral Topological Surface States on a Finite Square Photonic Crystal Bounded by Air

Chiral, topologically protected, photonic surface states can be found at the boundary between gyrotropic photonic crystals where a changing magnetic field induces different topology across the interface. Typically, photonic crystals with either a suitable band structure on both sides of the interface to provide a band gap and evanescent decay of the surface states away from the interface, or an outer layer with engineered material properties is required. In this paper, we show the emergence of topological, unidirectional surface states at the termination of finite gyrotropic photonic crystals with a simple square lattice and C 4 rotational symmetry bounded by a vacuum, eliminating the need for an outside layer to enable chiral surface modes. Here, we start from an infinite, time-reversal-symmetry-breaking photonic crystal with a band gap associated with bands with nonzero Chern numbers, different from all-zero Chern numbers in air. We then modify the photonic crystal to move this band gap below the light line, while maintaining the Chern-number discontinuities. Band-structure calculations for a supercell approximating a photonic crystal finite in the direction normal to the surface demonstrate the existence, dispersion, and chirality of the surface mode. Extensive direct scattering calculations for a point source and spatial Fourier analysis further reveal a unidirectional free-space topological surface state, which propagates counterclockwise around the surface of a finite photonic crystal, providing a nearly foolproof way to cross-check the surface-mode band structure unaffected by backscattering from local defects. Additionally, scattering simulations allow an independent characterization of the state dispersion and unveil the robustness of the topological plasmonic mode propagation around the 90° bends of the structure, being due to only radiation leakage. In contrast to buried topological surface states, the observed surface modes at the photonic crystal–air interface have the advantage of being accessible to the outside world, allowing one to take advantage of the defect-tolerant backscattering-free surface modes to engineer emission from photonic crystal surfaces into arbitrary free-space beam shapes and directions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗