Laser-TIG–hybridihitsauksen soveltuvuustutkimus

Tässä tutkimuksessa selvitetään ilman hitsauslisäainetta tapahtuvan laser–TIG–hybridihitsausprosessin soveltuvuus 6 mm ja 8 mm paksujen päittäisliitettyjen S355 K2 ja Laser 355 MC rakenneterästen hitsaukseen. Hitsien tarkastelussa huomio kiinnitetään hitsausnopeuteen, hitsien tunkeumaan, liittämistehokkuuteen, hitsien kovuuteen ja hitsausliitoksen ulkonäköön. Muita tutkittavia asioita ovat laser-TIG-hybridihitsattujen levyjen muodonmuutokset ja suuresta hitsausnopeudesta sekä pienestä t8/5 jäähtymisajasta johtuvat mahdolliset kylmähalkeamat. Laser-TIG-hybridihitsejä verrataan robotti-MAG- ja käsin MAG-hitseihin sekä kaarihitsausstandardin SFS-EN ISO 5817 hitsiluokkien mukaisiin raja-arvoihin. Laser-TIG-hybridihitsausprosessissa TIG-valokaari mahdollistaa tasaisen ja lähes roiskeettoman hitsin ja lasersäde aikaansaa syvän tunkeuman sekä tasalaatuisen juurihitsin. Laser-TIG-hybridihitsausprosessilla 6 mm paksut S355 K2 rakenneteräslevyt on mahdollista hitsata levyn yhdeltä puolelta kerralla valmiiksi. Paksummat 8 mm levyt voidaan hitsata levyn yhdeltä tai molemmilta puolilta suoritettavalla laser-TIG-hybridihitsauksella. Laser-TIG-hybridihitsausprosessilla hitsatut hitsit ovat hyvin siistejä ja lähes roiskeettomia. Verrattaessa laser-TIG-hybridihitsausprosessia muihin hitsausprosesseihin sen voidaan todeta olevan erittäin kilpailukykyinen 6 mm paksujen päittäisliitettyjen rakenneterästen hitsaamisessa, mutta se soveltuu myös 8 mm paksujen rakenneterästen hitsaamiseen. Tutkitut hitsit täyttävät kaarihitsausstandardin SFS-EN ISO 5817 B- ja D-hitsiluokkien mukaiset raja-arvot. Vertailukokeet 6 mm paksulla S355 rakenneteräksellä osoittavat, että yhdeltä puolelta suoritettavan laser-TIG-hybridihitsauksen hitsausnopeus on robotti-MAG-hitsaukseen verrattuna yli nelinkertainen ja MAG-käsinhitsaukseen verrattuna yli viisinkertainen. Laser-TIG-hybridihitsauksessa liittämistehokkuus on noin viisinkertainen robotti-MAGhitsaukseen verrattuna. Molemmilta puolilta suoritettavalla laser-TIG-hybridihitsauksella voidaan 8 mm paksulla S355 rakenneteräksellä saavuttaa noin kolminkertainen hitsausnopeus ja liittämistehokkuus robotti-MAG-hitsaukseen verrattuna. Laser-TIG-hybridihitsauksessa TIG-kaaren tuoman lisälämmön ansiosta suurillakin hitsausnopeuksilla (1 m/min) voidaan saavuttaa edulliset kovuusarvot. Kovuusmittausten tulosten perusteella 6 mm ja 8 mm paksujen S355 K2 ja Laser 355 MC rakenneterästen hitsit eivät ylittäneet kaarihitsausstandardin määrittelemää 350 HV kovuuden enimmäisrajaa. Laser-TIG-hybridihitsauksen edullisesta lämmöntuonnista johtuen levyjen pituus- ja poikittaissuuntaiset muodonmuutokset ovat noin 80 prosenttia pienemmät kuin käsin suoritettavassa MAG-hitsauksessa. Laser-TIG-hybridihitsausprosessilla käytetään I-railoa, mutta robotti-MAG- ja käsin MAG-hitsausprosesseilla joudutaan käyttämään V-railoa, jolloin lämmöntuonti ja siitä johtuvat muodonmuutokset ovat suuremmat. Korkea liittämistehokkuus ja edullinen lämmöntuonti merkitsevät vähäisempiä muodonmuutoksia ja siten merkittäviä säästöjä työ-, materiaali- ja energiakustannuksissa. 8 mm ja sitä paksummilla S355 rakenneteräksillä levyn yhdeltä puolelta suoritettava päittäisliitoksen hitsaaminen on laser-TIG hybridihitsauksella haastavaa, koska yli 200 A:n TIG-kaarivirralla suuri metallisula aiheuttaa avaimenreiän sulkeutumisen ja avaimenreiän alaosaan muodostuu kaasukuplia. Tästä voidaan tehdä sellainen johtopäätös, että päittäisliitettävien levyjen ilmarakoa pitäisi kasvattaa niin suureksi, että avaimenreiän sulavirtaus ei pääse estymään. Yli 0,25 mm:n ilmarako edellyttää lasersäteen vaaputusta tai säteen halkaisijan kasvattamista. Ilmaraon kasvattaminen edellyttää myös lisäaineen käyttöä. Tutkimustulosten perusteella laser-TIG-hybridihitsausprosessilla voidaan saavuttaa merkittäviä etuja ja kustannussäästöjä, joten sen hyödyntämistä kannattaa harkita 8 mm ja sitä ohuempien päittäisliitettävien tuotteiden konepaja- ja tehdastuotannossa. Laser-TIGhybridihitsausprosessi soveltuu esimerkiksi seuraavien tuotteiden hitsaamiseen: päittäisliitettävät levyt, palkit, koneenosat, putket, säiliöt ja erilaiset pyörähdyskappaleet.

