The surface- and size-dependent elastic moduli of nanostructures

A theoretical model is presented to investigate the size-dependent elastic moduli of nanostructures with the effects of the surface relaxation surface energy taken into consideration. At nanoscale, due to the large ratios of the surface-to-volume, the surface effects, which include surface relaxation surface energy, etc., can play important roles. Thus, the elastic moduli of nanostructures become surface- and size-dependent. In the research, the three-dimensional continuum model of the nanofilm with the surface effects is investigated. The analytical expressions of five nonzero elastic moduli of the nanofilm are derived, and then the dependence of the elastic moduli is discussed on the surface effects and the characteristic dimensions of nanofilms.

Mechanisms Responsible for the Trophic Effect of Beta-Adrenoceptors on the I-to Current Density in Type 1 Diabetic Rat Cardiomyocytes

Background/Aims: In diabetic ventricular myocytes, transient outward potassium current (I-to) amplitude is severely reduced because of the impaired catecholamine release that characterizes diabetic autonomic neuropathy. Sympathetic nervous system exhibits a trophic effect on I-to since incubation of myocytes with noradrenaline restores current amplitude via beta-adrenoceptor (beta AR) stimulation. Here, we investigate the intracellular signalling pathway though which incubation of diabetic cardiomyocytes with the beta AR agonist isoproterenol recovers I-to amplitude to normal values. Methods: Experiments were performed in ventricular myocytes isolated from streptozotocin-diabetic rats. I-to current was recorded by using the patch-clamp technique. Kv4 channel expression was determined by immunofluorescence. Protein-protein interaction was determined by coimmunoprecipitation. Results: Stimulation of beta AR activates first a G alpha s protein, adenylyl cyclase and Protein Kinase A. PKA-phosphorylated receptor then switches to the G alpha i protein. This leads to the activation of the beta AR-Kinase-1 and further receptor phosphorylation and arrestin dependent internalization. The internalized receptor-arrestin complex recruits and activates cSrc and the MAPK cascade, where Ras, c-Raf1 and finally ERK1/2 mediate the increase in Kv4.2 and Kv4.3 protein abundance in the plasma membrane. Conclusion: beta(2)AR stimulation activates a G alpha s and G alpha i protein dependent pathway where the ERK1/2 modulates the Ito current amplitude and the density of the Kv4.2 and Kv4.2 channels in the plasma membrane upon sympathetic stimulation in diabetic heart.

Critical current anisotropy in high temperature superconductors

After nearly 15 years of research effort, High Temperature Superconductors (HTS) are finding a wide range of practical applications. A clear understanding of the factors controlling the current carrying capacity of these materials is a prerequisite to their successful technological development. The critical current density (Jc) in HTS is directly dependent on the structure and pinning of the Flux Line Lattice (FLL) in these materials. This thesis presents an investigation of the Jc anisotropy in HTS. The use of thin films grown on off c-axis (vicinal) substrates allowed the effect of current directions outside the cuprate planes to be studied. With this experimental geometry Berghuis, et al. (Phys. Rev. Lett. 79, 12, pg. 2332) observed a striking flux channelling effect in vicinal YBa2Cu3O7-δ (YBCO) films. By confirming, and extending, this observation, it is demonstrated that this is an intrinsic effect. The results obtained, appear to fit well with the predictions of a field angle dependent cross-over from a three dimensional rectilinear FLL to a kinked lattice of strings and pancakes. The pinning force density for movement of strings inside the cuprate planes is considerably less than that on vortex pancake elements. When the FLL is entirely string-like this reduced pinning leads to the observed channelling minima. It is observed that anti-phase boundaries enhance the Jc in vicinal YBCO films by strongly pinning vortex strings. The effect on the FLL structure cross-over of increasing anisotropy has been elucidated using de-oxygenated vicinal YBCO films. Intriguingly, the counter intuitive prediction that the range of applied field angle for which the kinked lattice is fully developed reduces with increasing anisotropy, appears to be confirmed. Although vortex channelling cannot be observed in c-axis YBCO films, the pinning force density for vortex string channelling has been extracted by observing string dragging. By studying the effect of rotating the applied field at a constant angle to the cuprate planes, it is possible to observe the cross-over into the string pancake regime in c-axis films. In the 3D region, the observed behaviour is well explained by the anisotropic Ginzburg-Landau model. Measurements were also made on thin films of the much more anisotropic Bi 2Sr2CaCu2O8+x material, grown on vicinal substrates. The absence of any flux channelling effect and clear adherence to the expected Kes-Law behaviour in the observed Jc characteristics does not provide evidence for the existence of the predicted ‘crossing lattice’ in Bi 2Sr2CaCu2O8+x .

