Relating thermodynamics to information theory : the equality of free energy and mutual information

In this thesis we uncover a new relation which links thermodynamics and information theory. We consider time as a channel and the detailed state of a physical system as a message. As the system evolves with time, ever present noise insures that the "message" is corrupted. Thermodynamic free energy measures the approach of the system toward equilibrium. Information theoretical mutual information measures the loss of memory of initial state. We regard the free energy and the mutual information as operators which map probability distributions over state space to real numbers. In the limit of long times, we show how the free energy operator and the mutual information operator asymptotically attain a very simple relationship to one another. This relationship is founded on the common appearance of entropy in the two operators and on an identity between internal energy and conditional entropy. The use of conditional entropy is what distinguishes our approach from previous efforts to relate thermodynamics and information theory.

Non-linear dispersive waves with a small dissipation

The general theory of Whitham for slowly-varying non-linear wavetrains is extended to the case where some of the defining partial differential equations cannot be put into conservation form. Typical examples are considered in plasma dynamics and water waves in which the lack of a conservation form is due to dissipation; an additional non-conservative element, the presence of an external force, is treated for the plasma dynamics example. Certain numerical solutions of the water waves problem (the Korteweg-de Vries equation with dissipation) are considered and compared with perturbation expansions about the linearized solution; it is found that the first correction term in the perturbation expansion is an excellent qualitative indicator of the deviation of the dissipative decay rate from linearity.

A method for deriving necessary and sufficient conditions for the existence of a general uniform wavetrain solution is presented and illustrated in the plasma dynamics problem. Peaking of the plasma wave is demonstrated, and it is shown that the necessary and sufficient existence conditions are essentially equivalent to the statement that no wave may have an amplitude larger than the peaked wave.

A new type of fully non-linear stability criterion is developed for the plasma uniform wavetrain. It is shown explicitly that this wavetrain is stable in the near-linear limit. The nature of this new type of stability is discussed.

Steady shock solutions are also considered. By a quite general method, it is demonstrated that the plasma equations studied here have no steady shock solutions whatsoever. A special type of steady shock is proposed, in which a uniform wavetrain joins across a jump discontinuity to a constant state. Such shocks may indeed exist for the Korteweg-de Vries equation, but are barred from the plasma problem because entropy would decrease across the shock front.

Finally, a way of including the Landau damping mechanism in the plasma equations is given. It involves putting in a dissipation term of convolution integral form, and parallels a similar approach of Whitham in water wave theory. An important application of this would be towards resolving long-standing difficulties about the "collisionless" shock.

Conditional independence in quantum many-body systems

In this thesis, I will discuss how information-theoretic arguments can be used to produce sharp bounds in the studies of quantum many-body systems. The main advantage of this approach, as opposed to the conventional field-theoretic argument, is that it depends very little on the precise form of the Hamiltonian. The main idea behind this thesis lies on a number of results concerning the structure of quantum states that are conditionally independent. Depending on the application, some of these statements are generalized to quantum states that are approximately conditionally independent. These structures can be readily used in the studies of gapped quantum many-body systems, especially for the ones in two spatial dimensions. A number of rigorous results are derived, including (i) a universal upper bound for a maximal number of topologically protected states that is expressed in terms of the topological entanglement entropy, (ii) a first-order perturbation bound for the topological entanglement entropy that decays superpolynomially with the size of the subsystem, and (iii) a correlation bound between an arbitrary local operator and a topological operator constructed from a set of local reduced density matrices. I also introduce exactly solvable models supported on a three-dimensional lattice that can be used as a reliable quantum memory.

Predicting the effects of climate change on sea turtle nesting habitat in Florida