The effects of location, feedstock availability, and supply-chain logistics on the greenhouse gas emissions of forest-biomass energy utilization in Finland

Forest biomass represents a geographically distributed feedstock, and geographical location affects the greenhouse gas (GHG) performance of a given forest-bioenergy system in several ways. For example, biomass availability, forest operations, transportation possibilities and the distances involved, biomass end-use possibilities, fossil reference systems, and forest carbon balances all depend to some extent on location. The overall objective of this thesis was to assess the GHG emissions derived from supply and energy-utilization chains of forest biomass in Finland, with a specific focus on the effect of location in relation to forest biomass’s availability and the transportation possibilities. Biomass availability and transportation-network assessments were conducted through utilization of geographical information system methods, and the GHG emissions were assessed by means of lifecycle assessment. The thesis is based on four papers in which forest biomass supply on industrial scale was assessed. The feedstocks assessed in this thesis include harvesting residues, smalldiameter energy wood and stumps. The principal implication of the findings in this thesis is that in Finland, the location and availability of biomass in the proximity of a given energyutilization or energy-conversion plant is not a decisive factor in supply-chain GHG emissions or the possible GHG savings to be achieved with forest-biomass energy use. Therefore, for the greatest GHG reductions with limited forest-biomass resources, energy utilization of forest biomass in Finland should be directed to the locations where most GHG savings are achieved through replacement of fossil fuels. Furthermore, one should prioritize the types of forest biomass with the lowest direct supply-chain GHG emissions (e.g., from transport and comminution) and the lowest indirect ones (in particular, soil carbon-stock losses), regardless of location. In this respect, the best combination is to use harvesting residues in combined heat and power production, replacing peat or coal.

Mitigation possibilities to Greenhouse gas emissions from power production in India

