Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “parameterization”

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 919 records · Page 51
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

An ensemble constrained variational analysis of atmospheric forcing data and its application to evaluate clouds in CAM5: Ensemble 3DCVA and Its Application

Large-scale atmospheric forcing data can greatly impact the simulations of atmospheric process models including Large Eddy Simulations (LES), Cloud Resolving Models (CRMs) and Single-Column Models (SCMs), and impact the development of physical parameterizations in global climate models. This study describes the development of an ensemble variationally constrained objective analysis of atmospheric large-scale forcing data and its application to evaluate the cloud biases in the Community Atmospheric Model (CAM5). Sensitivities of the variational objective analysis to background data, error covariance matrix and constraint variables are described and used to quantify the uncertainties in the large-scale forcing data. Application of the ensemble forcing in the CAM5 SCM during March 2000 intensive operational period (IOP) at the Southern Great Plains (SGP) of the Atmospheric Radiation Measurement (ARM) program shows systematic biases in the model simulations that cannot be explained by the uncertainty of large-scale forcing data, which points to the deficiencies of physical parameterizations. The SCM is shown to overestimate high clouds and underestimate low clouds. These biases are found to also exist in the global simulation of CAM5 when it is compared with satellite data.

58 GEOSCIENCES↗
DOE OSTI2016

Potential sea salt aerosol sources from frost flowers in the pan-Arctic region

In order to better represent observed wintertime aerosol concentrations at Barrow, Alaska, we implemented an observationally-based parameterization for estimating sea salt production from frost flowers in the Community Earth System Model (CESM). In this work, we evaluate the potential influence of this sea salt source on the pan-Arctic (60ºN-90ºN) climate. Results show that frost flower salt emissions substantially increase the modeled surface sea salt aerosol concentration in the winter months when new sea ice and frost flowers are present. The parameterization reproduces both the magnitude and seasonal variation of the observed submicron sea salt aerosol concentration at surface in Barrow during winter much better than the standard CESM simulation without a frost-flower salt particle source. Adding these frost flower salt particle emissions increases aerosol optical depth by 10% and results in a small cooling at surface. Here, the increase in salt particle mass concentrations of a factor of 8 provides nearly two times the cloud condensation nuclei concentration, as well as 10% increases in cloud droplet number and 40% increases in liquid water content near coastal regions adjacent to continents. These cloud changes reduce longwave cloud forcing by 3% and cause a small surface warming, increasing the downward longwave flux at the surface by 2 W m -2 in the pan-Arctic under the present-day climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Energy considerations in the Community Atmosphere Model (CAM)

An error in the energy formulation in the Community Atmosphere Model (CAM) is identified and corrected. Ten year AMIP simulations are compared using the correct and incorrect energy formulations. Statistics of selected primary variables all indicate physically insignificant differences between the simulations, comparable to differences with simulations initialized with rounding sized perturbations. The two simulations are so similar mainly because of an inconsistency in the application of the incorrect energy formulation in the original CAM. CAM used the erroneous energy form to determine the states passed between the parameterizations, but used a form related to the correct formulation for the state passed from the parameterizations to the dynamical core. If the incorrect form is also used to determine the state passed to the dynamical core the simulations are significantly different. In addition, CAM uses the incorrect form for the global energy fixer, but that seems to be less important. The difference of the magnitude of the fixers using the correct and incorrect energy definitions is very small.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Complex functionality with minimal computation: Promise and pitfalls of reduced‐tracer ocean biogeochemistry models

Abstract Earth System Models increasingly include ocean biogeochemistry models in order to predict changes in ocean carbon storage, hypoxia, and biological productivity under climate change. However, state‐of‐the‐art ocean biogeochemical models include many advected tracers, that significantly increase the computational resources required, forcing a trade‐off with spatial resolution. Here, we compare a state‐of‐the art model with 30 prognostic tracers (TOPAZ) with two reduced‐tracer models, one with 6 tracers (BLING), and the other with 3 tracers (miniBLING). The reduced‐tracer models employ parameterized, implicit biological functions, which nonetheless capture many of the most important processes resolved by TOPAZ. All three are embedded in the same coupled climate model. Despite the large difference in tracer number, the absence of tracers for living organic matter is shown to have a minimal impact on the transport of nutrient elements, and the three models produce similar mean annual preindustrial distributions of macronutrients, oxygen, and carbon. Significant differences do exist among the models, in particular the seasonal cycle of biomass and export production, but it does not appear that these are necessary consequences of the reduced tracer number. With increasing CO 2 , changes in dissolved oxygen and anthropogenic carbon uptake are very similar across the different models. Thus, while the reduced‐tracer models do not explicitly resolve the diversity and internal dynamics of marine ecosystems, we demonstrate that such models are applicable to a broad suite of major biogeochemical concerns, including anthropogenic change. These results are very promising for the further development and application of reduced‐tracer biogeochemical models that incorporate “sub‐ecosystem‐scale” parameterizations.

