Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TEMPERATURE DEPENDENCE”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Calculation of far wing of allowed spectra: The water continuum

A far-wing line shape theory based on the binary collision and quasistatic approximations that is applicable for both the low- and high-frequency wings of allowed vibrational-rotational lines has been developed. This theory has been applied in order to calculate the frequency and temperature dependence of the continuous absorption coefficient for frequencies up to 10,000 cm(exp -1) for pure H2O and for H2O-N2 mixtures. The calculations are made assuming an interaction potential consisting of an isotropic Lennard-Jones part and the leading long-range anisotropic part, and utilizing the measured line strengths and transition frequencies. The results compare well with existing data, both in magnitude and in temperature dependence. This leads us to the conclusion that although dimer and collision-induced absorptions are present, the primary mechanism responsible for the observed water continuum is the far-wing absorption of allowed lines. Recent progress on near-wing corrections to the theory and validations with recent laboratory measurements are discussed briefly.

Tipping, R. H.↗
DOE OSTI2024

X-ray sources for in situ wavelength calibration of x-ray imaging crystal spectrometers

X-ray sources for a range of wavelengths are being considered for in situ calibration of X-ray Imaging Crystal Spectrometers (XICSs) and for monitoring line shifts due to changes in the crystal temperature, which can vary during experimental operation over a day [A. Ince-Cushman et al., Rev. Sci. Instrum. 79, 10E302 (2008), L. Delgado-Aparicio et al., Plasma Phys. Control. Fusion 55, 125011 (2013)]. Such crystal temperature dependent shifts, if not accounted for, could be erroneously interpreted as Doppler shifts leading to errors in plasma flow-velocity measurements. The x-ray sources encompass characteristic x-ray lines falling within the wavelength range of 0.9–4.0 Å, relevant for the XICSs on present and future fusion devices. Several technological challenges associated with the development of x-ray sources for in situ calibration are identified and are being addressed in the design of multiple x-ray tubes, which will be installed inside the spectrometer housing of the XICS for the JT-60SA tokamak. These x-ray sources will be especially useful for in situ calibration between plasma discharges. Here, in this paper, laboratory experiments are described that were conducted with a Cu x-ray source, a heated quartz (102) crystal, and a pixelated Pilatus detector to measure the temperature dependent shifts of the Cu Kα 1 and Kα 2 lines at 1.5405 and 1.5443 Å, respectively, and to evaluate the 2d-lattice constant for the Bragg reflecting crystal planes as a function of temperature, which, in the case of in situ wavelength calibration, would have to be used for numerical analysis of the x-ray spectra from the plasma.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Leveraging Polymorphism in YbCuBi to Map Transport and Elastic Properties

AMX Zintl compounds with the hexagonal ZrBeSi structure have gained significant attention for their remarkable vacancy tolerance and low thermal conductivity. Their 2D honeycomb sublattice, composed of M–X covalent bonds, is believed to contribute to high anharmonicity and unusual thermal transport properties. In this study, we explore the temperature-dependent polymorphism of YbCuBi as a model system to investigate the relationship between the structure and elastic and thermal transport properties in AMX Zintls. YbCuBi undergoes a structural transition from the “flat” Cu–Bi layers in the ZrBeSi structure to corrugated layers in the LiGaGe structure below 410 K, resulting in a distortion of its centrosymmetric structure. To probe the effects of this crystallographic transition, we employ inelastic neutron scattering and temperature-dependent resonant ultrasound spectroscopy. These experimental findings, coupled with first-principles calculations and thermal conductivity measurements, allow us to elucidate a direct relationship between corrugation of the honeycomb lattice and the observed changes in elastic and thermal transport properties. These insights can be extended to other Zintl phases with similar structure types, providing a platform for the rational design of functional materials with tailored thermal properties.

Chemical structure↗
NASA NTRS1978

On the possibility of negative activation energies in bimolecular reactions

The temperature dependence of the rate constants for model reacting systems was studied to understand some recent experimental measurements which imply the existence of negative activation energies. A collision theory model and classical trajectory calculations are used to demonstrate that the reaction probability can vary inversely with collision energy for bimolecular reactions occurring on attractive potential energy surfaces. However, this is not a sufficient condition to ensure that the rate constant has a negative temperature dependence. On the basis of these calculations, it seems unlikely that a true bimolecular reaction between neutral molecules will have a negative activation energy.