Global warming is assertively the greatest environmental challenge for humans of 21st century. It is primarily caused by the anthropogenic greenhouse gas (GHG) that trap heat in the atmosphere. Because of which, the GHG emission mitigation, globally, is a critical issue in the political agenda of all high-profile nations. India, like other developing countries, is facing this threat of climate change while dealing with the challenge of sustaining its rapid economic growth. India’s economy is closely connected to its natural resource base and climate sensitive sectors like water, agriculture and forestry. Due to Climate change the quality and distribution of India’s natural resources may transform and lead to adverse effects on livelihood of its people. Therefore, India is expected to face a major threat due to the projected climate change. This study proposes possible solutions for GHG emission mitigation that are specific to the power sector of India. The methods discussed here will take Indian power sector from present coal dominant ideology to a system, centered with renewable energy sources. The study further proposes a future scenario for 2050, based on the present Indian government policies and global energy technologies advancements.

Power to Gas - sähköenergiasta luonnonkaasua, liiketoimintamahdollisuudet Suomessa

Uusiutuvan sähköntuotannon osuuden kasvaessa kasvaa tarve tasata sähköntuotannon ja kulutuksen vaihteluita varastoimalla sähköä. Power to Gas (PtG) - sähköenergiasta luonnonkaasua tarjoaa yhden mahdollisuuden varastoida sähköä. Sähköä käytetään veden elektrolyysiin, jossa syntynyt vety käytetään metanoinissa yhdessä hiilidioksidin kanssa muodostamaan korvaavaa luonnonkaasua. Näin syntynyttä korvaava luonnonkaasua sähköstä kutsutaan e-SNG-kaasuksi. Tässä työssä tutkitaan PtG-laitoksen investointi, käyttö- ja kunnossapitokuluja. Työssä luodaan laskentamalli, jolla lasketaan PtG-laitoksen neljälle käyttötapaukselle kannattavuuslaskelma. Käyttötapauksille lasketaan myös herkkyystarkasteluja. Kannattavuuslaskelmien perusteella päätellään PtG-laitoksen liiketoimintamahdollisuudet Suomessa. Työssä laskettujen kannattavuuslaskelmien perusteella PtG-laitoksen perustapausten liiketoimintamahdollisuudet ovat huonot. Laskettujen herkkyystarkastelujen perusteella havaittiin, että investointikulut, laitoksen ajoaika ja lisätulot hapesta ja lämmöstä ovat kannattavuuden kannalta kriittisimmät menestystekijät.

Photosynthetic responses of four tropical tree species grown under gap and understorey conditions in a semi-deciduous forest

Leaf CO2 assimilation (A) as a function of photosynthetic photon flux density (Q) or intercellular CO2 concentration (Ci) and chlorophyll fluorescence measurements were carried out on four tropical woody species growing in forest gap and understorey (Bauhinia forficata Link. and Guazuma ulmifolia Lam. as pioneers, and Hymenaea courbaril L. and Esenbeckia leiocarpa Engl. as non-pioneers). Chlorophyll fluorescence indicated similar acclimation capacities of photochemical apparatus to contrasting light environments irrespective to plant species. Maximum CO2 assimilation and quantum yield derived from A/Q curves indicated higher photosynthetic capacity in pioneer than in non-pioneer species in forest gap. However, the differences among species did not show a straightforward relation with their successional status regarding data derived from A/Q curves under understorey conditions. Both successional groups are able to sustain positive carbon balance under contrasting natural light availabilities, modifying photochemical and biochemical photosynthetic traits with similar phenotypic plasticity capacity.

Potential for Greenhouse Gas Emission Reductions by Using Biomethane as a Road Transportation Fuel