58 GEOSCIENCES↗
DOE OSTI2015

An improved hindcast approach for evaluation and diagnosis of physical processes in global climate models

Here, we present an improved procedure of generating initial conditions (ICs) for climate model hindcast experiments with specified sea surface temperature and sea ice. The motivation is to minimize errors in the ICs and lead to a better evaluation of atmospheric parameterizations' performance in the hindcast mode. We apply state variables (horizontal velocities, temperature and specific humidity) from the operational analysis/reanalysis for the atmospheric initial states. Without a data assimilation system, we apply a two-step process to obtain other necessary variables to initialize both the atmospheric (e.g., aerosols and clouds) and land models (e.g., soil moisture). First, we nudge only the model horizontal velocities towards operational analysis/reanalysis values, given a 6-hour relaxation time scale, to obtain all necessary variables. Compared to the original strategy in which horizontal velocities, temperature and specific humidity are nudged, the revised approach produces a better representation of initial aerosols and cloud fields which are more consistent and closer to observations and model's preferred climatology. Second, we obtain land ICs from an offline land model simulation forced with observed precipitation, winds, and surface fluxes. This approach produces more realistic soil moisture in the land ICs. With this refined procedure, the simulated precipitation, clouds, radiation, and surface air temperature over land are improved in the Day 2 mean hindcasts. Following this procedure, we propose a “Core” integration suite which provides an easily repeatable test allowing model developers to rapidly assess the impacts of various parameterization changes on the fidelity of modelled cloud-associated processes relative to observations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

On the relationships among cloud cover, mixed‐phase partitioning, and planetary albedo in GCMs

Abstract In this study, it is shown that CMIP5 global climate models (GCMs) that convert supercooled water to ice at relatively warm temperatures tend to have a greater mean‐state cloud fraction and more negative cloud feedback in the middle and high latitude Southern Hemisphere. We investigate possible reasons for these relationships by analyzing the mixed‐phase parameterizations in 26 GCMs. The atmospheric temperature where ice and liquid are equally prevalent (T5050) is used to characterize the mixed‐phase parameterization in each GCM. Liquid clouds have a higher albedo than ice clouds, so, all else being equal, models with more supercooled liquid water would also have a higher planetary albedo. The lower cloud fraction in these models compensates the higher cloud reflectivity and results in clouds that reflect shortwave radiation (SW) in reasonable agreement with observations, but gives clouds that are too bright and too few. The temperature at which supercooled liquid can remain unfrozen is strongly anti‐correlated with cloud fraction in the climate mean state across the model ensemble, but we know of no robust physical mechanism to explain this behavior, especially because this anti‐correlation extends through the subtropics. A set of perturbed physics simulations with the Community Atmospheric Model Version 4 (CAM4) shows that, if its temperature‐dependent phase partitioning is varied and the critical relative humidity for cloud formation in each model run is also tuned to bring reflected SW into agreement with observations, then cloud fraction increases and liquid water path (LWP) decreases with T5050, as in the CMIP5 ensemble.

58 GEOSCIENCES↗
DOE OSTI2016

Rain evaporation and moist patches in tropical boundary layers

Moist patches are areas in the subcloud layer characterized by a positive water vapor anomaly compared to the environment and are considered important in triggering new convective cells. A correct understanding of the origin of the water vapor in these patches is, thus, essential to improving existing convective parameterizations. Recent studies have addressed this problem and have shown that contrary to what was previously thought, the main source of water vapor in moist patches are surface latent heat fluxes, instead of rain evaporation. This paper offers a different perspective to the topic, focusing on the origin only of the water vapor that makes moist patches anomalously moist when compared to the environment. Finally, it is found that near the surface, rain evaporation contributes half as much as latent heat fluxes, implying that a parameterization of the thermodynamic forcing should be more sensitive to environmental variables, like relative humidity, than recently suggested.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Model representations of aerosol layers transported from North America over the Atlantic Ocean during the Two-Column Aerosol Project

