Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamic temperature”

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

Semimetallicity and electron-hole liquid in two-dimensional C and BN based compounds

Insulating-metallic transition mediated by substitutional atoms is predicted in a series of two-dimensional carbon-based structures. Introducing Si atoms in selected sites of tetrahexcarbon [Carbon 137 (2018) 266] according to rational chemical rules, metallicity by trivial band inversion without band gap opening is induced. Additional substitution of remaining C atoms by BN dimers introduces no changes in the metallic properties. A series of isomorphous two-dimensional materials with isoelectronic structures derived by exchanging group IV elements exhibiting various band gaps is obtained. Dynamical stability is verified with phonon analysis and beyond the harmonic approximation with molecular dynamics up to room temperature. The semi-metallic compounds have well-nested pockets of carriers and are good candidates for the formation of an excitonic insulator.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Shear-based deformation processing of age-hardened aluminum alloy for single-step sheet production

Shear-based deformation processing by hybrid cutting-extrusion and free machining are used to make continuous strip, of thickness up to one millimeter, from low-workability AA6013-T6 in a single deformation step. The intense shear can impose effective strains as large as 2 in the strip without pre-heating of the workpiece. The creation of strip in a single step is facilitated by three factors inherent to the cutting deformation zone: highly confined shear deformation, in situ plastic deformation-induced heating and high hydrostatic pressure. The hybrid cutting-extrusion, which employs a second die located across from the primary cutting tool to constrain the chip geometry, is found to produce strip with smooth surfaces (S a < 0.4 μm) that is similar to cold-rolled strip. The strips show an elongated grain microstructure that is inclined to the strip surfaces – a shear texture – that is quite different from rolled sheet. Furthermore, this shear texture (inclination) angle is determined by the deformation path. Through control of the deformation parameters such as strain and temperature, a range of microstructures and strengths could be achieved in the strip. When the cutting-based deformation was done at room temperature, without workpiece pre-heating, the starting T6 material was further strengthened by as much as 30% in a single step. In elevated-temperature cutting-extrusion, dynamic recrystallization was observed, resulting in a refined grain size in the strip. Implications for deformation processing of age-hardenable Al alloys into sheet form, and microstructure control therein, are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Network Uncertainty Quantification for Analysis of Multi-Component Systems

To impact physical mechanical system design decisions and realize the full promise of high-fidelity computational tools, simulation results must be integrated at the earliest stages of the design process. This is particularly challenging when dealing with uncertainty and optimizing for system-level performance metrics, as full-system models (often notoriously expensive and time-consuming to develop) are generally required to propagate uncertainties to system-level quantities of interest. Methods for propagating parameter and boundary condition uncertainty in networks of interconnected components hold promise for enabling design under uncertainty in real-world applications. These methods avoid the need for time consuming mesh generation of full-system geometries when changes are made to components or subassemblies. Additionally, they explicitly tie full-system model predictions to component/subassembly validation data which is valuable for qualification. These methods work by leveraging the fact that many engineered systems are inherently modular, being comprised of a hierarchy of components and subassemblies that are individually modified or replaced to define new system designs. By doing so, these methods enable rapid model development and the incorporation of uncertainty quantification earlier in the design process. The resulting formulation of the uncertainty propagation problem is iterative. We express the system model as a network of interconnected component models, which exchange solution information at component boundaries. We present a pair of approaches for propagating uncertainty in this type of decomposed system and provide implementations in the form of an open-source software library. We demonstrate these tools on a variety of applications and demonstrate the impact of problem-specific details on the performance and accuracy of the resulting UQ analysis. This work represents the most comprehensive investigation of these network uncertainty propagation methods to date.

42 ENGINEERING↗
DOE OSTI2020

Optically pumped spin polarization as a probe of many-body thermalization

Disorder and many body interactions are known to impact transport and thermalization in competing ways, with the dominance of one or the other giving rise to fundamentally different dynamical phases. Here we investigate the spin diffusion dynamics of 13 C in diamond, which we dynamically polarize at room temperature via optical spin pumping of engineered color centers. We focus on low-abundance, strongly hyperfine-coupled nuclei, whose role in the polarization transport we expose through the integrated impact of variable radio-frequency excitation on the observable bulk 13 C magnetic resonance signal. Unexpectedly, we find good thermal contact throughout the nuclear spin bath, virtually independent of the hyperfine coupling strength, which we attribute to effective carbon-carbon interactions mediated by the electronic spin ensemble. In particular, observations across the full range of hyperfine couplings indicate the nuclear spin diffusion constant takes values up to two orders of magnitude greater than that expected from homo-nuclear spin couplings.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Hydrodynamic instabilities and heat transfer characteristics in the duct flow of a fluid in the supercritical thermodynamic regime