The aim of this thesis is to study whether the use of biomethane as a transportation fuel is reasonable from climate change perspective. In order to identify potentials and challenges for the reduction of greenhouse gas (GHG) emissions, this dissertation focuses on GHG emission comparisons, on feasibility studies and on the effects of various calculation methodologies. The GHG emissions calculations are carried out by using life cycle assessment (LCA) methodologies. The aim of these LCA studies is to figure out the key parameters affecting the GHG emission saving potential of biomethane production and use and to give recommendations related to methodological choices. The feasibility studies are also carried out from the life cycle perspective by dividing the biomethane production chain for various operators along the life cycle of biomethane in order to recognize economic bottlenecks. Biomethane use in the transportation sector leads to GHG emission reductions compared to fossil transportation fuels in most cases. In addition, electricity and heat production from landfill gas, biogas or biomethane leads to GHG reductions as well. Electricity production for electric vehicles is also a potential route to direct biogas or biomethane energy to transportation sector. However, various factors along the life cycle of biomethane affect the GHG reduction potentials. Furthermore, the methodological selections have significant effects on the results. From economic perspective, there are factors related to different operators along the life cycle of biomethane, which are not encouraging biomethane use in the transportation sector. To minimize the greenhouse gas emissions from the life cycle of biomethane, waste feedstock should be preferred. In addition, energy consumption, methane leakages, digestate utilization and the current use of feedstock or biogas are also key factors. To increase the use of biomethane in the transportation sector, political steering is needed to improve the feasibility for the operators. From methodological perspective, it is important to recognize the aim of the life cycle assessment study. The life cycle assessment studies can be divided into two categories: 1.) To produce average GHG information of biomethane to evaluate the acceptability of biomethane use compared to fossil transportation fuels. 2.) To produce GHG information of biomethane related to actual decision-making situations. This helps to figure out the actual GHG emission changes in cases when feedstock, biogas or biomethane are already in other use. For example directing biogas from electricity production to transportation use does not necessarily lead to additional GHG emission reductions. The use of biomethane seems to have a lot of potential for the reduction of greenhouse gas emissions as a transportation fuel. However, there are various aspects related to production processes, to the current use of feedstock or biogas and to the feasibility that have to be taken into account.

Application and development of numerical methods for the modelling of innovative gas cooled fission reactors

Innovative gas cooled reactors, such as the pebble bed reactor (PBR) and the gas cooled fast reactor (GFR) offer higher efficiency and new application areas for nuclear energy. Numerical methods were applied and developed to analyse the specific features of these reactor types with fully three dimensional calculation models. In the first part of this thesis, discrete element method (DEM) was used for a physically realistic modelling of the packing of fuel pebbles in PBR geometries and methods were developed for utilising the DEM results in subsequent reactor physics and thermal-hydraulics calculations. In the second part, the flow and heat transfer for a single gas cooled fuel rod of a GFR were investigated with computational fluid dynamics (CFD) methods. An in-house DEM implementation was validated and used for packing simulations, in which the effect of several parameters on the resulting average packing density was investigated. The restitution coefficient was found out to have the most significant effect. The results can be utilised in further work to obtain a pebble bed with a specific packing density. The packing structures of selected pebble beds were also analysed in detail and local variations in the packing density were observed, which should be taken into account especially in the reactor core thermal-hydraulic analyses. Two open source DEM codes were used to produce stochastic pebble bed configurations to add realism and improve the accuracy of criticality calculations performed with the Monte Carlo reactor physics code Serpent. Russian ASTRA criticality experiments were calculated. Pebble beds corresponding to the experimental specifications within measurement uncertainties were produced in DEM simulations and successfully exported into the subsequent reactor physics analysis. With the developed approach, two typical issues in Monte Carlo reactor physics calculations of pebble bed geometries were avoided. A novel method was developed and implemented as a MATLAB code to calculate porosities in the cells of a CFD calculation mesh constructed over a pebble bed obtained from DEM simulations. The code was further developed to distribute power and temperature data accurately between discrete based reactor physics and continuum based thermal-hydraulics models to enable coupled reactor core calculations. The developed method was also found useful for analysing sphere packings in general. CFD calculations were performed to investigate the pressure losses and heat transfer in three dimensional air cooled smooth and rib roughened rod geometries, housed inside a hexagonal flow channel representing a sub-channel of a single fuel rod of a GFR. The CFD geometry represented the test section of the L-STAR experimental facility at Karlsruhe Institute of Technology and the calculation results were compared to the corresponding experimental results. Knowledge was gained of the adequacy of various turbulence models and of the modelling requirements and issues related to the specific application. The obtained pressure loss results were in a relatively good agreement with the experimental data. Heat transfer in the smooth rod geometry was somewhat under predicted, which can partly be explained by unaccounted heat losses and uncertainties. In the rib roughened geometry heat transfer was severely under predicted by the used realisable k − epsilon turbulence model. An additional calculation with a v2 − f turbulence model showed significant improvement in the heat transfer results, which is most likely due to the better performance of the model in separated flow problems. Further investigations are suggested before using CFD to make conclusions of the heat transfer performance of rib roughened GFR fuel rod geometries. It is suggested that the viewpoints of numerical modelling are included in the planning of experiments to ease the challenging model construction and simulations and to avoid introducing additional sources of uncertainties. To facilitate the use of advanced calculation approaches, multi-physical aspects in experiments should also be considered and documented in a reasonable detail.