The ability of the Weather Research and Forecasting model with chemistry (WRF-Chem) version 3.7 and the Community Atmosphere Model version 5.3 (CAM5) in simulating profiles of aerosol properties is quantified using extensive in situ and remote sensing measurements from the Two Column Aerosol Project (TCAP) conducted during July of 2012. TCAP was supported by the U.S. Department of Energy’s Atmospheric Radiation Measurement program and was designed to obtain observations within two atmospheric columns; one fixed over Cape Cod, Massachusetts and the other several hundred kilometers over the ocean. The performance is quantified using most of the available aircraft and surface measurements during July, and two days are examined in more detail to identify the processes responsible for the observed aerosol layers. The higher resolution WRF-Chem model produced more aerosol mass in the free troposphere than the coarser resolution CAM5 model so that the fraction of aerosol optical thickness above the residual layer from WRF-Chem was more consistent with lidar measurements. We found that the free troposphere layers are likely due to mean vertical motions associated with synoptic-scale convergence that lifts aerosols from the boundary layer. The vertical displacement and the time period associated with upward transport in the troposphere depend on the strength of the synoptic system and whether relatively high boundary layer aerosol concentrations are present where convergence occurs. While a parameterization of subgrid scale convective clouds applied in WRF-Chem modulated the concentrations of aerosols aloft, it did not significantly change the overall altitude and depth of the layers.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2017

Dynamical conditions of ice supersaturation and ice nucleation in convective systems: A comparative analysis between in situ aircraft observations and WRF simulations: Ice Supersaturation in Convective System

Occurrence frequency and dynamical conditions of ice supersaturation (ISS, where relative humidity with respect to ice (RHi) > 100%) are examined in the upper troposphere around convective activity. Comparisons are conducted between in situ airborne observations and the Weather Research and Forecasting model simulations using four double-moment microphysical schemes at temperatures ≤-40°C. All four schemes capture both clear-sky and in-cloud ISS conditions. However, the clear-sky (in-cloud) ISS conditions are completely (significantly) limited to the RHi thresholds of the Cooper parameterization. In all of the simulations, ISS occurrence frequencies are higher by ~3–4 orders of magnitude at higher updraft speeds (>1 m s -1 ) than those at the lower updraft speeds when ice water content (IWC) > 0.01 g m -3 , while observations show smaller differences up to ~1–2 orders of magnitude. The simulated ISS also occurs less frequently at weaker updrafts and downdrafts than observed. These results indicate that the simulations have a greater dependence on stronger updrafts to maintain/generate ISS at higher IWC. At lower IWC (≤0.01 g m -3 ), simulations unexpectedly show lower ISS frequencies at stronger updrafts. Overall, the Thompson aerosol-aware scheme has the closest magnitudes and frequencies of ISS >20% to the observations, and the modified Morrison has the closest correlations between ISS frequencies and vertical velocity at higher IWC and number density. The Cooper parameterization often generates excessive ice crystals and therefore suppresses the frequency and magnitude of ISS, indicating that it should be initiated at higher ISS (e.g., ≥25%).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Using large eddy simulations to reveal the size, strength, and phase of updraft and downdraft cores of an Arctic mixed-phase stratocumulus cloud

Three-dimensional large eddy simulations (LES) are used to analyze a springtime Arctic mixed-phase stratocumulus observed on 26 April 2008 during the Indirect and Semi-Direct Aerosol Campaign. Two subgrid-scale turbulence parameterizations are compared. The first scheme is a 1.5-order turbulent kinetic energy (1.5-TKE) parameterization that has been previously applied to boundary layer cloud simulations. The second scheme, Cloud Layers Unified By Binormals (CLUBB), provides higher-order turbulent closure with scale awareness. The simulations, in comparisons with observations, show that both schemes produce the liquid profiles within measurement variability but underpredict ice water mass and overpredict ice number concentration. The simulation using CLUBB underpredicted liquid water path more than the simulation using the 1.5-TKE scheme, so the turbulent length scale and horizontal grid box size were increased to increase liquid water path and reduce dissipative energy. The LES simulations show this stratocumulus cloud to maintain a closed cellular structure, similar to observations. The updraft and downdraft cores self-organize into a larger meso-γ-scale convective pattern with the 1.5-TKE scheme, but the cores remain more isotropic with the CLUBB scheme. Additionally, the cores are often composed of liquid and ice instead of exclusively containing one or the other. Furthermore, these results provide insight into traditionally unresolved and unmeasurable aspects of an Arctic mixed-phase cloud. From analysis, this cloud's updraft and downdraft cores appear smaller than other closed-cell stratocumulus such as midlatitude stratocumulus and Arctic autumnal mixed-phase stratocumulus due to the weaker downdrafts and lower precipitation rates.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Broadband optical properties of biomass-burning aerosol and identification of brown carbon chromophores