Resveratrol prevents inflammation-dependent hepatic melanoma metastasis by inhibiting the secretion and effects of interleukin-18

Background: Implantation and growth of metastatic cancer cells at distant organs is promoted by inflammation-dependent mechanisms. A hepatic melanoma metastasis model where a majority of metastases are generated via interleukin-18-dependent mechanisms was used to test whether anti-inflammatory properties of resveratrol can interfere with mechanisms of metastasis. Methods: Two experimental treatment schedules were used: 1) Mice received one daily oral dose of 1 mg/kg resveratrol after cancer cell injection and the metastasis number and volume were determined on day 12. 2) Mice received one daily oral dose of 1 mg/kg resveratrol along the 5 days prior to the injection of cancer cells and both interleukin-18 (IL-18) concentration in the hepatic blood and microvascular retention of luciferase-transfected B16M cells were determined on the 18(th) hour. In vitro, primary cultured hepatic sinusoidal endothelial cells were treated with B16M-conditioned medium to mimic their in vivo activation by tumor-derived factors and the effect of resveratrol on IL-18 secretion, on vascular cell adhesion molecule-1 (VCAM-1) expression and on tumor cell adhesion were studied. The effect of resveratrol on melanoma cell activation by IL-18 was also studied. Results: Resveratrol remarkably inhibited hepatic retention and metastatic growth of melanoma cells by 50% and 75%, respectively. The mechanism involved IL-18 blockade at three levels: First, resveratrol prevented IL-18 augmentation in the blood of melanoma cell-infiltrated livers. Second, resveratrol inhibited IL-18-dependent expression of VCAM-1 by tumor-activated hepatic sinusoidal endothelium, preventing melanoma cell adhesion to the microvasculature. Third, resveratrol inhibited adhesion-and proliferation-stimulating effects of IL-18 on metastatic melanoma cells through hydrogen peroxide-dependent nuclear factor-kappaB translocation blockade on these cells. Conclusions: These results demonstrate multiple sites for therapeutic intervention using resveratrol within the prometastatic microenvironment generated by tumor-induced hepatic IL-18, and suggest a remarkable effect of resveratrol in the prevention of inflammation-dependent melanoma metastasis in the liver.

Scale-similarity model for Lagrangian velocity correlations in isotropic and stationary turbulence

A scale-similarity model for Lagrangian two-point, two-time velocity correlations LVCs in isotropic turbulence is developed from the Kolmogorov similarity hypothesis. It is a second approximation to the isocontours of LVCs, while the Smith-Hay model is only a first approximation. This model expresses the LVC by its space correlation and a dispersion velocity. We derive the analytical expression for the dispersion velocity from the Navier-Stokes equations using the quasinormality assumption. The dispersion velocity is dependent on enstrophy spectra and shown to be smaller than the sweeping velocity for the Eulerian velocity correlation. Therefore, the Lagrangian decorrelation process is slower than the Eulerian decorrelation process. The data from direct numerical simulation of isotropic turbulence support the scale-similarity model: the LVCs for different space separations collapse into a universal form when plotted against the separation axis defined by the model.

Microscopic Model Of Nano-Scale Particles Removal In High Pressure Co2-Based Solvents

A physical model is presented to describe the kinds of static forces responsible for adhesion of nano-scale copper metal particles to silicon surface with a fluid layer. To demonstrate the extent of particle cleaning, Received in revised form equilibrium separation distance (ESD) and net adhesion force (NAF) of a regulated metal particle with different radii (10-300 nm) on the silicon surface in CO2-based cleaning systems under different pressures were simulated. Generally, increasing the pressure of the cleaning system decreased the net adhesion force between spherical copper particle and silicon surface entrapped with medium. For CO2 + isopropanol cleaning system, the equilibrium separation distance exhibited a maximum at temperature 313.15 K in the Equilibrium separation distance regions of pressure space (1.84-8.02 MPa). When the dimension of copper particle was given, for example, High pressure 50 nm radius particles, the net adhesion force decreased and equilibrium separation distance increased with increased pressure in the CO2 + H2O cleaning system at temperature 348.15 K under 2.50-12.67 MPa pressure range. However, the net adhesion force and equilibrium separation distance both decreased with an increase in surfactant concentration at given pressure (27.6 or 27.5 MPa) and temperature (318 or 298 K) for CO2 + H2O with surfactant PFPE COO-NH4+ or DiF(8)-PO4-Na+. (C) 2008 Elsevier B.V. All rights reserved.

