Studies of phonon anharmonicity in solids

Today our understanding of the vibrational thermodynamics of materials at low temperatures is emerging nicely, based on the harmonic model in which phonons are independent. At high temperatures, however, this understanding must accommodate how phonons interact with other phonons or with other excitations. We shall see that the phonon-phonon interactions give rise to interesting coupling problems, and essentially modify the equilibrium and non-equilibrium properties of materials, e.g., thermodynamic stability, heat capacity, optical properties and thermal transport of materials. Despite its great importance, to date the anharmonic lattice dynamics is poorly understood and most studies on lattice dynamics still rely on the harmonic or quasiharmonic models. There have been very few studies on the pure phonon anharmonicity and phonon-phonon interactions. The work presented in this thesis is devoted to the development of experimental and computational methods on this subject.

Modern inelastic scattering techniques with neutrons or photons are ideal for sorting out the anharmonic contribution. Analysis of the experimental data can generate vibrational spectra of the materials, i.e., their phonon densities of states or phonon dispersion relations. We obtained high quality data from laser Raman spectrometer, Fourier transform infrared spectrometer and inelastic neutron spectrometer. With accurate phonon spectra data, we obtained the energy shifts and lifetime broadenings of the interacting phonons, and the vibrational entropies of different materials. The understanding of them then relies on the development of the fundamental theories and the computational methods.

We developed an efficient post-processor for analyzing the anharmonic vibrations from the molecular dynamics (MD) calculations. Currently, most first principles methods are not capable of dealing with strong anharmonicity, because the interactions of phonons are ignored at finite temperatures. Our method adopts the Fourier transformed velocity autocorrelation method to handle the big data of time-dependent atomic velocities from MD calculations, and efficiently reconstructs the phonon DOS and phonon dispersion relations. Our calculations can reproduce the phonon frequency shifts and lifetime broadenings very well at various temperatures.

To understand non-harmonic interactions in a microscopic way, we have developed a numerical fitting method to analyze the decay channels of phonon-phonon interactions. Based on the quantum perturbation theory of many-body interactions, this method is used to calculate the three-phonon and four-phonon kinematics subject to the conservation of energy and momentum, taking into account the weight of phonon couplings. We can assess the strengths of phonon-phonon interactions of different channels and anharmonic orders with the calculated two-phonon DOS. This method, with high computational efficiency, is a promising direction to advance our understandings of non-harmonic lattice dynamics and thermal transport properties.

These experimental techniques and theoretical methods have been successfully performed in the study of anharmonic behaviors of metal oxides, including rutile and cuprite stuctures, and will be discussed in detail in Chapters 4 to 6. For example, for rutile titanium dioxide (TiO₂), we found that the anomalous anharmonic behavior of the B_1g mode can be explained by the volume effects on quasiharmonic force constants, and by the explicit cubic and quartic anharmonicity. For rutile tin dioxide (SnO₂), the broadening of the B_2g mode with temperature showed an unusual concave downwards curvature. This curvature was caused by a change with temperature in the number of down-conversion decay channels, originating with the wide band gap in the phonon dispersions. For silver oxide (Ag₂O), strong anharmonic effects were found for both phonons and for the negative thermal expansion.

Chesapeake and Delaware Canal Affects Environmental: presented at American Society of Civil Engineers National Water Resources Engineering Meeting, Phoenix, Arizona, 14 January 1971

The Chesapeake and Delaware Canal is a man-made waterway connecting the upper Chesapeake Bay with the Delaware Bay. It started in 1829 as a private barge canal with locks, two at the Delaware end, and one at the Chesapeake end. For the most part, natural tidal and non-tidal waterways were connected by short dredged sections to form the original canal. In 1927, the C and D Canal was converted to a sea-level canal, with a controlling depth of 14 feet, and a width of 150 feet. In 1938 the canal was deepened to 27 feet, with a channel width of 250 feet. Channel side slopes were dredged at 2.5:1, thus making the total width of the waterway at least 385 feet in those segments representing new cuts or having shore spoil area dykes rising above sea level. In 1954 Congress authorized a further enlargement of the Canal to a depth of 35 feet and a channel width of 450 feet. (pdf contains 27 pages)

Oyster spat set on natural clutch in the Maryland portion of the Chesapeake Bay (1939 - 1975)

Basically this report is an attempt to document trends in oyster recruitment since 1939 and to relate those trends to the actual oyster harvest throughout the Maryland portion of the Chesapeake Bay. It is also hoped that the data as well as the charts compiled in this report will serve as a reference to aid in future studies on Chesapeake Bay oysters. A few if the major biological factors that affect the natural reproduction of the oyster and environmental degradations that may possibly affect oyster reproduction or harvest in the Chesapeake Bay are also briefly discussed. (PDF contains 32 pages)