The behavior of fluids at supercritical thermodynamic conditions is inherently complex due to large variations in thermodynamic and transport properties. Recent numerical and experimental investigations illustrate ongoing interest for these fluids, especially supercritical CO 2 and supercritical water, for a variety of applications. For example, supercritical water reactors (SCWR) operate in this extreme condition of high-pressure and temperature, resulting in highly dynamic flow fields and unexpected heat transfer regimes. The potential heat transfer benefits in this regime are directly associated with the extreme variations in thermodynamic and transport properties, which occur at, and above, the critical point. This work characterizes the hydrodynamic instabilities that arise for fluids at supercritical thermodynamic conditions when buoyancy forces are significant. Two specific configurations are considered, a natural convection cavity flow, and a mixed convection, heated, horizontal channel flow. Natural convection flow in a cavity is a classical configuration with expected behavior below the critical point. This configuration aids in characterizing the effect of the variable properties in the supercritical thermodynamic regime. Further, limited studies in the existing literature have been conducted for low-Reynolds and intermediate-Rayleigh numbers, mixed-convection channel flows for supercritical water, which is the focus of the channel flow configuration. To investigate the thermally driven hydrodynamic instabilities in this regime, a high-order fully-implicit numerical method is used. Such strong variations in thermophysical properties (in particular, density) are difficult to simulate and an altogether compressible framework is needed. Therefore, the compressible Navier-Stokes equations are solved without any additional assumptions. The fully implicit, high-order in space and time, reconstructed discontinuous Galerkin method as implemented within the multi-physics code called ALE3D (Arbitrary Lagrangian and Eulerian in 2D and 3D), developed at Lawrence Livermore National Laboratory (LLNL), is used. This fully implicit, L-stable method accurately captures the compressible nature of the ow in the limit of very low Mach number. It has been widely accepted that above the critical point, only one phase is observed. However, recent research has indicated the existence of the distinct gas-like and liquid-like regions separated by the Widom line, the locus of the maxima of the specific heat. Along the Widom line, density decreases 6-fold, viscosity drops by a factor of 2, while specific heat spikes by an order of magnitude. These variations, specifically in density and viscosity, produce a thick pseudo-interface and flow dynamics behavior akin to film boiling. A pseudo-film at the heated wall of the cavity and the horizontal channel is observed where buoyancy forces induce mixing through the specific configurations. Further the local Rayleigh and Richardson numbers provide maps of the flow field and the buoyancy forces driving the microscopic mixing. In the first chapter, I describe a background of supercritical fluid and the various applications. The second chapter focuses on the mathematical model and numerical method used for simulations, where a description of the equation of state for supercritical water is described. The third chapter focuses on the natural convection cavity with a heated bottom wall. In this cavity a gas-like and a liquid-like flow within the supercritical thermodynamic regime are observed. The fourth chapter focuses on a forced convection, horizontal channel, distinguishing between the gas-like, liquid-like, and mixed flow regimes. Mixed convection flow, with the addition of gravitational forces in the horizontal channel show the influence of variable properties on the hydrodynamic development, heat transfer, and rising instabilities. The last chapter of this research focuses on characterizing the unstable hydrodynamics through time-averaging processes and analysis of the movement of energy through the developing plumes.

42 ENGINEERING↗
DOE OSTI2025

Evaluation of the Energy, Hygrothermal, and Thermal Capacity Performance of Cross-Laminated Timber