Rising global temperatures threaten the survival of many plant and animal species. Having already risen at an unprecedented rate in the past century, temperatures are predicted to rise between 0.3 and 7.5C in North America over the next 100 years (Hawkes et al. 2007). Studies have documented the effects of climate warming on phenology (timing of seasonal activities), with observations of early arrival at breeding grounds, earlier ends to the reproductive season, and delayed autumnal migrations (Pike et al. 2006). In addition, for species not suited to the physiological demands of cold winter temperatures, increasing temperatures could shift tolerable habitats to higher latitudes (Hawkes et al. 2007). More directly, climate warming will impact thermally sensitive species like sea turtles, who exhibit temperature-dependent sexual determination. Temperatures in the middle third of the incubation period determine the sex of sea turtle offspring, with higher temperatures resulting in a greater abundance of female offspring. Consequently, increasing temperatures from climate warming would drastically change the offspring sex ratio (Hawkes et al. 2007). Of the seven extant species of sea turtles, three (leatherback, Kemp’s ridley, and hawksbill) are critically endangered, two (olive ridley and green) are endangered, and one (loggerhead) is threatened. Considering the predicted scenarios of climate warming and the already tenuous status of sea turtle populations, it is essential that efforts are made to understand how increasing temperatures may affect sea turtle populations and how these species might adapt in the face of such changes. In this analysis, I seek to identify the impact of changing climate conditions over the next 50 years on the availability of sea turtle nesting habitat in Florida given predicted changes in temperature and precipitation. I predict that future conditions in Florida will be less suitable for sea turtle nesting during the historic nesting season. This may imply that sea turtles will nest at a different time of year, in more northern latitudes, to a lesser extent, or possibly not at all. It seems likely that changes in temperature and precipitation patterns will alter the distribution of sea turtle nesting locations worldwide, provided that beaches where the conditions are suitable for nesting still exist. Hijmans and Graham (2006) evaluate a range of climate envelope models in terms of their ability to predict species distributions under climate change scenarios. Their results suggested that the choice of species distribution model is dependent on the specifics of each individual study. Fuller et al. (2008) used a maximum entropy approach to model the potential distribution of 11 species in the Arctic Coastal Plain of Alaska under a series of projected climate scenarios. Recently, Pike (in press) developed Maxent models to investigate the impacts of climate change on green sea turtle nest distribution and timing. In each of these studies, a set of environmental predictor variables (including climate variables), for which ‘current’ conditions are available and ‘future’ conditions have been projected, is used in conjunction with species occurrence data to map potential species distribution under the projected conditions. In this study, I will take a similar approach in mapping the potential sea turtle nesting habitat in Florida by developing a Maxent model based on environmental and climate data and projecting the model for future climate data. (PDF contains 5 pages)

Electronic structure and phonon thermodynamics of iron alloys

Inelastic neutron scattering (INS) and nuclear-resonant inelastic x-ray scattering (NRIXS) were used to measure phonon spectra of FeV as a B2- ordered compound and as a bcc solid solution. Contrary to the behavior of ordering alloys studied to date, the phonons in the B2-ordered phase are softer than in the solid solution. Ordering increases the vibrational entropy, which stabilizes the ordered phase to higher temperatures. Ab initio calculations show that the number of electronic states at the Fermi level increases upon ordering, enhancing the screening between ions, and reducing the interatomic force constants. The effect of screening is larger at the V atomic sites than at the Fe atomic sites.

The phonon spectra of Au-rich alloys of fcc Au-Fe were also measured. The main effect on the vibrational entropy of alloying comes from a stiffening of the Au partial phonon density of states (DOS) with Fe concentration that increases the miscibility gap temperature. The magnitude of the effect is non- linear and it is reduced at higher Fe concentrations. Force constants were calculated for several compositions and show a local stiffening of Au–Au bonds close to Fe atoms, but Au–Au bonds that are farther away do not show this effect. Phonon DOS curves calculated from the force constants reproduced the experimental trends. The Au–Fe bond is soft and favors ordering, but a charge transfer from the Fe to the Au atoms stiffens the Au–Au bonds enough to favor unmixing. The stiffening is attributed to two main effects comparable in magnitude: an increase in electron density in the free-electron-like states, and stronger sd-hybridization.