Methods for melting temperature calculation

Melting temperature calculation has important applications in the theoretical study of phase diagrams and computational materials screenings. In this thesis, we present two new methods, i.e., the improved Widom's particle insertion method and the small-cell coexistence method, which we developed in order to capture melting temperatures both accurately and quickly.

We propose a scheme that drastically improves the efficiency of Widom's particle insertion method by efficiently sampling cavities while calculating the integrals providing the chemical potentials of a physical system. This idea enables us to calculate chemical potentials of liquids directly from first-principles without the help of any reference system, which is necessary in the commonly used thermodynamic integration method. As an example, we apply our scheme, combined with the density functional formalism, to the calculation of the chemical potential of liquid copper. The calculated chemical potential is further used to locate the melting temperature. The calculated results closely agree with experiments.

We propose the small-cell coexistence method based on the statistical analysis of small-size coexistence MD simulations. It eliminates the risk of a metastable superheated solid in the fast-heating method, while also significantly reducing the computer cost relative to the traditional large-scale coexistence method. Using empirical potentials, we validate the method and systematically study the finite-size effect on the calculated melting points. The method converges to the exact result in the limit of a large system size. An accuracy within 100 K in melting temperature is usually achieved when the simulation contains more than 100 atoms. DFT examples of Tantalum, high-pressure Sodium, and ionic material NaCl are shown to demonstrate the accuracy and flexibility of the method in its practical applications. The method serves as a promising approach for large-scale automated material screening in which the melting temperature is a design criterion.

We present in detail two examples of refractory materials. First, we demonstrate how key material properties that provide guidance in the design of refractory materials can be accurately determined via ab initio thermodynamic calculations in conjunction with experimental techniques based on synchrotron X-ray diffraction and thermal analysis under laser-heated aerodynamic levitation. The properties considered include melting point, heat of fusion, heat capacity, thermal expansion coefficients, thermal stability, and sublattice disordering, as illustrated in a motivating example of lanthanum zirconate (La₂Zr₂O₇). The close agreement with experiment in the known but structurally complex compound La₂Zr₂O₇ provides good indication that the computation methods described can be used within a computational screening framework to identify novel refractory materials. Second, we report an extensive investigation into the melting temperatures of the Hf-C and Hf-Ta-C systems using ab initio calculations. With melting points above 4000 K, hafnium carbide (HfC) and tantalum carbide (TaC) are among the most refractory binary compounds known to date. Their mixture, with a general formula Ta_xHf_1-xC_y, is known to have a melting point of 4215 K at the composition Ta₄HfC₅, which has long been considered as the highest melting temperature for any solid. Very few measurements of melting point in tantalum and hafnium carbides have been documented, because of the obvious experimental difficulties at extreme temperatures. The investigation lets us identify three major chemical factors that contribute to the high melting temperatures. Based on these three factors, we propose and explore a new class of materials, which, according to our ab initio calculations, may possess even higher melting temperatures than Ta-Hf-C. This example also demonstrates the feasibility of materials screening and discovery via ab initio calculations for the optimization of "higher-level" properties whose determination requires extensive sampling of atomic configuration space.

Dendrites Inhibition in Rechargeable Lithium Metal Batteries

The specific high energy and power capacities of rechargeable lithium metal (Li⁰) batteries are ideally suited to portable devices and are valuable as storage units for intermittent renewable energy sources. Lithium, the lightest and most electropositive metal, would be the optimal anode material for rechargeable batteries if it were not for the fact that such devices fail unexpectedly by short-circuiting via the dendrites that grow across electrodes upon recharging. This phenomenon poses a major safety issue because it triggers a series of adverse events that start with overheating, potentially followed by the thermal decomposition and ultimately the ignition of the organic solvents used in such devices.

In this thesis, we developed experimental platform for monitoring and quantifying the dendrite populations grown in a Li battery prototype upon charging under various conditions. We explored the effects of pulse charging in the kHz range and temperature on dendrite growth, and also on loss capacity into detached “dead” lithium particles.

Simultaneously, we developed a computational framework for understanding the dynamics of dendrite propagation. The coarse-grained Monte Carlo model assisted us in the interpretation of pulsing experiments, whereas MD calculations provided insights into the mechanism of dendrites thermal relaxation. We also developed a computational framework for measuring the dead lithium crystals from the experimental images.

Mortality, growth and growth rate variability of striped bass larvae in Chesapeake Subestuaries

In this report we develop age-length keys and derive age-frequency data. We estimate striped bass and white perch mortality and growth rates, based on the otolith-aging analysis. We also report on hatch-date frequencies of striped bass and white perch larvae, and we discuss environmental effects on recruitment potential.