Jaffe, R. L.↗
NASA NTRS1992

Activated transport in AMTEC electrodes

Transport of alkali, metal atoms through porous cathodes of alkali metal thermal-to-electric converter (AMTEC) cells is responsible for significant reducible losses in the electrical performance of these cells. Experimental evidence for activated transport of metal atoms at grain surfaces and boundaries within some AMTEC electrodes has been derived from temperature dependent studies as well as from analysis of the detailed frequency dependence of ac impedance results for other electrodes, including thin, mature molybdenum electrodes which exhibit transport dominated by free molecular flow of sodium gas at low frequencies or dc conditions. Activated surface transport will almost always exist in parallel with free molecular flow transport, and the process of alkali atom adsorption/desorption from the electrode surface will invariably be part of the transport process, and possibly a dominant part in some cases. The temperature dependence of the diffusion coefficient of the alkali metal through the electrode in several cases provides an activation energy and preexponential, but at least two activated processes may be operative, and the activation parameters should be expected to depend on the alkali metal activity gradient that the electrode experiences. In the case of Pt/W/Mn electrodes operated for 2500 hours, limiting currents varied with electrode thickness, and the activation parameters could be assigned primarily to the surface/grain boundary diffusion process.

Williams, R. M.↗
DOE OSTI2020

Transient magnetic field diffusion considerations relevant to magnetically assisted indirect drive inertial confinement fusion

We report application of a magnetic field to an indirect drive inertial confinement fusion target requires diffusion of the field through the high-Z and electrically conducting Hohlraum. The onset of the external field generates eddy currents in the Hohlraum wall that result in (1) a reduction of the peak field at the capsule, (2) heating of the Hohlraum wall through Ohmic dissipation, and (3) wall movement due to the inward force from the eddy current interacting with the field. Heating of the wall causes an increase in blackbody radiation which can preheat the capsule and frozen deuterium–tritium fuel, while wall motion leads to potential misalignment of the lasers at the Hohlraum wall. Limiting these detrimental effects sets requirements on the tolerable magnitude of each effect. We present a nonlinear model for B-field diffusion through an infinitely long thin-walled cylinder with a temperature dependent resistivity, to show that a 15 μm thick wall of pure gold fails to meet these requirements. A new Hohlraum material made from an alloy of Au and Ta has a measured resistivity of ≥60 times that of Au and is shown with the nonlinear model to meet the requirements for magnetization. We compare the nonlinear model to simulations of the actual Hohlraum target using a finite element code which includes temperature-dependent Hohlraum resistivity.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2006

Ozone Observations by the Gas and Aerosol Measurement Sensor during SOLVE II

The Gas and Aerosol Measurement Sensor (GAMS) was deployed aboard the NASA DC-8 aircraft during the second SAGE III Ozone Loss and Validation Experiment (SOLVE II). GAMS acquired line-of-sight (LOS) direct solar irradiance spectra during the sunlit portions of ten science flights of the DC-8 between 12 January and 4 February 2003. Differential line-of-sight (DLOS) optical depth spectra are produced from the GAMS raw solar irradiance spectra. Then, DLOS ozone number densities are retrieved from the GAMS spectra using a multiple linear regression spectral fitting technique. Both the DLOS optical depth spectra and retrieved ozone data are compared with coincident measurements from two other solar instruments aboard the DC-8 platform to demonstrate the robustness and stability of the GAMS data. The GAMS ozone measurements are then utilized to evaluate the quality of the Wulf band ozone cross sections, a critical component of the SAGE III aerosol, water vapor, and temperature/pressure retrievals. Results suggest the ozone cross section compilation of Shettle and Anderson currently used operationally in SAGE III data processing may be in error by as much as 10-20% in theWulf bands, and their lack of reported temperature dependence is a significant deficiency. A second, more recent, cross section database compiled for the SCIAMACHY satellite mission appears to be of much better quality in the Wulf bands, but still may have errors as large as 5% near the Wulf band absorption peaks, which is slightly larger than their stated uncertainty. Additional laboratory measurements of the Wulf band cross sections should be pursued to further reduce their uncertainty and better quantify their temperature dependence.