INS and NRIXS measurements were performed at elevated temperatures on B2-ordered FeTi and NRIXS measurements were performed at high pressures. The high-pressure behavior is quasi- harmonic. The softening of the phonon DOS curves with temperature is strongly nonharmonic. Calculations of the force constants and Born-von Karman fits to the experimental data show that the bonds between second nearest neighbors (2nn) are much stiffer than those between 1nn, but fits to the high temperature data show that the former softens at a faster rate with temperature. The Fe–Fe bond softens more than the Ti–Ti bond. The unusual stiffness of the 2nn bond is explained by the calculated charge distribution, which is highly aspherical and localized preferentially in the t2g orbitals. Ab initio molecular dynamics (AIMD) simulations show a charge transfer from the t2g orbitals to the eg orbitals at elevated temperatures. The asphericity decreases linearly with temperature and is more severe at the Fe sites.

Reactivity of titanocene methylidene with metal halides, alkene sulfides, and alkene oxides

Titanocene metallacyclobutanes show a wide variety of reactivites with organic and inorganic reagents. Their reactions include methylene transfer to organic carbonyls, formation of enolates, electron transfer from activated alkyl chlorides, olefin metathesis, ring opening polymerization. Recently, preparations of heterobinuclear µ-methylene complexes were reported. In this thesis, mechanistic, synthetic, and structural studies of the heterobinuclear µ-methylene complexes will be described. Also, the reaction of titanocene methylidene trimethylphosphine complex with alkene sulfide and styrene sulfide will be presented.

Heterobinuclear µ-methylene-µ-methyl complexes C_(p2)Ti(µ-CH_2)( µ-CH_3)M(1,5-COD) have been prepared (M = Rh, Ir). X-ray crystallography showed that the methyl group of the complex was bonded to the rhodium and bridges to the titanium through an agostic bond. The ^(1)H,^(13)CNMR, IR spectra along with partial deuteration studies supported the structure in both solution and solid state. Activation of the agostic bond is demonstrated by the equilibration of the µ-CH_3 and µ-CH_2 groups. A nonlinear Arrhenius plot, an unusually large kinetic isotope effect (24(5)), and a large negative activation entropy (-64(3)eu) can be explained by the quantum-mechanical tunneling. Calculated rate constants with Bell-type barrier fitted well with the observed one. This equilibration was best explained by a 4e-4c mechanism (or σ bond metathesis) with the character of quantum-mechanical tunneling.

Heterobinuclear µ-methylene-µ-phenyl complexes were synthesized. Structural study of C_(p2)Ti(µ-CH_(2))(µ-p-Me_(2)NC_(6)H_(4))Rh(l,5-COD) showed that the two metal atoms are bridged by the methylene carbon and the ipso carbon of the p-N,N-dimethylarninophenyl group. The analogous structure of C_(p2))Ti(µ-CH_(2))(µ-o-MeOC_(6)H_(4))Rh(1,5-COD) has been verified by the differential NOE. The aromaticity of the phenyl group observed by ^(1)H NMR, was confirmed by the comparison of the C-C bond lengths in the crystallographic structure. The unusual downfield shifts of the ipso carbon in the ^(13)C NMR are assumed to be an indication of the interaction between the ipso carbon and electron-deficient titanium.

Titanium-platinum heterobinuclear µ-methylene complexes C_(p2)Ti(µ-CH_(2))(µ -X)Pt(Me)(PM_(2)Ph) have been prepared (X= Cl, Me). Structural studies indicate the following:(1) the Ti-CH2 bond possesses residual double bond character, (2) there is a dative Pt→Ti interaction which may be regarded as a π back donation from the platinum atom to the 'Ti=CH_(2)'' group, and (3) the µ-CH_3 group is bound to the titanium atom through a three-center, two-electron agostic bond.