The radiative effects of biomass burning aerosols on regional and global scale is substantial. Accurate modeling of the radiative effects of smoke aerosols require wavelength-dependent measurements and parameterizations of their optical properties in the UV and visible spectral ranges along with improved description of their chemical composition. To address this issue, we used a recently developed approach to retrieve the time- and spectral-dependent optical properties of ambient biomass burning aerosols between 300 and 650 nm wavelength during a regional bonfire festival in Israel. We report that during the biomass burning event, the overall absorption at 400 nm increased by about two orders of magnitude, changing the size-weighted single scattering albedo from a background level of 0.95 to 0.7. Based on the new retrieval method, we provide parameterizations of the wavelength-dependent effective complex refractive index from 350 to 650 nm for freshly emitted and aged biomass burning aerosols. In addition, PM 2.5 filter samples were collected for detailed off-line chemical analysis of the water soluble organics that contribute to light absorption. Nitrophenols were identified as the main organic species responsible for the increased absorption at 400-500 nm. These include species such as 4- nitrocatechol, 4-nitrophenol, nitro-syringol and nitro-guaiacol; oxidation-nitration products of methoxyphenols, known products of lignin pyrolysis. Our findings emphasize the importance of both primary and secondary organic aerosol from biomass burning in absorption of solar radiation and in effective radiative forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Quantification of marine aerosol subgrid variability and its correlation with clouds based on high‐resolution regional modeling

Abstract One limitation of most global climate models (GCMs) is that with the horizontal resolutions they typically employ, they cannot resolve the subgrid variability (SGV) of clouds and aerosols, adding extra uncertainties to the aerosol radiative forcing estimation. To inform the development of an aerosol subgrid variability parameterization, here we analyze the aerosol SGV over the southern Pacific Ocean simulated by the high‐resolution Weather Research and Forecasting model coupled to Chemistry. We find that within a typical GCM grid, the aerosol mass subgrid standard deviation is 15% of the grid‐box mean mass near the surface on a 1 month mean basis. The fraction can increase to 50% in the free troposphere. The relationships between the sea‐salt mass concentration, meteorological variables, and sea‐salt emission rate are investigated in both the clear and cloudy portion. Under clear‐sky conditions, marine aerosol subgrid standard deviation is highly correlated with the standard deviations of vertical velocity, cloud water mixing ratio, and sea‐salt emission rates near the surface. It is also strongly connected to the grid box mean aerosol in the free troposphere (between 2 km and 4 km). In the cloudy area, interstitial sea‐salt aerosol mass concentrations are smaller, but higher correlation is found between the subgrid standard deviations of aerosol mass and vertical velocity. Additionally, we find that decreasing the model grid resolution can reduce the marine aerosol SGV but strengthen the correlations between the aerosol SGV and the total water mixing ratio (sum of water vapor, cloud liquid, and cloud ice mixing ratios).

Lin, Guangxing↗
DOE OSTI2017

Toward low‐cloud‐permitting cloud superparameterization with explicit boundary layer turbulence