Working fluid selection and design of small-scale waste heat recovery systems based on organic Rankine cycles

Demand for the use of energy systems, entailing high efficiency as well as availability to harness renewable energy sources, is a key issue in order to tackling the threat of global warming and saving natural resources. Organic Rankine cycle (ORC) technology has been identified as one of the most promising technologies in recovering low-grade heat sources and in harnessing renewable energy sources that cannot be efficiently utilized by means of more conventional power systems. The ORC is based on the working principle of Rankine process, but an organic working fluid is adopted in the cycle instead of steam. This thesis presents numerical and experimental results of the study on the design of small-scale ORCs. Two main applications were selected for the thesis: waste heat re- covery from small-scale diesel engines concentrating on the utilization of the exhaust gas heat and waste heat recovery in large industrial-scale engine power plants considering the utilization of both the high and low temperature heat sources. The main objective of this work was to identify suitable working fluid candidates and to study the process and turbine design methods that can be applied when power plants based on the use of non-conventional working fluids are considered. The computational work included the use of thermodynamic analysis methods and turbine design methods that were based on the use of highly accurate fluid properties. In addition, the design and loss mechanisms in supersonic ORC turbines were studied by means of computational fluid dynamics. The results indicated that the design of ORC is highly influenced by the selection of the working fluid and cycle operational conditions. The results for the turbine designs in- dicated that the working fluid selection should not be based only on the thermodynamic analysis, but requires also considerations on the turbine design. The turbines tend to be fast rotating, entailing small blade heights at the turbine rotor inlet and highly supersonic flow in the turbine flow passages, especially when power systems with low power outputs are designed. The results indicated that the ORC is a potential solution in utilizing waste heat streams both at high and low temperatures and both in micro and larger scale appli- cations.

Thermal reactions of the major hydrocarbon components of biomass gasification gas

Gasification of biomass is an efficient method process to produce liquid fuels, heat and electricity. It is interesting especially for the Nordic countries, where raw material for the processes is readily available. The thermal reactions of light hydrocarbons are a major challenge for industrial applications. At elevated temperatures, light hydrocarbons react spontaneously to form higher molecular weight compounds. In this thesis, this phenomenon was studied by literature survey, experimental work and modeling effort. The literature survey revealed that the change in tar composition is likely caused by the kinetic entropy. The role of the surface material is deemed to be an important factor in the reactivity of the system. The experimental results were in accordance with previous publications on the subject. The novelty of the experimental work lies in the used time interval for measurements combined with an industrially relevant temperature interval. The aspects which are covered in the modeling include screening of possible numerical approaches, testing of optimization methods and kinetic modelling. No significant numerical issues were observed, so the used calculation routines are adequate for the task. Evolutionary algorithms gave a better performance combined with better fit than the conventional iterative methods such as Simplex and Levenberg-Marquardt methods. Three models were fitted on experimental data. The LLNL model was used as a reference model to which two other models were compared. A compact model which included all the observed species was developed. The parameter estimation performed on that model gave slightly impaired fit to experimental data than LLNL model, but the difference was barely significant. The third tested model concentrated on the decomposition of hydrocarbons and included a theoretical description of the formation of carbon layer on the reactor walls. The fit to experimental data was extremely good. Based on the simulation results and literature findings, it is likely that the surface coverage of carbonaceous deposits is a major factor in thermal reactions.