Asymptotic scaling in the two-dimensional O(3) Nonlinear sigma model: a Monte Carlo study on parallel computers

We investigate the 2d O(3) model with the standard action by Monte Carlo simulation at couplings β up to 2.05. We measure the energy density, mass gap and susceptibility of the model, and gather high statistics on lattices of size L ≤ 1024 using the Floating Point Systems T-series vector hypercube and the Thinking Machines Corp.'s Connection Machine 2. Asymptotic scaling does not appear to set in for this action, even at β = 2.10, where the correlation length is 420. We observe a 20% difference between our estimate m/Λ^─_(Ms) = 3.52(6) at this β and the recent exact analytical result . We use the overrelaxation algorithm interleaved with Metropolis updates and show that decorrelation time scales with the correlation length and the number of overrelaxation steps per sweep. We determine its effective dynamical critical exponent to be z' = 1.079(10); thus critical slowing down is reduced significantly for this local algorithm that is vectorizable and parallelizable.

We also use the cluster Monte Carlo algorithms, which are non-local Monte Carlo update schemes which can greatly increase the efficiency of computer simulations of spin models. The major computational task in these algorithms is connected component labeling, to identify clusters of connected sites on a lattice. We have devised some new SIMD component labeling algorithms, and implemented them on the Connection Machine. We investigate their performance when applied to the cluster update of the two dimensional Ising spin model.

Finally we use a Monte Carlo Renormalization Group method to directly measure the couplings of block Hamiltonians at different blocking levels. For the usual averaging block transformation we confirm the renormalized trajectory (RT) observed by Okawa. For another improved probabilistic block transformation we find the RT, showing that it is much closer to the Standard Action. We then use this block transformation to obtain the discrete β-function of the model which we compare to the perturbative result. We do not see convergence, except when using a rescaled coupling β_E to effectively resum the series. For the latter case we see agreement for m/ Λ^─_(Ms) at , β = 2.14, 2.26, 2.38 and 2.50. To three loops m/Λ^─_(Ms) = 3.047(35) at β = 2.50, which is very close to the exact value m/ Λ^─_(Ms) = 2.943. Our last point at β = 2.62 disagrees with this estimate however.

Three-dimensional models of mantle convection : influence of phase transitions and temperature-dependent viscosity

Two of the most important questions in mantle dynamics are investigated in three separate studies: the influence of phase transitions (studies 1 and 2), and the influence of temperature-dependent viscosity (study 3).

(1) Numerical modeling of mantle convection in a three-dimensional spherical shell incorporating the two major mantle phase transitions reveals an inherently three-dimensional flow pattern characterized by accumulation of cold downwellings above the 670 km discontinuity, and cylindrical 'avalanches' of upper mantle material into the lower mantle. The exothermic phase transition at 400 km depth reduces the degree of layering. A region of strongly-depressed temperature occurs at the base of the mantle. The temperature field is strongly modulated by this partial layering, both locally and in globally-averaged diagnostics. Flow penetration is strongly wavelength-dependent, with easy penetration at long wavelengths but strong inhibition at short wavelengths. The amplitude of the geoid is not significantly affected.

(2) Using a simple criterion for the deflection of an upwelling or downwelling by an endothermic phase transition, the scaling of the critical phase buoyancy parameter with the important lengthscales is obtained. The derived trends match those observed in numerical simulations, i.e., deflection is enhanced by (a) shorter wavelengths, (b) narrower up/downwellings (c) internal heating and (d) narrower phase loops.

(3) A systematic investigation into the effects of temperature-dependent viscosity on mantle convection has been performed in three-dimensional Cartesian geometry, with a factor of 1000-2500 viscosity variation, and Rayleigh numbers of 10^5-10^7. Enormous differences in model behavior are found, depending on the details of rheology, heating mode, compressibility and boundary conditions. Stress-free boundaries, compressibility, and temperature-dependent viscosity all favor long-wavelength flows, even in internally heated cases. However, small cells are obtained with some parameter combinations. Downwelling plumes and upwelling sheets are possible when viscosity is dependent solely on temperature. Viscous dissipation becomes important with temperature-dependent viscosity.