Cross-laminated timber (CLT) construction is gaining momentum in the US because it offers multiple advantages over traditional construction methods. Benefits that have received the most attention focus on constructability, the environment, and protection (e.g., blast resistance), although CLT construction is likely to offer other benefits, as well. Still, these have not been studied at length because such evaluations are costly, requiring long-term assessments in an actual building and specialized technical knowledge. Among the possible benefits, CLT construction likely provides a higher-performing building envelope. Using CLT panels to enclose a building means fewer joints in the opaque envelope than what is required in traditional stick-framed construction. Fewer joints mean fewer locations where the air- and water-resistive barrier (WRB) could be compromised; thus, a CLT building enclosure may require less maintenance and have a longer lifespan than a traditionally built structure because of fewer air and water leaks. In addition, CLT’s thermal mass moderates indoor temperatures, allowing the heating, ventilation, and air conditioning (HVAC) system to operate more efficiently during peak hours, reducing operational energy consumption throughout the lifetime of the CLT building (Salonvaara et al., 2022). Furthermore, more stable indoor temperatures can increase occupant comfort. The CLT’s thermal mass can also reduce energy costs by adjusting to utility time-of-use pricing without affecting occupant comfort. The ability of CLT buildings to bridge periods without HVAC operation prepares them for future grid interaction and provides a certain level of resilience against power outages. Researchers have attempted to quantify these benefits; however, their work is based on simplified simulations with numerous assumptions. To correctly understand the benefits, an actual building must be monitored. Therefore, information needs to be gathered on indoor and outdoor temperatures, HVAC energy consumption, thermostat setpoints, temperatures, and thermal transport in CLT components to comprehend how these parameters are affected by the CLT’s thermal mass. These data are needed to reduce the number of assumptions and calibrate simulation models to optimize HVAC controls to minimize overall energy consumption, reduce energy use and higher fees during peak demand, and maintain occupant comfort. Additionally, the calibrated simulation model allows the optimization exercise to be repeated in various US climates. Potential benefits can be tailored to buildings in various locations, and decisions can be made on where CLT construction could be most advantageous. Furthermore, monitoring and simulation results are needed to evaluate the durability of the CLT structures in different climates. This project’s researchers gathered information to help understand and quantify the benefits of CLT buildings concerning operational energy, moderated indoor temperatures, and comfort; the dynamic operation to provide grid services; and resilience in times of power outage. Through the corroboration of simulation models with real-world measurements, this study paves the way for extrapolating findings to other climatic zones and building typologies, thereby broadening the understanding of CLT’s multifaceted benefits and reinforcing its position as a material of choice in sustainable construction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2026

Planetary Boundary-Layer Height (PBLHT) Value-Added Product: Remote-Sensing Retrievals

The planetary boundary layer (PBL) is fundamental to numerous atmospheric processes, including aerosol mixing and transport, cloud evolution, and precipitation formation. A critical parameter in these studies is the PBL height (PBLHT). This vertical depth is essential for characterizing PBL structures in numerical simulations and serves as a primary metric for estimating flux exchanges between the Earth’s surface and the atmosphere. Radiosonde (SONDE) observations provide high-vertical-resolution measurements of temperature and moisture profiles and are widely used to estimate PBLHT (Liu and Liang 2010, Seidel et al. 2010). The U.S. Department of Energy Atmospheric Radiation Measurement (ARM) User Facility’s PBLHT value-added product (VAP) for radiosonde measurements, known as PBLHTSONDE, applies three commonly used methods—the Heffter (1980) method, the Liu and Liang (2010) method, and the bulk Richardson number approach (Seibert et al. 2000)—to derive PBLHT. The PBLHTSONDE VAP operates routinely at ARM observatories and mobile facilities, with data available from the ARM Data Center shortly after sounding observations are collected (Sivaraman et al. 2013). However, radiosonde observations are limited by their low temporal resolution. Most stations launch soundings only twice daily, which constrains the ability to investigate and characterize the temporal evolution of the PBL using radiosonde data alone. The use of continuous remote-sensing observations provides high temporal resolution of PBLHT estimates. These observations include aerosol lidars (Dang et al. 2019, Su et al. 2020), Doppler lidar (DL; Tucker et al. 2009, Krishnamurthy et al. 2021), and water vapor and/or temperature lidars and radiometers (Turner et al. 2014). These observations provide valuable data on the PBL’s thermodynamic properties (e.g., water vapor and/or temperature lidars and radiometers), dynamic properties (e.g., DL), and distribution of tracer substances (e.g., aerosol lidars), all of which can be used to estimate PBLHT. ARM developed PBLHT estimates from the micropulse lidar (MPL; PBLHTMPL), Doppler lidar (PBLHTDL), and combined Raman lidar (RL)/atmospheric emitted radiance interferometer (AERI) thermodynamic profiles (PBLHTTHERMO). Each estimate captures different physical characteristics of the boundary layer—aerosol tracers, vertical velocity turbulence, and thermodynamic structure—and exhibits distinct strengths and limitations depending on the PBL regime and time of day. In addition, the ARM ceilometer (CEIL) provides three potential PBLHT candidates derived from the vendor's built-in algorithm. Building on these individual retrievals, ARM developed the PBLHTBEML VAP, which combines the four remote-sensing-based estimates with ancillary meteorological variables using the machine learning approach of Zhang et al. (2025) to produce a best-estimate PBLHT at 10-minute resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