Potential for synthetic gas in Finnish industry

Finland, other Nordic countries and European Union aim to decarbonize their energy production by 2050. Decarbonization requires large scale implementation of non-emission energy sources, i.e. renewable energy and nuclear power. Stochastic renewable energy sources present a challenge to balance the supply and demand for energy. Energy storages, non-emissions fuels in mobility and industrial processes are required whenever electrification is not possible. Neo-Carbon project studies the decarbonizing the energy production and the role of synthetic gas in it. This thesis studies the industrial processes in steel production, oil refining, cement manufacturing and glass manufacturing, where natural gas is already used or fuel switch to SNG is possible. The technical potential for fuel switching is assessed, and economic potential is necessary after this. All studied processes have potential for fuel switching, but total decarbonization of steel production, oil refining requires implementation of other zero-emission technologies.

Tandem-MAG-welding of high strength steels for shipbuilding applications

Tämän työn tavoitteena oli hitsata tandem MAG –laitteistolla 25 mm paksua Ruukin E500 TMCP terästä. Työssä oli tarkoituksena vähentää railotilavuutta mahdollisimman paljon sekä suorittaa testihitsaukset 0.8 kJ/mm sekä 2.5 kJ/mm lämmöntuonneilla. Teoriaosuudessa käsiteltiin Tandem MAG-hitsaukseen, sen tuottavuuteen ja laatukysymyksiin liittyviä asioita sekä siinä perehdyttiin suurlujuusteräksien käyttöön hitsauksessa sekä laivanrakennuksessa. Kokeellisessa osuudessa perehdyttiin hitsauksessa huomattuihin etuihin, ongelmiin sekä ongelmien ratkaisumahdollisuuksiin. Hitsausliitoksen mekaaniset ominaisuudet tutkittiin rikkomattomin sekä rikkovin menetelmin. Alustavat hitsausohjeet luotiin kummallekin lämmöntuonnille. Testaukset aloitettiin 30 º railokulmalla pienentäen kulmaa mahdollisuuksien mukaan. Testauksissa ei saatu hitsattua onnistuneesti alle 30 º railokulmalla. Hitsaustestien aikana huomattiin magneettisen puhalluksen vaikutus hitsaustapahtumaan. Kaasunvirtausnopeuden tuli olla tietyn suuruinen jotta palkokerrokset onnistuivat ilman huokoisuusongelmaa. Pienemmällä lämmöntuonnilla hitsattaessa kaasunvirtausnopeudet olivat tärkeämpiä hitsatessa ylempiä palkokerroksia. Kääntämällä hitsauspoltinta sivuttaissuunnassa 7-10 astetta auttoi ehkäisemään reunahaavan syntymistä. Rikkovista menetelmistä testitulokset olivat hyväksyttyjä kaikkien muiden paitsi päittäishitsin sivutaivutuskokeen osalta.

Economics of Power-to-Gas integration to wastewater treatment plant

Solar and wind power produce electricity irregularly. This irregular power production is problematic and therefore production can exceed the need. Thus sufficient energy storage solutions are needed. Currently there are some storages, such as flywheel, but they are quite short-term. Power-to-Gas (P2G) offers a solution to store energy as a synthetic natural gas. It also improves nation’s energy self-sufficiency. Power-to-Gas can be integrated to an industrial or a municipal facility to reduce production costs. In this master’s thesis the integration of Power-to-Gas technologies to wastewater treatment as a part of the VTT’s Neo-Carbon Energy project is studied. Power-to-Gas produces synthetic methane (SNG) from water and carbon dioxide with electricity. This SNG can be considered as stored energy. Basic wastewater treatment technologies and the production of biogas in the treatment plant are studied. The utilisation of biogas and SNG in heat and power production and in transportation is also studied. The integration of the P2G to wastewater treatment plant (WWTP) is examined mainly from economic view. First the mass flows of flowing materials are calculated and after that the economic impact based on the mass flows. The economic efficiency is evaluated with Net Present Value method. In this thesis it is also studied the overall profitability of the integration and the key economic factors.