The sensitivity of mantle flow and structure to these various complexities illustrates the importance of performing mantle convection calculations with rheological and thermodynamic properties matching as closely as possible those of the Earth.

State-dependent modulation of neuronal circuits in C. elegans sleep

C. elegans is a compact system of 302 neurons with identifiable and mapped connections that makes it ideal for systems analysis. This work is a demonstration of what I have been able to learn about the nature of state-specific modulation and reversibility during a state called lethargus, a sleep-like state in the worm. I begin with description about the nervous system of the worm, the nature of sleep in the worm, the questions about behavior and its apparent circuit properties, the tools available and used to manipulate the nervous system, and what I have been able to learn from these studies. I end with clues that the physiology helps to teach us about the dynamics of state specific modulation, what makes sleep so different from other states, and how we can use these measurements to understand which modulators, neurotransmitters, and channels can be used to create different dynamics in a simple model system.

The deep ocean density structure at the Last Glacial Maximum : What was it and why?

The search for reliable proxies of past deep ocean temperature and salinity has proved difficult, thereby limiting our ability to understand the coupling of ocean circulation and climate over glacial-interglacial timescales. Previous inferences of deep ocean temperature and salinity from sediment pore fluid oxygen isotopes and chlorinity indicate that the deep ocean density structure at the Last Glacial Maximum (LGM, approximately 20,000 years BP) was set by salinity, and that the density contrast between northern and southern sourced deep waters was markedly greater than in the modern ocean. High density stratification could help explain the marked contrast in carbon isotope distribution recorded in the LGM ocean relative to that we observe today, but what made the ocean's density structure so different at the LGM? How did it evolve from one state to another? Further, given the sparsity of the LGM temperature and salinity data set, what else can we learn by increasing the spatial density of proxy records?

We investigate the cause and feasibility of a highly and salinity stratified deep ocean at the LGM and we work to increase the amount of information we can glean about the past ocean from pore fluid profiles of oxygen isotopes and chloride. Using a coupled ocean--sea ice--ice shelf cavity model we test whether the deep ocean density structure at the LGM can be explained by ice--ocean interactions over the Antarctic continental shelves, and show that a large contribution of the LGM salinity stratification can be explained through lower ocean temperature. In order to extract the maximum information from pore fluid profiles of oxygen isotopes and chloride we evaluate several inverse methods for ill-posed problems and their ability to recover bottom water histories from sediment pore fluid profiles. We demonstrate that Bayesian Markov Chain Monte Carlo parameter estimation techniques enable us to robustly recover the full solution space of bottom water histories, not only at the LGM, but through the most recent deglaciation and the Holocene up to the present. Finally, we evaluate a non-destructive pore fluid sampling technique, Rhizon samplers, in comparison to traditional squeezing methods and show that despite their promise, Rhizons are unlikely to be a good sampling tool for pore fluid measurements of oxygen isotopes and chloride.

Nonstationary response of structures and its application to earthquake engineering

This thesis presents a simplified state-variable method to solve for the nonstationary response of linear MDOF systems subjected to a modulated stationary excitation in both time and frequency domains. The resulting covariance matrix and evolutionary spectral density matrix of the response may be expressed as a product of a constant system matrix and a time-dependent matrix, the latter can be explicitly evaluated for most envelopes currently prevailing in engineering. The stationary correlation matrix of the response may be found by taking the limit of the covariance response when a unit step envelope is used. The reliability analysis can then be performed based on the first two moments of the response obtained.

The method presented facilitates obtaining explicit solutions for general linear MDOF systems and is flexible enough to be applied to different stochastic models of excitation such as the stationary models, modulated stationary models, filtered stationary models, and filtered modulated stationary models and their stochastic equivalents including the random pulse train model, filtered shot noise, and some ARMA models in earthquake engineering. This approach may also be readily incorporated into finite element codes for random vibration analysis of linear structures.

A set of explicit solutions for the response of simple linear structures subjected to modulated white noise earthquake models with four different envelopes are presented as illustration. In addition, the method has been applied to three selected topics of interest in earthquake engineering, namely, nonstationary analysis of primary-secondary systems with classical or nonclassical dampings, soil layer response and related structural reliability analysis, and the effect of the vertical components on seismic performance of structures. For all the three cases, explicit solutions are obtained, dynamic characteristics of structures are investigated, and some suggestions are given for aseismic design of structures.