First Use of a Longitudinally Polarized Target with CLAS12

Run Group C comprises eight experiments utilizing the CLAS12 detector system in Hall B at Jefferson Lab to study the multidimensional partonic structure of nucleons. The experiments scatter electrons from polarized protons and neutrons in samples of solid NH$_3$ and ND$_3$, dynamically polarized at a temperature of 1 K in a 5 T magnetic field. After a brief description of the target system, the current status of the target and preliminary results of its performance in the experiment are presented.

Keith, Christopher↗
DOE OSTI2023

Characterization of Interconnects

Combined-accelerated stress testing (C-AST) simultaneously combines stress factors of the natural environment (including UV radiation, temperature, humidity, electrical current, and external mechanical force) into a single test that requires fewer modules, fewer chambers, and makes it possible to discover weaknesses in new designs that are not known a-priori. C-AST reduces risk, accelerates time to market, and improves bankability by reducing costly overdesign using test levels not exceeding those seen in the natural environment. This project seeks to evaluate strengths and weaknesses of modern cell interconnect designs with combined-accelerated stress testing (C-AST) to screen multiple climates along with finite element and failure analysis to determine the potential for a 50-year life. Supporting this, we seek to develop methods of characterization to assess the degradation of interconnects, enable predictive rate models with material forensics and FEA and show C-AST's ability to find interconnect failures and benchmark it relative to other accelerated stress test methods like temperature cycling and cyclic dynamic mechanical loading. Examples of modern interconnect designs covered are Canadian Solar Hetero Technology ribbon enabling closer cell interconnections, SmartWire (eg. Meyer-Berger technology), and shingled cells with new interconnect materials.

C-AST↗
DOE OSTI2024

Thin, Highly Selective Polymer Membrane Separators for Advanced Liquid Alkaline Water Electrolysis

Our team will develop novel, dense, polymer electrolyte separators based on polybenzimidazole (PBI) derivatives that will enable tuning of water/KOH uptake in mechanically robust films for advanced liquid alkaline water electrolysis (LAWE). These material advances will enable thinner separators with increased selectivity (hydroxide conductivity/hydrogen permeability), thereby enabling higher operating currents, increased dynamic operating high-temperature stability, and catalyst-coated membranes (zero-gap design advantages) in direct alignment with stated lab call goals.

HYDROGEN↗
DOE OSTI2024

Thermally anisotropic composites for thermal management in building environments

An improved system for thermal management is provided. The system includes thermally anisotropic composites coupled with a thermal loop to re-direct, reduce, and shape heat flows through a building envelope, having the potential to (1) significantly reduce envelope-generated heating and cooling loads and (2) provide grid services such as decreasing peak loads and shaping energy use. In one embodiment, the thermal management system includes an anisotropic composite that consists of alternating layers of thermal insulation and thermally conductive materials that are immediately adjacent to each other, including polyisocyanurate foam boards and aluminum sheets. The thermal management system also includes a thermal loop along the long edge or the entire the perimeter of the anisotropic composite, the thermal loop having dynamically controlled or floating temperature that is maintained at lower than an outdoor ambient temperature (for cooling). An interior wall structure is inwardly adjacent to the anisotropic composite.

Shrestha, Som S.↗
DOE OSTI2020

Heterogeneous Polymer Dynamics Explored Using Static 1H NMR Spectra

