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. Most importantly, 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. Because an unlimited number of structures not included in the potential parameterization are encountered, 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 demonstrate that our potential reasonably captures the property trends of important carbon phases. 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.

Spatial resolution of the electrical conductance of ionic fluids using a Green-Kubo method.

We present a Green-Kubo method to spatially resolve transport coefficients in compositionally heterogeneous mixtures. We develop the underlying theory based on well-known results from mixture theory, Irving-Kirkwood field estimation, and linear response theory. Then, using standard molecular dynamics techniques, we apply the methodology to representative systems. With a homogeneous salt water system, where the expectation of the distribution of conductivity is clear, we demonstrate the sensitivities of the method to system size, and other physical and algorithmic parameters. Then we present a simple model of an electrochemical double layer where we explore the resolution limit of the method. In this system, we observe significant anisotropy in the wall-normal vs. transverse ionic conductances, as well as near wall effects. Finally, we discuss extensions and applications to more realistic systems such as batteries where detailed understanding of the transport properties in the vicinity of the electrodes is of technological importance.

Melt-growth dynamics in CdTe crystals.

We use a new, quantum-mechanics-based bond-order potential (BOP) to reveal melt growth dynamics and fine scale defect formation mechanisms in CdTe crystals. Previous molecular dynamics simulations of semiconductors have shown qualitatively incorrect behavior due to the lack of an interatomic potential capable of predicting both crystalline growth and property trends of many transitional structures encountered during the melt→crystal transformation. Here, we demonstrate successful molecular dynamics simulations of melt growth in CdTe using a BOP that significantly improves over other potentials on property trends of different phases. Our simulations result in a detailed understanding of defect formation during the melt growth process. Equally important, we show that the new BOP enables defect formation mechanisms to be studied at a scale level comparable to empirical molecular dynamics simulation methods with a fidelity level approaching quantum-mechanical methods.

Accuracy of existing atomic potentials for the CdTe semiconductor compound.

CdTe and CdTe-based Cd(1-x)Zn(x)Te (CZT) alloys are important semiconductor compounds that are used in a variety of technologies including solar cells, radiation detectors, and medical imaging devices. Performance of such systems, however, is limited due to the propensity of nano- and micro-scale defects that form during crystal growth and manufacturing processes. Molecular dynamics simulations offer an effective approach to study the formation and interaction of atomic scale defects in these crystals, and provide insight on how to minimize their concentrations. The success of such a modeling effort relies on the accuracy and transferability of the underlying interatomic potential used in simulations. Such a potential must not only predict a correct trend of structures and energies of a variety of elemental and compound lattices, defects, and surfaces but also capture correct melting behavior and should be capable of simulating crystalline growth during vapor deposition as these processes sample a variety of local configurations. In this paper, we perform a detailed evaluation of the performance of two literature potentials for CdTe, one having the Stillinger-Weber form and the other possessing the Tersoff form. We examine simulations of structures and the corresponding energies of a variety of elemental and compound lattices, defects, and surfaces compared to those obtained from ab initio calculations and experiments. We also perform melting temperature calculations and vapor deposition simulations. Our calculations show that the Stillinger-Weber parameterization produces the correct lowest energy structure. This potential, however, is not sufficiently transferrable for defect studies. Origins of the problems of these potentials are discussed and insights leading to the development of a more transferrable potential suitable for molecular dynamics simulations of defects in CdTe crystals are provided.

Engineering size-scaling of plastic deformation in nanoscale asperities.

Size-dependent plastic flow behavior is manifested in nanoindentation, microbending, and pillar-compression experiments and plays a key role in the contact mechanics and friction of rough surfaces. Recent experiments using a hard flat plate to compress single-crystal Au nano-pyramids and others using a Berkovich indenter to indent flat thin films show size scaling into the 100-nm range where existing mechanistic models are not expected to apply. To bridge the gap between single-dislocation nucleation at the 1-nm scale and dislocation-ensemble plasticity at the 1-microm scale, we use large-scale molecular dynamics (MD) simulations to predict the magnitude and scaling of hardness H versus contact size l(c) in nano-pyramids. Two major results emerge: a regime of near-power-law size scaling H approximately l(c)(-eta) exists, with eta(MD) approximately 0.32 compared with eta(expt) approximately 0.75, and unprecedented quantitative and qualitative agreement between MD and experiments is achieved, with H(MD) approximately 4 GPa at l(c) = 36 nm and H(expt) approximately 2.5 GPa at l(c) = 100 nm. An analytic model, incorporating the energy costs of forming the geometrically necessary dislocation structures that accommodate the deformation, is developed and captures the unique magnitude and size scaling of the hardness at larger MD sizes and up to experimental scales while rationalizing the transition in scaling between MD and experimental scales. The model suggests that dislocation-dislocation interactions dominate at larger scales, whereas the behavior at the smallest MD scales is controlled by nucleation over energy barriers. These results provide a basic framework for understanding and predicting size-dependent plasticity in nanoscale asperities under contact conditions in realistic engineered surfaces.

Reflex inhibition of crural diaphragmatic activity by esophageal distention in cats.