Theory of ozone isotopic effects and various electron transfer reactions

Three separate topics, each stimulated by experiments, are treated theoretically in this dessertation: isotopic effects of ozone, electron transfer at interfaces, and intramolecular directional electron transfer in a supramolecular system.

The strange mass-independent isotope effect for the enrichment of ozone, which has been a puzzle in the literature for some 20 years, and the equally puzzling unconventional strong mass-dependent effect of individual reaction rate constants are studied as different aspects of a symmetry-driven behavior. A statistical (RRKM-based) theory with a hindered-rotor transition state is used. The individual rate constant ratios of recombination reactions at low pressures are calculated using the theory involving (1) small deviation from the statistical density of states for symmetric isotopomers, and (2) weak collisions for deactivation of the vibrationally excited ozone molecules. The weak collision and partitioning among exit channels play major roles in producing the large unconventional isotope effect in "unscrambled" systems. The enrichment studies reflect instead the non-statistical effect in "scrambled" systems. The theoretical results of low-pressure ozone enrichments and individual rate constant ratios obtained from these calculations are consistent with the corresponding experimental results. The isotopic exchange rate constant for the reaction ^(16)O + ^(18)O ^(18)O→+ ^(16)O ^(18)O + ^(18)O provides information on the nature of a variationally determined hindered-rotor transition state using experimental data at 130 K and 300 K. Pressure effects on the recombination rate constant, on the individual rate constant ratios and on the enrichments are also investigated. The theoretical results are consistent with the experimental data. The temperature dependence of the enrichment and rate constant ratios is also discussed, and experimental tests are suggested. The desirability of a more accurate potential energy surface for ozone in the transition state region is also noted.

Electron transfer reactions at semiconductor /liquid interfaces are studied using a tight-binding model for the semiconductors. The slab method and a z-transform method are employed in obtaining the tight-binding electronic structures of semiconductors having surfaces. The maximum electron transfer rate constants at Si/viologen^(2-/+) and InP /Me_(2)Fc^(+/O) interfaces are computed using the tight-binding type calculations for the solid and the extended-Huckel for the coupling to the redox agent at the interface. These electron transfer reactions are also studied using a free electron model for the semiconductor and the redox molecule, where Bardeen's method is adapted to calculate the coupling matrix element between the molecular and semiconductor electronic states. The calculated results for maximum rate constant of the electron transfer from the semiconductor bulk states are compared with the experimentally measured values of Lewis and coworkers, and are in reasonable agreement, without adjusting parameters. In the case of InP /liquid interface, the unusual current vs applied potential behavior is additionally interpreted, in part, by the presence of surface states.

Photoinduced electron transfer reactions in small supramolecular systems, such as 4-aminonaphthalimide compounds, are interesting in that there are, in principle, two alternative pathways (directions) for the electron transfer. The electron transfer, however, is unidirectional, as deduced from pH-dependent fluorescence quenching studies on different compounds. The role of electronic coupling matrix element and the charges in protonation are considered to explain the directionality of the electron transfer and other various results. A related mechanism is proposed to interpret the fluorescence behavior of similar molecules as fluorescent sensors of metal ions.

Intense laser beam guiding in self-induced electron cavitation channel in underdense plasmas

In underdense plasmas, the transverse ponderomotive force of an intense laser beam with Gaussian transverse profile expels electrons radially, and it can lead to an electron cavitation. An improved cavitation model with charge conservation constraint is applied to the determination of the width of the electron cavity. The envelope equation for laser spot size derived by using source-dependent expansion method is extended to including the electron cavity. The condition for self-guiding is given and illuminated by an effective potential for the laser spot size. The effects of the laser power, plasma density and energy dissipation on the self-guiding condition are discussed.

Magnetohydrodynamic turbulence

We simulate incompressible, MHD turbulence using a pseudo-spectral code. Our major conclusions are as follows.

1) MHD turbulence is most conveniently described in terms of counter propagating shear Alfvén and slow waves. Shear Alfvén waves control the cascade dynamics. Slow waves play a passive role and adopt the spectrum set by the shear Alfvén waves. Cascades composed entirely of shear Alfvén waves do not generate a significant measure of slow waves.

2) MHD turbulence is anisotropic with energy cascading more rapidly along k_⊥ than along k_∥, where k_⊥ and k_∥ refer to wavevector components perpendicular and parallel to the local magnetic field. Anisotropy increases with increasing k_⊥ such that excited modes are confined inside a cone bounded by k_∥ ∝ k^γ_⊥ where γ less than 1. The opening angle of the cone, θ(k_⊥) ∝ k^-(1-γ)_⊥, defines the scale dependent anisotropy.