Titanocene (η^(2)-thioformaldehyde)•PMe_3 was prepared from C_(p2)Ti=CH_(2)•PMe_3 and sulfur-containing organic compounds (e.g. alkene sulfide, triphenylphosphine sulfide) including elemental sulfur. Mechanistic studies utilizing trans-styrene sulfide-d_1 suggested the stepwise reaction to explain equimolar mixture of trans- and cis-styrene-d_1 as by-products. The product reacted with methyl iodide to produce cationic titanocene (η_(2)-thiomethoxymethyl) complex. Complexes having less coordinating anion like BF_4 or BPh_4 could be obtained through metathesis. Together with structural analyses, the further reactivities of the complexes have been explored.

The complex C_(p2)TiOCH_(2)CH(Ph)CH_2 was prepared from the compound C_(p2)Ti=CH_(2)-PMe_3 and styrene oxide. The product was characterized with ^(1)H-^(1)H correlated 2-dimensional NMR, selective decoupling of ^(1)H NMR, and differential NOE. Stereospecificity of deuterium in the product was lost when trans-styrene oxide-d_1 was allowed to react. Relative rates of the reaction were measured with varying substituents on the phenyl ring. Better linearity (r = -0.98, p^(+) = -0.79) was observed with σ_(p)^(+)than σ(r = -0.87, p = -1.26). The small magnitude of p^+ value and stereospecificity loss during the formation of product were best explained by the generation of biradicals, but partial generation of charge cannot be excluded. Carbonylation of the product followed by exposure to iodine yields the corresponding β-phenyl γ-lactone.

Icthyofauna of Oramiriukwa River in Imo State, Nigeria

Oramiriukwa one of the rivers in the Imo River drainage basin runs approximately a 14.0km course from its headwaters in Abba Isu (Nwangele LGA) and discharge into Otamiri River. Based on the fish catches of local fishermen using a variety of fishing gear in a sampling station in Emekuku (Owerri LGA), a total of 17 families of fish consisting of 27 species were identified in the river. In the January sample the most numerically abundant species were Tilapia zilli (12.5%). Calamoichthys calabaricus (8.7%) and T. guineensis (8.6%). In the July samples the most numerically abundant species were Mastacemblus loenbergii (9.2%). T. zilli (7.2%) and Hemichromis fasciatus (6.6%). The Shannon Weiner Diversity index for Oranmiriukwa River was high (1.305) implying a relatively diversified and equilibrated fish community

Fishery resource surveys: beyond length-weight relationship models

Fisheries resource surveys are regular management tools for rational exploitation of commercial fisheries. In a growing number of cases, the use of these resource surveys has been largely restricted to assessment of the relative well being of fish stocks and the potential yields of such fisheries. This paper seeks to demonstrate that the data from such surveys can also be easily used to evaluate species diversity of such fisheries, both in terms of species richness and equitability of distribution. Using published data on two freshwater and two marine fisheries as case studies, Shannon-Wiener Diversity Function and Simpson's Index were computed for each of these fisheries. These biodiversity indices gave a deeper insight into the environmental status of each of these fisheries, beyond what the length-weight relationship models can reveal. Generally, while the marine fisheries showed more species richness, the freshwater fisheries apparently had more stable and equilibrated fish communities

Fermionic quantum systems. Part I: Phase transitions in quantum dots. Part II: Nuclear matter on a lattice

In the first part I perform Hartree-Fock calculations to show that quantum dots (i.e., two-dimensional systems of up to twenty interacting electrons in an external parabolic potential) undergo a gradual transition to a spin-polarized Wigner crystal with increasing magnetic field strength. The phase diagram and ground state energies have been determined. I tried to improve the ground state of the Wigner crystal by introducing a Jastrow ansatz for the wave function and performing a variational Monte Carlo calculation. The existence of so called magic numbers was also investigated. Finally, I also calculated the heat capacity associated with the rotational degree of freedom of deformed many-body states and suggest an experimental method to detect Wigner crystals.