Abstract Systematic biases in the representation of boundary layer (BL) clouds are a leading source of uncertainty in climate projections. A variation on superparameterization (SP) called “ultraparameterization” (UP) is developed, in which the grid spacing of the cloud‐resolving models (CRMs) is fine enough (250 × 20 m) to explicitly capture the BL turbulence, associated clouds, and entrainment in a global climate model capable of multiyear simulations. UP is implemented within the Community Atmosphere Model using 2° resolution (∼14,000 embedded CRMs) with one‐moment microphysics. By using a small domain and mean‐state acceleration, UP is computationally feasible today and promising for exascale computers. Short‐duration global UP hindcasts are compared with SP and satellite observations of top‐of‐atmosphere radiation and cloud vertical structure. The most encouraging improvement is a deeper BL and more realistic vertical structure of subtropical stratocumulus (Sc) clouds, due to stronger vertical eddy motions that promote entrainment. Results from 90 day integrations show climatological errors that are competitive with SP, with a significant improvement in the diurnal cycle of offshore Sc liquid water. Ongoing concerns with the current UP implementation include a dim bias for near‐coastal Sc that also occurs less prominently in SP and a bright bias over tropical continental deep convection zones. Nevertheless, UP makes global eddy‐permitting simulation a feasible and interesting alternative to conventionally parameterized GCMs or SP‐GCMs with turbulence parameterizations for studying BL cloud‐climate and cloud‐aerosol feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Performance modeling of microsecond scale biological molecular dynamics simulations on heterogeneous architectures

Performance improvements in biomolecular simulations based on molecular dynamics (MD) codes are widely desired. Unfortunately, the factors, which allowed past performance improvements, particularly the microprocessor clock frequencies, are no longer increasing. Hence, novel software and hardware solutions are being explored for accelerating performance of widely used MD codes. In this paper, we describe our efforts on porting, optimizing and tuning of Large-scale Atomic/Molecular Massively Parallel Simulator, a popular MD framework, on heterogeneous architectures: multi-core processors with graphical processing unit (GPU) accelerators. Our implementation is based on accelerating the most computationally expensive non-bonded interaction terms on the GPUs and overlapping the computation on the CPU and GPUs. This functionality is built on top of message passing interface that allows multi-level parallelism to be extracted even at the workstation level with the multi-core CPUs and allows extension of the implementation on GPU-enabled clusters. We hypothesize that the optimal benefit of heterogeneous architectures for applications will come by utilizing all possible resources (for example, CPU-cores and GPU devices on GPU-enabled clusters). Benchmarks for a range of biomolecular system sizes are provided, and an analysis is performed on four generations of NVIDIA's GPU devices. On GPU-enabled Linux clusters, by overlapping and pipelining computation and communication, we observe up to 10-folds application acceleration in multi-core and multi-GPU environments illustrating significant performance improvements. Detailed analysis of the implementation is presented that allows identification of bottlenecks in algorithm, indicating that code optimization and improvements on GPUs could allow microsecond scale simulation throughput on workstations and inexpensive GPU clusters, putting widely desired biologically relevant simulation time-scales within reach of a large user community. In order to systematically optimize simulation throughput and to enable performance prediction, we have developed a parameterized performance model that will allow developers and users to explore the performance potential of future heterogeneous systems for biological simulations.

Agarwal, Pratul K.↗
DOE OSTI2018

Photodegradation accelerates ecosystem N cycling in a simulated California grassland

Abstract Photodegradation accelerates litter decay in arid grasslands where plant growth and litter decay are strongly controlled by precipitation and evapotranspiration. However, the effects of photodegradation on ecosystem C and N dynamics are not well understood. We examined the effects using an ecosystem biogeochemical model DayCent‐ UV with photodegradation explicitly represented and validated. The model was parameterized for a California grassland where photodegradation was documented to release CO 2 from litter. The model was parameterized with an inverse modeling approach using an extensive data set of six years of daily observed carbon and water gas exchange (gross primary production, ecosystem respiration, and evapotranspiration), soil temperature, and soil moisture. DayCent‐ UV correctly simulated the seasonal patterns of the observed gas exchange and closely simulated the inter‐annual variation in the gas exchange and biomass production rates. The simulations suggested that the inter‐annual variation is driven more by actual evapotranspiration than by precipitation because a large portion of precipitation is lost as runoff during wet years. Photodegradation in DayCent‐ UV accelerated C and N cycling, decreasing system C and N by 9.2% and 9.5% and C and N residence times by 9.4% and 18.2%. Accelerated N cycling made a greater fraction of system N available for plants, increasing net N mineralization and plant production for a given amount of system N. Increased net N mineralization was due to decreased immobilization by microbes in the aboveground organic matter. Photodegradation did not alter the control on plant production by evapotranspiration. These results suggest that at the ecosystem level, the central effect of photodegradation is to suppress microbial activity. We conclude that photodegradation accelerates N cycling at the expense of microbes in this grassland, making it more efficient in supporting plant growth for a given amount of N in the system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Investigating old-growth ponderosa pine physiology using tree-rings, δ 13 C, δ 18 O, and a process-based model