Effects of suspended sediments on the development of eggs and larvae of striped bass and white perch

The possible ecological effects of suspended sediments are manifold. Briefly, suspended sediments may cause an increased surface for microorganism growth, fewer temperature fluctuations, chemical adsorption or absorption, blanketing, mechanical-abrasive actions, and light penetration reduction (Cairns, 1968). Sherk and Cronin (1970) have pointed out that the above effects have been little studied in the estuarine environment. The ecological effects of suspended sediments on fish eggs and larvae may be of prime importance t o the C and D Canal area, an important spawning and primary nursery area for a variety of estuary: e species (Johnson,1972). This section discusses the effects of suspended sediment on the eggs and larvae of striped bass and white perch.

Eggs, Larvae and Young of The Striped Bass, Roccus saxatilis

Detailed descriptions of the early development of the striped bass, Roccus saxitilis (Walbaum), with emphasis on variation in size and morphology, sequence of fin formation, changes in body form, and attainment of the full complement of maristic numbers, are presented and illustrated for the first time. The egg is spherical, transparent, non-adhesive and relatively large. It is pelagic and buoyant, although it sinks in quiet fresh water. When unfertilized, it averages 1.3 mm, in diameter, but is 3.4 mm. when fertilized and water-hardened. The granular yolk sac, green when alive and whitish-yellow when preserved, averages 1.2 mm., and the single amber-colored oil globule is about 0.6 mm. in diameter. Newly hatched striped bass prolarvae, which range from 2.9-3.7 mm. in total length, are relatively undeveloped and nearly transparent, with no mouth opening, unpigmented eyes, and a greatly enlarged yolk sac with the large oil globule projecting beyond the head. When 5-6 mm. long the yolk sac and oil globule are assimilated and the postlarvae I show advanced development of the internal anatomy. Although the fish is still transparent, scattered melanophores are found on the head and body and chromatophores in the eyes and the ventro-posterior edge of the body. Postlarvae transform to young between 7 and 10 mm. in length when the finfolds are lost except in the dorsal, anal and caudal regions. The largest fish in this group possess a well-formed skeleton with a full complement of 25 vertebrae. Between 10 and 20 mm. in length all fish are fully transformed, muscular tissue renders most of the internal structure obscure, and the myotomes, which generally correspond in number with the vertebrae, are no longer visible. At fish lengths of 20-30 mm. scales are found on all specimens, and with the exception of the pectoral fin-rays, a full complement of meristic structures is present in all other fins. At this stage the body is pigmented uniformly with small spots. Linear regressions between several dependent variables and the , independent variable of standard length indicate that the rate of development of head, eye. and snout to anus lengths is proportional to the length of the larvae and young. Body depth and standard length are non-linear among newly-hatched larvae. Hatchery-reared striped bass demonstrated a slow rate of growth, and were regarded as "stunted," when compared to growth rates observed in another study and field collections. Observations were also made on abnormal eggs and teratological larvae and young. Blue-sac disease is tentatively identified and described for the first time in larvae and pugnosed larvae and young are also described for the first time in striped bass.

Changes in Fish Communities in the Upper Patuxent River from 1966 to 1977

Ten year comparison of fish survey's with respect to diversity evenness and composition of fish communities. The upper Patuxent River was divided into Piedmont Plateau and Coastal Plain regions, not only for geographical purposes, but also because of the clustering of sewage treatment plants in the Coastal Plain region. In the Piedmont Plateau region, the fish species diversity changed very little from 1966 to 1977 ( Little Patuxent -- 2.82 to 2.66; Middle Patuxent -- 2.86 to 2.83; and main stem -- 2.46 to 2.63), except in a section of Little Patuxent River at and below the City of Columbia where the species diversity index showed a significant reduction from 2.97 to 1.99, and in a section of the main stem Patuxent River immediately downstream from the Brighton Dam of the Triadelphia Reservoir where the index increased significantly from 1.66 to 3.20. In the Coastal Plain region, a significant reduction in the fish species diversity index occurred between 1966 and 1977 below the two sewage treatment plant outfalls : Savage -- 2.69 to 0 and Patuxent-Crofton -- 3.06 to 1.33. Also, the substantial reduction in the species diversity index which had already occurred in 1966 below the six other plant: outfalls of Fort Meade No. I, Fort Meade No. 2, Maryland House of Correction, Maryland City , Parkway and Bowie, remained depressed in 1977. On the other hand, below the Horsepen Sewage Treatment Plant (a tertiary plant practicing dechlorination) the species diversity index increased from 1.91 to 2.8. (PDF contains 48 pages)