NMR spectroscopy continues to provide important molecular level details of dynamics in different polymer materials, ranging from rubbers to highly crosslinked composites. It has been argued that thermoset polymers containing dynamic and chemical heterogeneities can be fully cured at temperatures well below the final glass transition temperature (Tg). In this paper, we described the use of static solid-state 1H NMR spectroscopy to measure the activation of different chain dynamics as a function of temperature. Near Tg, increasing polymer segmental chain fluctuations lead to dynamic averaging of the local homonuclear proton-proton (1H-1H) dipolar couplings, as reflected in the reduction of the NMR line shape second moment (M2) when motions are faster than the magnitude of the dipolar coupling. In general, for polymer systems, distributions in the dynamic correlation times are commonly expected. To help identify the limitations and pitfalls of M2 analyses, the impact of activation energy or, equivalently, correlation time distributions, on the analysis of 1H NMR M2 temperature variations is explored. It is shown by using normalized reference curves that the distributions in dynamic activation energies can be measured from the M2 temperature behavior. An example of the M2 analysis for a series of thermosetting polymers with systematically varied dynamic heterogeneity is presented and discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Fragile-to-strong transition in liquid As 2 S 3 under pressure: The effect of melt metallization

The well-known classification of glass-forming melts into fragile and strong liquids has several notable exceptions, including water, silica, and certain phase-change materials (PCMs). These exceptional fluid systems exhibit a fragile-to-strong transition (FST) upon cooling: a transformation from a high-temperature liquid with fast atomic dynamics, low viscosity, and low flow activation energy, to a viscous supercooled melt with high energy barriers near the glass transition temperature T g . This behavior is critically important for non-volatile memories, photonic tensor cores, reconfigurable metamaterials, and other devices, that use PCMs, enabling nanosecond-scale crystallization in the fragile regime and long data retention in the strong regime near or below T g . A significant structural transformation is expected between these two viscosity regimes, along with a semiconductor-metal (SC-M) transition upon heating, driven by high internal pressure and associated density increase. By applying high external pressure to the canonical low-conducting chalcogenide melt As 2 S 3 , we observed both the FST and the SC-M transition, occurring simultaneously within the same domain of the P, T−phase space. These findings suggest that the FST is not limited to a few exceptional liquids but is a common phenomenon, at least in systems that exhibit melt metallization within specific regions of their P, T−phase diagrams.

first-principles molecular dynamics↗
DOE OSTI2021

Structure of ice confined in silica nanopores

Observed anomalous thermodynamic properties of confined water such as deviations in the melting point and freezing point motivate the determination of the structure of confined water as a function of pore size and temperature. Here, we investigate the dynamic evolution of the structure of confined ice in SBA-15 porous materials with pore diameters of 4 nm, 6 nm, and 8 nm at temperatures ranging from 183 K to 300 K using in operando Wide-Angle X-Ray Scattering (WAXS) measurements, X-Ray Partial Distribution Function (PDF) measurements, and classical Molecular Dynamics (MD) simulations. Formation of hexagonal ice structures is noted in all the three pore sizes. In silica nanopores with diameters of 4 nm, cubic ice formation is noted in addition to hexagonal ice. Longer lasting hydrogen bonds and longer residence times of the water molecules in the first coordination shell contribute to observed crystalline organization of ice in confinement. Self-diffusion coefficients of confined liquid water, predicted from classical MD simulations, are four orders of magnitude higher compared to ice formed in confinement. These experimental and simulation results provide comprehensive insights underlying the organization of confined water and ice in silica nanopores and the underlying physico-chemical interactions that contribute to the observed structures.

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

Elevation distributed micro-climatology data in a coastal glaciated watershed