The second part of the thesis investigates infinite nuclear matter on a cubic lattice. The exact thermal formalism describes nucleons with a Hamiltonian that accommodates on-site and next-neighbor parts of the central, spin-exchange and isospin-exchange interaction. Using auxiliary field Monte Carlo methods, I show that energy and basic saturation properties of nuclear matter can be reproduced. A first order phase transition from an uncorrelated Fermi gas to a clustered system is observed by computing mechanical and thermodynamical quantities such as compressibility, heat capacity, entropy and grand potential. The structure of the clusters is investigated with the help two-body correlations. I compare symmetry energy and first sound velocities with literature and find reasonable agreement. I also calculate the energy of pure neutron matter and search for a similar phase transition, but the survey is restricted by the infamous Monte Carlo sign problem. Also, a regularization scheme to extract potential parameters from scattering lengths and effective ranges is investigated.

Estudo do efeito magnetocalórico em sistemas magnéticos com terras raras

O efeito magnetocalórico, base da refrigeração magnética, é caracterizado por duas quantidades: a variação isotérmica da entropia (ΔST) e a variação adiabática da temperatura (ΔTS); que são obtidas sob variações na intensidade de um campo magnético aplicado. Em sistemas que apresentam anisotropia magnética, pode‐se definir o efeito magnetocalórico anisotrópico, o qual, por definição, é calculado sob variações na direção de aplicação de um campo magnético cuja intensidade mantém‐se fixa, e é caracterizado por duas quantidades: a variação anisotrópico‐isotérmica da entropia (ΔSan) e a variação anisotrópico‐adiabática da temperatura (ΔTan). O efeito magnetocalórico e o efeito magnetocalórico anisotrópico foram estudados nos compostos intermetálicos formados por terras e outros materiais não magnéticos: RNi2, RNi5, RZn e Gd1‐nPrnAl2. Os cálculos foram feitos partindo de hamiltonianos modelo que incluem as interações de troca, Zeeman, de campo cristalino e quadrupolar.

Composição florística e estrutura da comunidade arbustivo-arbórea de um trecho de floresta submontana na vertente sudeste do Parque Estadual da Ilha Grande, Angra dos Reis / RJ