3) MHD turbulence is generically strong in the sense that the waves which comprise it suffer order unity distortions on timescales comparable to their periods. Nevertheless, turbulent fluctuations are small deep inside the inertial range. Their energy density is less than that of the background field by a factor θ² (k_⊥)≪1.

4) MHD cascades are best understood geometrically. Wave packets suffer distortions as they move along magnetic field lines perturbed by counter propagating waves. Field lines perturbed by unidirectional waves map planes perpendicular to the local field into each other. Shear Alfvén waves are responsible for the mapping's shear and slow waves for its dilatation. The amplitude of the former exceeds that of the latter by 1/θ(k_⊥) which accounts for dominance of the shear Alfvén waves in controlling the cascade dynamics.

5) Passive scalars mixed by MHD turbulence adopt the same power spectrum as the velocity and magnetic field perturbations.

6) Decaying MHD turbulence is unstable to an increase of the imbalance between the flux of waves propagating in opposite directions along the magnetic field. Forced MHD turbulence displays order unity fluctuations with respect to the balanced state if excited at low k by δ(t) correlated forcing. It appears to be statistically stable to the unlimited growth of imbalance.

7) Gradients of the dynamic variables are focused into sheets aligned with the magnetic field whose thickness is comparable to the dissipation scale. Sheets formed by oppositely directed waves are uncorrelated. We suspect that these are vortex sheets which the mean magnetic field prevents from rolling up.

8) Items (1)-(5) lend support to the model of strong MHD turbulence put forth by Goldreich and Sridhar (1995, 1997). Results from our simulations are also consistent with the GS prediction γ = 2/3. The sole not able discrepancy is that the 1D power law spectra, E(k_⊥) ∝ k^-∝_⊥, determined from our simulations exhibit ∝ ≈ 3/2, whereas the GS model predicts ∝ = 5/3.

Advanced density functional theory methods for materials science

In this work we chiefly deal with two broad classes of problems in computational materials science, determining the doping mechanism in a semiconductor and developing an extreme condition equation of state. While solving certain aspects of these questions is well-trodden ground, both require extending the reach of existing methods to fully answer them. Here we choose to build upon the framework of density functional theory (DFT) which provides an efficient means to investigate a system from a quantum mechanics description.

Zinc Phosphide (Zn₃P₂) could be the basis for cheap and highly efficient solar cells. Its use in this regard is limited by the difficulty in n-type doping the material. In an effort to understand the mechanism behind this, the energetics and electronic structure of intrinsic point defects in zinc phosphide are studied using generalized Kohn-Sham theory and utilizing the Heyd, Scuseria, and Ernzerhof (HSE) hybrid functional for exchange and correlation. Novel 'perturbation extrapolation' is utilized to extend the use of the computationally expensive HSE functional to this large-scale defect system. According to calculations, the formation energy of charged phosphorus interstitial defects are very low in n-type Zn₃P₂ and act as 'electron sinks', nullifying the desired doping and lowering the fermi-level back towards the p-type regime. Going forward, this insight provides clues to fabricating useful zinc phosphide based devices. In addition, the methodology developed for this work can be applied to further doping studies in other systems.

Accurate determination of high pressure and temperature equations of state is fundamental in a variety of fields. However, it is often very difficult to cover a wide range of temperatures and pressures in an laboratory setting. Here we develop methods to determine a multi-phase equation of state for Ta through computation. The typical means of investigating thermodynamic properties is via ’classical’ molecular dynamics where the atomic motion is calculated from Newtonian mechanics with the electronic effects abstracted away into an interatomic potential function. For our purposes, a ’first principles’ approach such as DFT is useful as a classical potential is typically valid for only a portion of the phase diagram (i.e. whatever part it has been fit to). Furthermore, for extremes of temperature and pressure quantum effects become critical to accurately capture an equation of state and are very hard to capture in even complex model potentials. This requires extending the inherently zero temperature DFT to predict the finite temperature response of the system. Statistical modelling and thermodynamic integration is used to extend our results over all phases, as well as phase-coexistence regions which are at the limits of typical DFT validity. We deliver the most comprehensive and accurate equation of state that has been done for Ta. This work also lends insights that can be applied to further equation of state work in many other materials.