Pitts, M. C.↗
NASA NTRSDate not supplied

Thermophysical Property Measurements of Indium Iodide Crystals

InI single crystals are a promising room temperature detector material for X-rays and γ-rays. To improve crystal growth of the material by simulations, knowledge of thermophysical properties is essential, and since InI is orthorhombic, the anisotropy has to be taken into account. The temperature dependence of several thermophysical properties have been measured for InI, including the anisotropic thermal expansion, specific heat, and the thermal diffusivity in the b direction. The anisotropic thermal expansion coefficients, determined by X-ray diffraction, were α11=1.03·10-5 K-1, α22= 3.77·10-5 K-1, and α33 = 6.26·10-5 K-1. The specific heat, measured by DSC, was 0.226 J g-1 K-1 at 335K, with a temperature dependence of 9.582·10-5 J g-1 K-2. In the course of the X-ray diffraction and DSC measurements, it was shown that supposed phase changes, reported in older literature, are actually not phase changes but oxidation effects. The thermal diffusivity in the b direction, measured by the Xenon Flash method, was 0.288·10-6 m2s-1 at RT, decreasing to 0.253·10-6 m2 s-1 at 450K. In addition, the volume increase upon melting and the thermal expansion of the melt have been determined.

Martin Volz↗
DOE OSTI2025

Self‐Trapped Hole Migration and Defect‐Mediated Thermal Quenching of Luminescence in α‐ and β‐Ga 2 O 3

Gallium oxide (Ga 2 O 3 ) is a promising ultrawide bandgap semiconductor for next-generation power electronics and optoelectronic devices. Here, temperature-dependent and polarization-resolved photoluminescence excitation spectroscopy data, complemented by hybrid-functional first-principles calculations, are presented, and a microscopic model is derived that explains the interplay of hole migration, defect trapping, and carrier recombination at defects underlying thermal quenching phenomena in α- and β-Ga 2 O 3 . In α-Ga 2 O 3 , the UV emission is attributed to self-trapped holes, while the blue luminescence arises from defect-related processes, including gallium split vacancies and their defect complexes. Calculations reveal an energy barrier of 88 meV for self-trapped hole migration in α-Ga 2 O 3 , consistent with activation energies from temperature-dependent photoluminescence. This enables efficient trapping by defects, enhancing blue luminescence and quenching UV emission. In β-Ga 2 O 3 , a higher migration barrier of 0.36 eV reduces the defect trapping, allowing the UV self-trapped hole emission to remain intense, with blue luminescence emerging only at elevated temperatures. These results establish a direct link between self-trapped hole migration, defect trapping, and thermal quenching of emission in both phases. The insights advance the understanding of carrier dynamics in ultrawide bandgap oxides and may guide defect engineering for high-performance functional devices.