O estudo objetivou avaliar a composição florística e estrutural dos componentes arbustivo-arbóreo da Floresta Ombrófila Densa submontana em diferentes estágios de regeneração natural, na vertente sudeste do Parque Estadual da Ilha Grande/RJ. Para o inventário florístico foram realizadas coletas assistemáticas em diferentes trechos nessa vertente. A complementação da lista de espécies foi feita a partir, da consulta às exsicatas dos herbários do Rio de Janeiro (FCAB, GUA, HB, HRJ, R, RB, RBR, RFA, RFFP e RUSU) e do inventário fitossociológico. Foi verificado o status de conservação das espécies inventariadas para a Flora Brasileira. Para o inventário fitossociológico foram estabelecidas 34 parcelas amostrais, totalizando 1,02 ha de área amostrada. Todos os indivíduos arbustivo-arbóreos com DAP ≥ 5 cm foram registrados e, após identificação, foram depositados no Herbário da Universidade do Estado do Rio de Janeiro (HRJ). O pacote estatístico FITOPAC 2.1. foi utilizado para a análise dos dados. A similaridade entre o remanescente investigado neste estudo e as outras quatorze áreas distintas do Rio de Janeiro, da própria Ilha Grande ou não, foi avaliada, utilizando-se o coeficiente de Similaridade de Sorensen; pelo critério de agrupamento por ligação média não ponderada (UPGMA) e pelo método de autorreamostragem para a estrutura de grupos; utilizados os programas PAST v1.34 e Multiv 2.4. A partir do levantamento em herbários e dos inventários florístico e fitossociológico realizados neste trabalho, foram analisados 3.470 registros, sendo 1.778 do levantamento de herbários, 1.536 do levantamento fitossociológico e 156 do inventário florístico. Esses registros corresponderam a 606 espécies ou morfo-espécies de Angiospermas e uma de Pteridófita. Os resultados obtidos revelaram a existência de 22 espécies ameaçadas de extinção para a Flora do Brasil. Dentre, as quais, sete são exclusivas da amostragem fitossociológica: Abarema cochliacarpos (Gomes) Barneby & J.W. Grimes, Chrysophyllum flexuosum Mart., Ficus pulchella Schott ex Spreng., Macrotorus utriculatus Perkins, Myrceugenia myrcioides (Cambess.) O.Berg, Rudgea interrupta Benth e Urbanodendron bahiense (Meisn.) Rohwer. No estudo fitossociológico, inventariou-se 1.536 indivíduos de 217 espécies, subordinadas a 53 famílias. O índice de diversidade de Shannon (H) calculado foi de 4,702 nats/ind e equabilidade (J) de 0,874. As 10 famílias com maior riqueza foram: Myrtaceae (31 spp.), Rubiaceae (21), Fabaceae (17), Lauraceae (12), Euphorbiaceae (11), Monimiaceae (8), Melastomataceae (7), Sapindaceae (7), Sapotaceae (6) e Annonaceae (6). Os 10 maiores Valores de Importância das espécies foram para Chrysophyllum flexuosum (3,43%), Lamanonia ternata Vell. (3,40%), Hyeronima alchorneoides Allemão (2,83%), Actinostemon verticillatus (Klotzsch) Baill. (2,55%), Psychotria brasiliensis Vell. (2,55%), Eriotheca pentaphylla (Vell.) A. Robyns (2,28%), Guatteria australis A. St.-Hil. (2,12%), Mabea brasiliensis Müll. Arg. (2,04%), Miconia prasina (Sw.) DC. (1,89%) e Rustia formosa (Cham. & Schltdl. ex DC.) Klotzsch (1,82%). Amostraram-se 27% de espécies representadas por apenas um indivíduo. As análises florísticas avaliadas a partir do Índice de Similaridade de Sorensen indicaram como principais variáveis para a formação dos blocos, os diferentes valores de diversidade para as áreas e a distribuição fitogeográfica das espécies. Os resultados obtidos junto aos dados dos grupos ecológicos, para os indivíduos da fitossociologia, indicaram maior percentual de indivíduos secundários tardios amostrados. Conclui-se que a área de estudo é uma floresta secundária em estágio intermediário de regeneração, com grande riqueza de espécies, muitas das quais de relevante importância ecológica.

Enhanced Thermoelectric Performance at the Superionic Phase Transitions of Mixed Ion-Electron Conducting Materials

The quality of a thermoelectric material is judged by the size of its temperature de- pendent thermoeletric-figure-of-merit (zT ). Superionic materials, particularly Zn₄Sb₃ and Cu₂Se, are of current interest for the high zT and low thermal conductivity of their disordered, superionic phase. In this work it is reported that the super-ionic materials Ag₂Se, Cu₂Se and Cu_1.97Ag_0.03Se show enhanced zT in their ordered, normal ion-conducting phases. The zT of Ag₂Se is increased by 30% in its ordered phase as compared to its disordered phase, as measured just below and above its first order phase transition. The zT ’s of Cu₂Se and Cu_1.97Ag_0.03Se both increase by more than 100% over a 30 K temperatures range just below their super-ionic phase transitions. The peak zT of Cu₂Se is 0.7 at 406 K and of Cu_1.97Ag_0.03Se is 1.0 at 400 K. In all three materials these enhancements are due to anomalous increases in their Seebeck coefficients, beyond that predicted by carrier concentration measurements and band structure modeling. As the Seebeck coefficient is the entropy transported per carrier, this suggests that there is an additional quantity of entropy co-transported with charge carriers. Such co-transport has been previously observed via co-transport of vibrational entropy in bipolaron conductors and spin-state entropy in Na_xCo₂O₄. The correlation of the temperature profile of the increases in each material with the nature of their phase transitions indicates that the entropy is associated with the thermodynamcis of ion-ordering. This suggests a new mechanism by which high thermoelectric performance may be understood and engineered.