Distention of the esophagus has been shown to result in selective inhibition of phasic inspiratory activity in the crural portion of the diaphragm, with no effect on costal diaphragmatic activity. The purpose of this study was to determine rigourously the afferent pathways that mediate this response. Bipolar EMG electrodes were placed in the costal and crural portions of the diaphragm in decerebrate, spontaneously breathing cats. Distention of the esophagus by inflation of a Foley catheter balloon with 20 ml of air resulted in a selective inhibition of crural hiatal EMG activity, while costal EMG activity was maintained at predistention levels. The distention was accompanied by a reduction in respiratory frequency. Transection of the spinal cord at the C8-T1 level did not obliterate the crural inhibition produced by inflation. Section of the C4-C8 dorsal roots also failed to abolish the response. However, after bilateral cervical vagotomy, esophageal distention no longer influenced diaphragmatic EMG activity. These results indicate that the crural inhibition observed with esophageal distention is vagally mediated and is not influenced importantly by intercostal or phrenic afferents. Records of activity of the phrenic nerve branch innervating the crural portion of the diaphragm showed a similar response pattern, confirming that the inhibition is central in origin and that the crural fibers inhibited by distention are only a fraction of the total population of crural phrenic motoneurons.

Patterns of neural and muscular electrical activity in costal and crural portions of the diaphragm.

To determine whether the central respiratory drives to costal and crural portions of the diaphragm differ from each other in response to chemical and mechanical feedbacks, activities of costal and crural branches of the phrenic nerve were recorded in decerebrate paralyzed cats, studied either with vagi intact and servo-ventilated in accordance with their phrenic nerve activity or vagotomized and ventilated conventionally. Costal and crural electromyograms (EMGs) were recorded in decerebrate spontaneously breathing cats. Hypercapnia and hypoxia resulted in significant increases in peak integrated costal, crural, and whole phrenic nerve activities when the vagi were either intact or cut. However, there were no consistent differences between costal and crural neural responses. Left crural EMG activity was increased significantly more than left costal EMG activity in response to hypercapnia and hypoxia. These results indicate that the central neural inputs to costal and crural portions of the diaphragm are similar in eupnea and in response to chemical and mechanical feedback in decerebrate paralyzed cats. The observed differences in EMG activities in spontaneously breathing animals must arise from modulation of central respiratory activity by mechanoreceptor feedback from respiratory muscles, likely the diaphragm itself.

Influence of extreme hypercapnia on respiratory motor nerve activity in cats.

Sedative drugs have been found to depress the respiratory activity of upper airway muscles more than that of the diaphragm. To determine whether CO2 at narcotic levels has a similar action, we recorded phrenic and hypoglossal nerve activities in decerebrate, vagotomized, paralyzed cats. T5 or T6 external intercostal nerve activity was also recorded in some animals. End-tidal CO2 concentration was raised progressively to over 30% or until depression of nerve activity was apparent. Respiratory frequency was reduced by severe hypercapnia in most cats. Hypoglossal nerve activity was consistently decreased more than that of the phrenic nerve. In most cases intercostal nerve activity was also more susceptible than phrenic nerve activity to hypercapnic depression. The results indicate that CO2 at narcotic levels interferes both with the central pattern generator for breathing movements and with the expression of the pattern in specific motor nerves.

Hypoxia inhibits abdominal expiratory nerve activity.

Our purpose was to examine the influence of steady-state changes in chemical stimuli, as well as discrete peripheral chemoreceptor stimulation, on abdominal expiratory motor activity. In decerebrate, paralyzed, vagotomized, and ventilated cats that had bilateral pneumothoraces, we recorded efferent activity from a phrenic nerve and from an abdominal nerve (cranial iliohypogastric nerve, L1). All cats showed phasic expiratory abdominal nerve discharge at normocapnia [end-tidal PCO2 38 +/- 2 Torr], but small doses (2-6 mg/kg) of pentobarbital sodium markedly depressed this activity. Hyperoxic hypercapnia consistently enhanced abdominal expiratory activity and shortened the burst duration. Isocapnic hypoxia caused inhibition of abdominal nerve discharge in 11 of 13 cats. Carotid sinus nerve denervation (3 cats) exacerbated the hypoxic depression of abdominal nerve activity and depressed phrenic motor output. Stimulation of peripheral chemoreceptors with NaCN increased abdominal nerve discharge in 7 of 10 cats, although 2 cats exhibited marked inhibition. Four cats with intact neuraxis, but anesthetized with ketamine, yielded qualitatively similar results. We conclude that when cats are subjected to steady-state chemical stimuli in isolation (no interference from proprioceptive inputs), hypercapnia potentiates, but hypoxia attenuates, abdominal expiratory nerve activity. Mechanisms to explain the selective inhibition of expiratory motor activity by hypoxia are proposed, and physiological implications are discussed.

The aging face.

Facial aging is a dynamic process which continues throughout adult life. Individuals are affected to a variable degree depending on facial motor habits, exposure and susceptibility to damaging ultraviolet radiation, smoking, and the microscopic tissue changes inherent to the aging process. Aside from avoiding smoking and utilizing sunscreens and other solar protection, there is little that can be done to prevent or retard the development of these changes. There is, however, a vast array of corrective procedures which may be utilized in a planned, sequential fashion to deal with age-related deformities as they occur.