Dynamic Modelling of Circulating Fluidized Bed Plant with Gas Turbine Repowering

Increasing amount of renewable energy source based electricity production has set high load control requirements for power grid balance markets. The essential grid balance between electricity consumption and generation is currently hard to achieve economically with new-generation solutions. Therefore conventional combustion power generation will be examined in this thesis as a solution to the foregoing issue. Circulating fluidized bed (CFB) technology is known to have sufficient scale to acts as a large grid balancing unit. Although the load change rate of the CFB unit is known to be moderately high, supplementary repowering solution will be evaluated in this thesis for load change maximization. The repowering heat duty is delivered to the CFB feed water preheating section by smaller gas turbine (GT) unit. Consequently, steam extraction preheating may be decreased and large amount of the gas turbine exhaust heat may be utilized in the CFB process to reach maximum plant electrical efficiency. Earlier study of the repowering has focused on the efficiency improvements and retrofitting to maximize plant electrical output. This study however presents the CFB load change improvement possibilities achieved with supplementary GT heat. The repowering study is prefaced with literature and theory review for both of the processes to maximize accuracy of the research. Both dynamic and steady-state simulations accomplished with APROS simulation tool will be used to evaluate repowering effects to the CFB unit operation. Eventually, a conceptual level analysis is completed to compare repowered plant performance to the state-of-the-art CFB performance. Based on the performed simulations, considerably good improvements to the CFB process parameters are achieved with repowering. Consequently, the results show possibilities to higher ramp rate values achieved with repowered CFB technology. This enables better plant suitability to the grid balance markets.

Phonon spectra and thermodynamic properties of rare gas solids based on empirical and semi-empirical (ab initio) two-body potentials : a comparative study /

We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.

Application of [Lambda] to the fourth perturbation theory in calculating the equation of state of rare gas solids and fcc metals

We have calculated the thermodynamic properties of monatomic fcc crystals from the high temperature limit of the Helmholtz free energy. This equation of state included the static and vibrational energy components. The latter contribution was calculated to order A4 of perturbation theory, for a range of crystal volumes, in which a nearest neighbour central force model was used. We have calculated the lattice constant, the coefficient of volume expansion, the specific heat at constant volume and at constant pressure, the adiabatic and the isothermal bulk modulus, and the Gruneisen parameter, for two of the rare gas solids, Xe and Kr, and for the fcc metals Cu, Ag, Au, Al, and Pb. The LennardJones and the Morse potential were each used to represent the atomic interactions for the rare gas solids, and only the Morse potential was used for the fcc metals. The thermodynamic properties obtained from the A4 equation of state with the Lennard-Jones potential, seem to be in reasonable agreement with experiment for temperatures up to about threequarters of the melting temperature. However, for the higher temperatures, the results are less than satisfactory. For Xe and Kr, the thermodynamic properties calculated from the A2 equation of state with the Morse potential, are qualitatively similar to the A 2 results obtained with the Lennard-Jones potential, however, the properties obtained from the A4 equation of state are in good agreement with experiment, since the contribution from the A4 terms seem to be small. The lattice contribution to the thermal properties of the fcc metals was calculated from the A4 equation of state, and these results produced a slight improvement over the properties calculated from the A2 equation of state. In order to compare the calculated specific heats and bulk moduli results with experiment~ the electronic contribution to thermal properties was taken into account~ by using the free electron model. We found that the results varied significantly with the value chosen for the number of free electrons per atom.