Macroscopically Dissipative Systems with Underlying Microscopic Dynamics : Properties and Limits of Measurement

While some of the deepest results in nature are those that give explicit bounds between important physical quantities, some of the most intriguing and celebrated of such bounds come from fields where there is still a great deal of disagreement and confusion regarding even the most fundamental aspects of the theories. For example, in quantum mechanics, there is still no complete consensus as to whether the limitations associated with Heisenberg's Uncertainty Principle derive from an inherent randomness in physics, or rather from limitations in the measurement process itself, resulting from phenomena like back action. Likewise, the second law of thermodynamics makes a statement regarding the increase in entropy of closed systems, yet the theory itself has neither a universally-accepted definition of equilibrium, nor an adequate explanation of how a system with underlying microscopically Hamiltonian dynamics (reversible) settles into a fixed distribution.

Motivated by these physical theories, and perhaps their inconsistencies, in this thesis we use dynamical systems theory to investigate how the very simplest of systems, even with no physical constraints, are characterized by bounds that give limits to the ability to make measurements on them. Using an existing interpretation, we start by examining how dissipative systems can be viewed as high-dimensional lossless systems, and how taking this view necessarily implies the existence of a noise process that results from the uncertainty in the initial system state. This fluctuation-dissipation result plays a central role in a measurement model that we examine, in particular describing how noise is inevitably injected into a system during a measurement, noise that can be viewed as originating either from the randomness of the many degrees of freedom of the measurement device, or of the environment. This noise constitutes one component of measurement back action, and ultimately imposes limits on measurement uncertainty. Depending on the assumptions we make about active devices, and their limitations, this back action can be offset to varying degrees via control. It turns out that using active devices to reduce measurement back action leads to estimation problems that have non-zero uncertainty lower bounds, the most interesting of which arise when the observed system is lossless. One such lower bound, a main contribution of this work, can be viewed as a classical version of a Heisenberg uncertainty relation between the system's position and momentum. We finally also revisit the murky question of how macroscopic dissipation appears from lossless dynamics, and propose alternative approaches for framing the question using existing systematic methods of model reduction.

Effects of self energy of the ions on the double layer structure and properties at the dielectric interface

Although numerous theoretical efforts have been put forth, a systematic, unified and predictive theoretical framework that is able to capture all the essential physics of the interfacial behaviors of ions, such as the Hofmeister series effect, Jones-Ray effect and the salt effect on the bubble coalescence remain an outstanding challenge. The most common approach to treating electrostatic interactions in the presence of salt ions is the Poisson-Boltzmann (PB) theory. However, there are many systems for which the PB theory fails to offer even a qualitative explanation of the behavior, especially for ions distributed in the vicinity of an interface with dielectric contrast between the two media (like the water-vapor/oil interface). A key factor missing in the PB theory is the self energy of the ion.

In this thesis, we develop a self-consistent theory that treats the electrostatic self energy (including both the short-range Born solvation energy and the long-range image charge interactions), the nonelectrostatic contribution of the self energy, the ion-ion correlation and the screening effect systematically in a single framework. By assuming a finite charge spread of the ion instead of using the point-charge model, the self energy obtained by our theory is free of the divergence problems and gives a continuous self energy across the interface. This continuous feature allows ions on the water side and the vapor/oil side of the interface to be treated in a unified framework. The theory involves a minimum set of parameters of the ion, such as the valency, radius, polarizability of the ions, and the dielectric constants of the medium, that are both intrinsic and readily available. The general theory is first applied to study the thermodynamic property of the bulk electrolyte solution, which shows good agreement with the experiment result for predicting the activity coefficient and osmotic coefficient.