Air temperature, ground temperature and relative humidity data were collected in a longitudinal transect of the Nooksack watershed at varying elevations from 500 to 1800 m above sea level. Data were collected by anchoring sensors from trees above winter snow levels and shaded from direct solar radiation. Paired sensors were also buried 3 cm under ground near each air temperature sensor to determine snow absence or presence. Select sites included relative humidity sensors to indicate whether precipitation was occurring. Data were collected every 3-4 h from December 2015 to Sept 2018 (with ongoing collection). Code for analysis of daily mean, minimum, maximum, and temperature change with elevation (lapse rates) are available on Github (https://doi.org/10.5281/zenodo.3239539). The sensor download and intermediate data products are available on HydroShare at (http://www.hydroshare.org/resource/222e832d3df24dea9bae9bbeb6f4219d) with publicly accessible visualization available from the Nooksack Observatory at data.cuahsi.org. Hydrologic models are generally structured with a single annual average lapse rate parameter which assumes a linear temperature gradient with elevation. The daily data (2016-2018) is used as part of ongoing studies on the non-linear dynamics and temporal variability of temperature with elevation to improve assessments of watershed function and salmon habitat.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Addition of Chloroform in a Solvent-in-Salt Electrolyte: Outcomes in the Microscopic Dynamics in Bulk and Confinement

Solvent-in-salt electrolytes (SISEs) are a promising alternative to the electrolytes currently used in commercial devices. Despite the SISEs’ advantages, their utilization is not yet realized due to the poor mobility of their chemical species. Here, we explore this problem by adding chloroform to a SISE formed by acetonitrile and a Li-salt. First, we performed illustrative cycling experiments to highlight the potential of this approach. Then, we focused on the description of the microscopic dynamics of the electrolytes and exposed the relevant aspects to be considered for their optimal performance. While the conductivity at low temperatures may be enhanced by the addition of chloroform, only subtle changes occur at room temperature. As revealed by molecular dynamics simulations and quasielastic neutron scattering (QENS) experiments, this effect is related to the preservation of the structure expected for a highly concentrated solution and promotion of the formation of ionic aggregates. These outcomes occur despite the increase in the overall mobility of the chemical species. The dynamics of the electrolytes in porous carbon was also investigated using QENS. In these circumstances, low concentrations of chloroform lead to diffusivities of the molecular species higher than those observed for the bulk electrolytes. As chloroform’s concentration increases, no further changes in the diffusivities are observed. Nonetheless, chloroform is mostly immobilized on the carbon surfaces and this behavior may be intensified at compositions closer to the eutectic mixture.

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

Ab initio molecular dynamics with enhanced sampling in heterogeneous catalysis

Ab initio molecular dynamics simulations combined with enhanced sampling techniques are becoming widespread methods to investigate chemical phenomena in catalytic systems. These techniques automatically include finite temperature effects, anharmonicity, and collective dynamics in their robust description of enthalpic and entropic contributions, which can have significant impact on reaction free energy landscapes. This contrasts with standard ab initio static approaches that are based on assessing reaction free energies from various coarse-grained descriptions of the reaction potential energy surface. Enhanced sampling ab initio molecular dynamics opens the way to first principles simulations of systems of increasing complexity like solid/liquid catalytic interfaces. Here, we aim at guiding the reader through the basis of these techniques, summarizing their fundamental theoretical and practical aspects, and reviewing the relevant literature in the field. After a brief introduction to the problem, we will illustrate the advantage of using molecular simulations to include finite temperature effects, examine the most common ab initio techniques currently in use, describe their application to solid state heterogeneous catalysts, and finally critically review the most popular enhanced sampling techniques used in computational catalysis.

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

The Debye-Waller factor of disordered solids: The case for amorphous silicon

This paper presents a first-principles study of the Debye-Waller factor and the Debye temperature for amorphous silicon (α-Si) from lattice-dynamical calculations and direct molecular-dynamics simulations using density-functional theory (DFT). The effects of temperature and structural disorder on the intensity of the diffraction maxima and the vibrational mean-square displacement (MSD) of Si atoms are studied in the harmonic approximation, with particular emphasis on the bond-length disorder, the presence of coordination defects, and microvoids in α-Si networks. It has been observed that the MSDs associated with tetrahedrally bonded Si atoms are considerably lower than their dangling-bond counterparts—originating from isolated and vacancy-induced clustered defects—and those on the surface of microvoids, leading to an asymmetric non-Gaussian tail in the distribution of atomic displacements. An examination of the effect of anharmonicity on the MSD at high temperatures using direct ab initio molecular-dynamics simulations (without the harmonic approximation) suggests that the vibrational motion in α-Si is practically unaffected by anharmonic effects at temperatures below 400 K, as far as the present DFT calculations are concerned. The Debye temperature of α-Si is found to be in the range of 488–541 K from specific-heat and MSD calculations using first-principles lattice-dynamical calculations in the harmonic approximation, which matches closely with the experimental value of 487–528 K obtained from specific-heat measurements of α-Si at low temperatures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