Hajizadeh, Nima [Leibniz-Institut im Forschungsver↗
NASA NTRS1994

Electrical characterization of 6H crystalline silicon carbide

Crystalline silicon carbide (SiC) substrates and epilayers, undoped as well as n- and p-doped, have been electrically characterized by performing Hall effect and resistivity measurements (van der Pauw) over the temperature range of approximately 85 K to 650 K (200 K to 500 K for p-type sample). By fitting the measured temperature dependent carrier concentration data to the single activation energy theoretical model: (1) the activation energy for the nitrogen donor ranged from 0.078 eV to 0.101 eV for a doping concentration range of 10(exp 17) cm(exp -3) to 10(exp 18) cm(exp -3) and (2) the activation energy for the aluminum acceptor was 0.252 eV for a doping concentration of 4.6 x 10(exp 18) cm(exp -3). By fitting the measured temperature dependent carrier concentration data to the double activation energy level theoretical model for the nitrogen donor: (1) the activation energy for the hexagonal site was 0.056 eV and 0.093 eV corresponding to doping concentrations of 3.33 x 10 (exp 17) cm(exp -3) and 1.6 x 10(exp 18) cm(exp -3) and (2) the activation energy for the cubic site was 0.113 and 0.126 eV corresponding to doping concentrations of 4.2 x 10(exp 17) cm(exp -3) and 5.4 x 10(exp 18) cm(exp -3).

Lempner, Stephen E.↗
DOE OSTI2021

Deep glassy state dynamic data challenge glass models: Configurational entropy models

This work was the result of the student's (Dongjie Chen) qualifying exam response. In it the experimental observation from our labs that the relaxation times of ultrastable glasses (a 20 million year old amber and a vapor deposited amorphous teflon) do not follow the expected Vogel-Fulcher or WLF types of divergence as the temperatures fall below the glass transition temperature. In the work we evaluated several theories by considering them in different ways. We used them as generally used and show that they do not capture the non-diverging behavior. Because the models were all entropy-based models, we also considered the possibility that the temperature dependence of the entropy used in the models is incorrect and chose two ways of estimating the entropy that do not follow the classical theoretical models. First we used the actually measured calorimetric response for the two ultrastable glasses, and we used the Milchev model, which does not go to zero above 0 Kelvin. Both of these approaches led to improved agreement between the theories and the experimental data. Finally, we also asked what the temperature dependence of the entropy would have to be so that the models agreed with the observed data. In some cases the entropy so-calculated was reasonable, while in others it was unphysical. In sum, the important point of the work is that the observation of non-diverging time-scales that suggest that there is no ideal glass transition can be reasonably captured by current theories if they take more realistic models of the entropy than is usually the case.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Spin-orbit exciton in a honeycomb lattice magnet CoTiO 3 : Revealing a link between magnetism in d - and f -electron systems

We carried out inelastic neutron scattering to study the spin-orbital (SO) exciton in a single crystal sample of CoTiO 3 as a function of temperature. CoTiO 3 is a honeycomb magnet with dominant XY-type magnetic interaction and an A-type antiferromagnetic order below T N ≈ 38 K. We found that the SO exciton becomes softer, but acquires a larger bandwidth in the paramagnetic phase, compared to that in the magnetically ordered phase. Moreover, an additional mode is only observed in the intermediate temperature range, as the sample is warmed up above the lowest accessible temperature below T N . Such an unusual temperature dependence observed in this material suggests that its ground states (an S eff = 1 2 doublet) and excited states multiplets are strongly coupled, and therefore cannot be treated independently, as often done in a pseudo-spin model. Furthermore, our observations can be explained by a multi-level theory within random phase approximation that explicitly takes into account both the ground and excited multiplets. The success of our theory, which is originally developed to explain temperature dependence of magnetic excitations in the rare-earth magnets, highlight the similarity between the magnetic excitations in rare-earth systems and those in transition metal systems with strong spin orbit coupling.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1997

Kinetics of the Reactions of Cl((sup 2)P(sub J)) and Br((sup 2)P(sub 3/2)) with O3

A laser flash photolysis-resonance fluorescence technique has been employed to study the kinetics of the important stratospheric reactions Cl((sup 2)P(sub J)) + O3 yields ClO + O2 and Br((sup 2)P(sub 3/2)) + O3 yields BrO + O2 as a function of temperature. The temperature dependence observed for the Cl((sup 2)P(sub J)) + O3 reaction is nonArrhenius, but can be adequately described by the following two Arrhenius expressions (units are cu cm/(molecule.s), errors are 2 sigma and represent precision only): k(sub 1)(T) = (1.19 +/- 0.21) x 10(exp -11) exp[(-33 +/- 37)/T] for T = 189-269 K and k(sub 1)(T) = (2.49 +/- 0.38) x 10(exp -11) exp[(-233 +/- 46)/T] for 269-385 K. At temperatures below 230 K, the rate coefficients determined in this study are faster than any reported previously. Incorporation of our values for k(sub 1)(T) into stratospheric models would increase calculated ClO levels and decrease calculated HCI levels; hence the calculated efficiency of ClO catalyzed ozone destruction would increase. The temperature dependence observed for the Br((sup 2)P(sub 3/2)) + O3 reaction is adequately described by the following Arrhenius expression (units are cu cm/(molecule.s), errors are 2 sigma and represent precision only): k(sub 2)(T) = (1.50 +/- 0.16) x 10(exp -11)exp[(-775 +/- 30)/T for 195-392 K. While not in quantitative agreement with Arrhenius parameters reported in most previous studies, our results almost exactly reproduce the average of all earlier studies and therefore will not affect the choice of k(sub 2)(T) for use in modeling stratospheric BrO2 chemistry.

Nicovich, J. M.↗
DOE OSTI2021

Structure-property correlations and scaling in the magnetic and magnetocaloric properties of GdCrO 3 particles

The structure, magnetic, and magnetocaloric (MC) properties of orthorhombic nanocrystalline GdCrO 3 with six particle sizes: < d > = 87, 103, 145, 224, 318, and 352 nm are reported. The particle size was tailored by annealing under different temperatures and estimated by scanning electron microscopy. With increase in < d >, Goldschmidt tolerance factor t, orthorhombic strain s, and out-of-plane Cr–O 1 –Cr bond angle first decrease, reaching minimum values for < d > = 224 nm, and then increase for sample with < d > = 318 nm and 352 nm, thus showing a V-shaped variation. Temperature dependence of the magnetization (M) reveals an antiferromagnetic transition at ${T}_{\text{N}}^{\text{Cr}}\sim 168$ K for < d > ≥ 224 nm and ${T}_{\text{N}}^{\text{Cr}}\sim 167$ K for < d > < 224 nm and an essentially d-independent spin-reorientation at T SR = 9 K. M measured at 5 K and 7 T first increases with increase in < d >, reaching maximum value for sample with < d > = 224 nm, and then decreases for samples with < d > = 318 nm and 352 nm, showing an inverted-V variation with < d >. Similar < d >-dependence is observed for the magnetic entropy change (MEC) and relative cooling power (RCP) showing a close relationship between the structural and magnetic properties of GdCrO 3 nanoparticles investigated here. The 224 nm sample with the minimum values of t, s, and Cr–O 1 –Cr bond angle exhibits the maximum value of MEC (–ΔS) = 37.8 J kg –1 K –1 at 5 K under a field variation (ΔH) of 7 T and its large estimated RCP of 623.6 J Kg –1 is comparable with those of typical MC materials. Both (–ΔS) and RCP are shown to scale with the saturation magnetization MS, suggesting that MS is the crucial factor controlling their magnitudes. Assuming (–ΔS) ~ (ΔH) n , the temperature dependence of n for the six samples are determined, n varying between 1.3 at 5 K to n = 2.2 at 130 K in line with its expected magnitudes based on mean-field theory. Furthermore, these results on structure-property correlations and scaling in GdCrO 3 suggest that its MC properties are tunable for potential low-temperature magnetic refrigeration applications.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Microscopic characterization of the magnetic properties of the itinerant antiferromagnet La 2 Ni 7 by La 139 NMR/NQR measurements

In this study, 139 La nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) measurements have been performed to investigate the magnetic properties of the itinerant magnet La 2 Ni 7 , which shows a series of antiferromagnetic (AFM) phase transitions at T N1 = 61 K, T N2 = 56 K, and T N3 = 42 K under zero magnetic field. Two distinct La NMR signals were observed due to the two crystallographically inequivalent La sites in La 2 Ni 7 (La1 and La2 in the La 2 Ni 4 and the LaNi 5 subunits of the La 2 Ni 7 unit cell, respectively). From the 139 La NQR spectrum in the AFM state below T N3 , the AFM state was revealed to be a commensurate state where Ni-ordered moments align along the crystalline c axis. Owing to the two different La sites, we were able to estimate the average values of the Ni-ordered moments (0.09 – 0.10 μB / Ni and 0.17 μB / Ni around La1 and La2, respectively) from 139 La NMR spectrum measurements in the AFM state below T N3 , suggesting a nonuniform distribution of the Ni-ordered moments in the AFM state. In contrast, a more uniform distribution of the Ni-ordered moments in the saturated paramagnetic state induced by the application of high magnetic fields is observed. The temperature dependence of the sublattice magnetization measured by the internal magnetic induction at the La2 site in the AFM state was reproduced by a local-moment model better than the self-consistent renormalization (SCR) theory for weak itinerant antiferromagnets. Given the small Ni-ordered moments in the magnetically ordered state, our results suggest that La 2 Ni 7 has characteristics of both itinerant nature and localized nature in its magnetism. With this in mind, it is noteworthy that the temperature dependence of nuclear spin-relaxation rates (1/T 1 ) in the paramagnetic state above T N1 measured at zero magnetic field can be explained qualitatively by both the SCR theory and the local-moment model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2016

Multi-Dimensional, Non-Pyrolyzing Ablation Test Problems

Non-pyrolyzingcarbonaceous materials represent a class of candidate material for hypersonic vehicle components providing both structural and thermal protection system capabilities. Two problems relevant to this technology are presented. The first considers the one-dimensional ablation of a carbon material subject to convective heating. The second considers two-dimensional conduction in a rectangular block subject to radiative heating. Surface thermochemistry for both problems includes finite-rate surface kinetics at low temperatures, diffusion limited ablation at intermediate temperatures, and vaporization at high temperatures. The first problem requires the solution of both the steady-state thermal profile with respect to the ablating surface and the transient thermal history for a one-dimensional ablating planar slab with temperature-dependent material properties. The slab front face is convectively heated and also reradiates to a room temperature environment. The back face is adiabatic. The steady-state temperature profile and steady-state mass loss rate should be predicted. Time-dependent front and back face temperature, surface recession and recession rate along with the final temperature profile should be predicted for the time-dependent solution. The second problem requires the solution for the transient temperature history for an ablating, two-dimensional rectangular solid with anisotropic, temperature-dependent thermal properties. The front face is radiatively heated, convectively cooled, and also reradiates to a room temperature environment. The back face and sidewalls are adiabatic. The solution should include the following 9 items: final surface recession profile, time-dependent temperature history of both the front face and back face at both the centerline and sidewall, as well as the time-dependent surface recession and recession rate on the front face at both the centerline and sidewall. The results of the problems from all submitters will be collected, summarized, and presented at a later conference.

Ablation↗
DOE OSTI2020

Vacancies and spin-phonon coupling in CrSi 0.8 Ge 0.1 Te 3

We report temperature-dependent Raman scattering and magnetization studies of van der Waals ferromagnetic compound CrSi 0.8 Ge 0.1 Te 3 . Magnetic susceptibility measurements revealed dominant ferromagnetic interactions below T C which shift to the lower values due to the presence of vacancies. A Raman active mode, additional to the ones predicted by symmetry in the parent compounds, has been observed. This A g symmetry mode most likely emerges as a consequence of the atomic vacancies on Si/Ge site. Presence of the strong spin–phonon coupling at temperature around 210 K is indicated by deviations from conventional phonon self-energy temperature dependence of all analysed modes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Toward Quantum Chemical Free Energy Simulations of Platinum Nanoparticles on Titania Support

Platinum nanoparticles (Pt-NPs) supported on titania surfaces are costly but indispensable heterogeneous catalysts because of their highly effective and selective catalytic properties. Therefore, it is vital to understand their physicochemical processes during catalysis to optimize their use and to further develop better catalysts. However, simulating these dynamic processes is challenging due to the need for a reliable quantum chemical method to describe chemical bond breaking and bond formation during the processes but, at the same time, fast enough to sample a large number of configurations required to compute the corresponding free energy surfaces. Density functional theory (DFT) is often used to explore Pt-NPs; nonetheless, it is usually limited to some minimum-energy reaction pathways on static potential energy surfaces because of its high computational cost. In this work, we report a combination of the density functional tight binding (DFTB) method as a fast but reliable approximation to DFT, the steered molecular dynamics (SMD) technique, and the Jarzynski equality to construct free energy surfaces of the temperature-dependent diffusion and growth of platinum particles on a titania surface. In particular, we present the parametrization for Pt-X (X = Pt, Ti, or O) interactions in the framework of the second-order DFTB method, using a previous parametrization for titania as a basis. The optimized parameter set was used to simulate the surface diffusion of a single platinum atom (Pt 1 ) and the growth of Pt 6 from Pt 5 and Pt 1 on the rutile (110) surface at three different temperatures (T = 400, 600, 800 K). The free energy profile was constructed by using over a hundred SMD trajectories for each process. We found that increasing the temperature has a minimal effect on the formation free energy; nevertheless, it significantly reduces the free energy barrier of Pt atom migration on the TiO 2 surface and the transition state (TS) of its deposition. In a concluding remark, the methodology opens the pathway to quantum chemical free energy simulations of Pt-NPs’ temperature-dependent growth and other transformation processes on the titania support.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