In dealing with predicted changes in environmental conditions outside those experienced today, forest managers and researchers rely on process-based models to inform physiological processes and predict future forest growth responses. The carbon and oxygen isotope ratios of tree-ring cellulose (d13Ccell, d18Ocell) reveal long-term, integrated physiological responses to environmental conditions. We incorporated a submodel of d18Ocell into the widely used Physiological Principles in Predicting Growth (3-PG) model for the first time, to complement a recently added d13Ccell submodel. We parameterized the model using previously reported stand characteristics and long-term trajectories of tree-ring growth, d13Ccell, and d18Ocell collected from the Metolius AmeriFlux site in central Oregon (upland trees). We then applied the parameterized model to a nearby set of riparian trees to investigate the physiological drivers of differences in observed basal area increment (BAI) and d13Ccell trajectories between upland and riparian trees. The model showed that greater available soil water and maximum canopy conductance likely explain the greater observed BAI and lower d13Ccell of riparian trees. Unexpectedly, both observed and simulated d18Ocell trajectories did not differ between the upland and riparian trees, likely due to similar d18O of source water isotope composition. The d18Ocell submodel with a Peclet effect improved model estimates of d18Ocell because its calculation utilizes 3-PG growth and allocation processes. Because simulated stand-level transpiration (E) is used in the d18O submodel, aspects of leaf-level anatomy such as the effective path length for transport of water from the xylem to the sites of evaporation could be estimated.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2019

Improving the SWAT forest module for enhancing water resource projections: A case study in the St. Croix River basin

Abstract Insufficiently calibrated forest parameters of the Soil & Water Assessment Tool (SWAT) may introduce uncertainties to water resource projections in forested watersheds. In this study, we improved SWAT forest parameterization and phosphorus cycling representations to better simulate forest ecosystems in the St. Croix River basin, and we further examined how those improvements affected model projections of streamflow, sediment, and nitrogen export under future climate conditions. Simulations with improved forest parameters substantially reduced model estimates of water, sediment, and nitrogen fluxes relative to those based on default parameters. Differences between improved and default projections can be attributed to the enhanced representation of forest water consumption, nutrient uptake, and protection of soil from erosion. Better representation of forest ecosystems in SWAT contributes to constraining uncertainties in water resource projections. Results of this study highlight the importance of improving SWAT forest ecosystem representations in projecting delivery of water, sediment, and nutrients from land to rivers in response to climate change, particularly for watersheds with large areas of forests. Improved forest parameters and the phosphorus weathering algorithms developed in this study are expected to help enhance future applications of SWAT to investigate hydrological and biogeochemical consequences of climate change.

Yang, Qichun↗
DOE OSTI2015

An analytical bond-order potential for carbon

Carbon is the most widely studied material today because it exhibits special properties not seen in any other materials when in nano dimensions such as nanotube and graphene. Reduction of material defects created during synthesis has become critical to realize the full potential of carbon structures. Molecular dynamics (MD) simulations, in principle, allow defect formation mechanisms to be studied with high fidelity, and can, therefore, help guide experiments for defect reduction. Such MD simulations must satisfy a set of stringent requirements. First, they must employ an interatomic potential formalism that is transferable to a variety of carbon structures. Second, the potential needs to be appropriately parameterized to capture the property trends of important carbon structures, in particular, diamond, graphite, graphene, and nanotubes. The potential must predict the crystalline growth of the correct phases during direct MD simulations of synthesis to achieve a predictive simulation of defect formation. An unlimited number of structures not included in the potential parameterization are encountered, thus the literature carbon potentials are often not sufficient for growth simulations. We have developed an analytical bond order potential for carbon, and have made it available through the public MD simulation package LAMMPS. We also demonstrate that our potential reasonably captures the property trends of important carbon phases. As a result, stringent MD simulations convincingly show that our potential accounts not only for the crystalline growth of graphene, graphite, and carbon nanotubes but also for the transformation of graphite to diamond at high pressure.

36 MATERIALS SCIENCE↗