Next, we address the effect of local Born solvation energy on the bulk thermodynamics and interfacial properties of electrolyte solution mixtures. We show that difference in the solvation energy between the cations and anions naturally gives rise to local charge separation near the interface, and a finite Galvani potential between two coexisting solutions. The miscibility of the mixture can either increases or decreases depending on the competition between the solvation energy and translation entropy of the ions. The interfacial tension shows a non-monotonic dependence on the salt concentration: it increases linearly with the salt concentration at higher concentrations, and decreases approximately as the square root of the salt concentration for dilute solutions, which is in agreement with the Jones-Ray effect observed in experiment.

Next, we investigate the image effects on the double layer structure and interfacial properties near a single charged plate. We show that the image charge repulsion creates a depletion boundary layer that cannot be captured by a regular perturbation approach. The correct weak-coupling theory must include the self-energy of the ion due to the image charge interaction. The image force qualitatively alters the double layer structure and properties, and gives rise to many non-PB effects, such as nonmonotonic dependence of the surface energy on concentration and charge inversion. The image charge effect is then studied for electrolyte solutions between two plates. For two neutral plates, we show that depletion of the salt ions by the image charge repulsion results in short-range attractive and long-range repulsive forces. If cations and anions are of different valency, the asymmetric depletion leads to the formation of an induced electrical double layer. For two charged plates, the competition between the surface charge and the image charge effect can give rise to like- charge attraction.

Then, we study the inhomogeneous screening effect near the dielectric interface due to the anisotropic and nonuniform ion distribution. We show that the double layer structure and interfacial properties is drastically affected by the inhomogeneous screening if the bulk Debye screening length is comparable or smaller than the Bjerrum length. The width of the depletion layer is characterized by the Bjerrum length, independent of the salt concentration. We predict that the negative adsorption of ions at the interface increases linearly with the salt concentration, which cannot be captured by either the bulk screening approximation or the WKB approximation. For asymmetric salt, the inhomogeneous screening enhances the charge separation in the induced double layer and significantly increases the value of the surface potential.

Finally, to account for the ion specificity, we study the self energy of a single ion across the dielectric interface. The ion is considered to be polarizable: its charge distribution can be self-adjusted to the local dielectric environment to minimize the self energy. Using intrinsic parameters of the ions, such as the valency, radius, and polarizability, we predict the specific ion effect on the interfacial affinity of halogen anions at the water/air interface, and the strong adsorption of hydrophobic ions at the water/oil interface, in agreement with experiments and atomistic simulations.

The theory developed in this work represents the most systematic theoretical technique for weak-coupling electrolytes. We expect the theory to be more useful for studying a wide range of structural and dynamic properties in physicochemical, colloidal, soft-matter and biophysical systems.

Theory of viscous and thermal attenuation of sound by small spheres

The problem is to calculate the attenuation of plane sound waves passing through a viscous, heat-conducting fluid containing small spherical inhomogeneities. The attenuation is calculated by evaluating the rate of increase of entropy caused by two irreversible processes: (1) the mechanical work done by the viscous stresses in the presence of velocity gradients, and (2) the flow of heat down the thermal gradients. The method is first applied to a homogeneous fluid with no spheres and shown to give the classical Stokes-Kirchhoff expressions. The method is then used to calculate the additional viscous and thermal attenuation when small spheres are present. The viscous attenuation agrees with Epstein's result obtained in 1941 for a non-heat-conducting fluid. The thermal attenuation is found to be similar in form to the viscous attenuation and, for gases, of comparable magnitude. The general results are applied to the case of water drops in air and air bubbles in water.

For water drops in air the viscous and thermal attenuations are camparable; the thermal losses occur almost entirely in the air, the thermal dissipation in the water being negligible. The theoretical values are compared with Knudsen's experimental data for fogs and found to agree in order of magnitude and dependence on frequency. For air bubbles in water the viscous losses are negligible and the calculated attenuation is almost completely due to thermal losses occurring in the air inside the bubbles, the thermal dissipation in the water being relatively small. (These results apply only to non-resonant bubbles whose radius changes but slightly during the acoustic cycle.)