Tag: Materials Science

Spring School Computational Tools: Day 1

Today our one-week spring school on computational tools for materials science kicked off. During this week, Kurt Lejaeghere and I host this spring school, which we have been busily organizing the last few months, intended to introduce materials scientists into the use of four major ab-initio codes (VASP, ABINIT, Gaussian and CP2K). During this first day, all participants are immersed in the theoretical background of molecular modeling and solid state physics.

springschool

Prof. Karen Hemelsoet presented a general introduction into molecular modeling, showing us which computational techniques are useful to treat problems of varying scales, both in space and time. With the focus going to the modeling of molecules she told us everything there is to know about the potential energy surface, (PES)  how to investigate it using different computational methods. She discussed the differences between localized (i.e. gaussian) and plane wave basis sets and taught us how to accurately sample the PES using both molecular dynamics and normal mode analysis. As a final topic she introduced us to the world of computational spectroscopy, showing how infrared spectra can be simulated, and the limitations of this type of simulations.

With the, somewhat mysterious, presentations of Prof. Stefaan Cottenier we moved from the  realm of molecules to that of solids. In his first session, he introduced density functional theory, a method ideally suited to treat extended systems at the quantum mechanical level. And showed that as much information is present in the electron density of a system as is in its wave function. In his second session, we fully plunged in the world of solids, and we were guided, step by step, towards a full understanding of the technical details generally found in the methods section of (ab-initio) computational materials science work. Throughout this session, NaCl was used as an ever present example, and we learned that our simple high-school picture of bonding in kitchen salt is a lie-to-children. In reality, Cl doesn’t gain an extra electron by stealing it away from Na, instead it is rather the Na 3s electron which is living to far away from the Na nucleus it belongs to.

De-activating an active atom.

It could be that I’ve perhaps found out a little bit about the structure
of atoms. You must not tell anyone anything about it. . .
–Niels Bohr (1885 – 1965),
in a letter to his brother (1912)

Getting the news that a paper got accepted for publication is exciting news, but it can also be a little bit sad since it indicates the end of a project. Little over a month ago we got this great news regarding our paper for the journal of chemical information and modeling. It was the culmination of a side project Goedele Roos and I had been working on, in an on-and-off fashion, over the last two years.

When we started the project each of us had his/her own goal in mind. In my case, it was my interest in showing that my Hirshfeld-I code could handle systems which are huge from the quantum mechanical calculation point of view. Goedele, on the other hand, was interested to see how good Hirshfeld-I charges behaved with increasing size of a molecular fraction. This is of interest for multiscale modeling approaches, for which Martin Karplus, Michael Levitt, and Arieh Warshel got the Nobel prize in chemistry in 2013. In such an approach, a large system, for example a solvated biomolecule containing tens of thousands of atoms, is split into several regions. The smallest central region, containing the part of the molecule one is interested in is studied quantum mechanically, and generally contains a few dozen up to a few hundred atoms. The second shell is much larger, and is described by force-field approaches (i.e. Newtonian mechanics) and can contain ten of thousands of atoms. Even further  from the quantum mechanically treated core a third region is described by continuum models.

What about the behavior of the charges? In a quantum mechanical approach, even though we still speak of electrons as-if referring to classical objects, we cannot point to a specific point in space to indicate: “There it is”. We only have a probability distribution in space indicating where the electron may be. As such, it also becomes hard to pinpoint an atom, and in an absolute sense measure/calculate it’s charge. However, because such concepts are so much more intuitive, many chemists and physicists have developed methods, with varying success, to split the electron probability distribution into atoms again. When applying such a scheme on the probability distributions of fractions of a large biomolecule, we would like the atoms at the center not to change to much when the fraction is made larger (i.e. contain more atoms). This would indicate that from some point onward you have included all atoms that interact with the central atoms. I think, you can already see the parallel with the multiscale modeling approach mentioned above; where that point would indicate the boundary between the quantum mechanical and the Newtonian shell.

Convergence of Hirshfeld-i charges for clusters of varying size of a biomolecule.

Convergence of Hirshfeld-I charges for clusters of varying size of a biomolecule. The black curves show the charge convergence of an active S atom, while the red curves indicate a deactivated S atom.

Although, we expected to merely be studying this convergence behavior, for the particular partitioning scheme I had implemented, we dug up an unexpected treasure. Of the set of central atoms we were interested all except one showed the nice (and boring) convergence behavior. The exception (a sulfur atom) showed a clear lack of convergence, it didn’t even show any intend toward convergence behavior even for our system containing almost 1000 atoms. However, unlike the other atoms we were checking, this S atom had a special role in the biomolecule: it was an active site, i.e. the atom where chemical reactions of the biomolecule with whatever else of molecule/atom are expected to occur.

Because this S atom had a formal charge of -1, we bound a H atom to it, and investigated this set of new fractions. In this case, the S atom, with the H atom bound to it, was no longer an active site. Lo and behold, the S atom shows perfect convergence like all other atoms of the central cluster. This shows us that an active site is more than an atom sitting at the right place at the right time. It is an atom which is reaching out to the world, interacting with other atoms over a very long range, drawing them in (>10 ångström=1 nm is very far on the atomic scale, imagine it like being able to touch someone who is standing >20 m away from you). Unfortunately, this is rather bad news for multiscale modeling, since this means that if you want to describe such an active site accurately you will need an extremely large central quantum mechanical region. When the active site is deactivated, on the other hand, a radius of ~0.5 nm around the deactivated site is already sufficient.

Similar  to Bohr, I have the feeling that “It could be that I’ve perhaps found out a little bit about the structure
of atoms.”, and it makes me happy.

Convergence of Atomic Charges with the Size of the Enzymatic Environment

Authors: Danny E. P. Vanpoucke, Julianna Oláh, Frank De Proft, Veronique Van Speybroeck, and Goedele Roos
Journal: J. Chem. Inf. Model. 55(3), 564-571 (2015)
doi: 10.1021/ci5006417
IF(2015): 3.657
export: bibtex
pdf: <J.Chem.Inf.Model.> 
Graphical Abstract: The influence of the cluster size and water presence on the atomic charge of active and inactive sites in Biomolecules.
Graphical Abstract: Graphical Abstract: The influence of the cluster size and water presence on the atomic charge of active and inactive sites in Bio-molecules.

Abstract

Atomic charges are a key concept to give more insight into the electronic structure and chemical reactivity. The Hirshfeld-I partitioning scheme applied to the model protein human 2-cysteine peroxiredoxin thioredoxin peroxidase B is used to investigate how large a protein fragment needs to be in order to achieve convergence of the atomic charge of both, neutral and negatively charged residues. Convergence in atomic charges is rapidly reached for neutral residues, but not for negatively charged ones. This study pinpoints difficulties on the road towards accurate modeling of negatively charged residues of large bio-molecular systems in a multiscale approach.

39th ICACC: day 3-5

The last three days of the conference, the virtual materials design session took place. This session was specifically focused on computational materials design. Because of this focus, the attendance was rather low, mainly computational materials scientists. Apparently, this type of specialized focus on computational work is the best way not to reach the general experimental public in the same field. As a computational scientist the only way to circumvent this, is by applying for a presentation in a relevant experimental session. This requires you to make a less technical presentation, but this is not a bad thing, since it forces you to think about your results and understand them in more simple terms.

An example of such a presentation was given Dr. Ong who discussed his high-throughput ab initio setup for designing solid state electrolytes. He showed that a material (Li10GeP2S12) which was thought to be a 1D conductor, actually is a 3D conductor, however, the Li-conductivity in the directions perpendicular to the 1D direction is ~100 times smaller, explaining why they were not noticed before. He also presented newly predicted materials for Li-transport, which lead to a standard experimental remark that such computational predictions mean very little, since they do not take into account temperature, and as such these structures may not be  stable. Although such remarks are “in theory” true, and are a nice example of a lack of understanding outside the computational community, in this case, the material the experimental researcher was referring to was recently synthesized and found to be stable.

On Thursday morning, I had the opportunity to present my contributed presentation. In contrast to my invited presentation, this presentation was solely focused on doped cerium dioxide. Using a three-step approach, we investigated all different contributions of the dopants to their modification of the mechanical properties of CeO2. In the first step, we look at group IV dopants, since Ce has an oxidation state of +IV in CeO2. Here we show that the character of the valence electrons (p or d) plays an important role with regard to stability and mechanical properties. In our second step dopants with an oxidation state different from +IV are considered, without the presence of oxygen vacancies. In this case, the same trends and behavior is observed as in the first step. In the third and last step also oxygen vacancies are included. We show that oxygen vacancies have a stabilizing influence of the doped system. Furthermore, the oxygen vacancies make CeO2 mechanically softer, i.e. reducing the bulk modulus.

In the afternoon, Prof. Frederico Rosei, of NanoFemto lab in Canada, gave an entertaining lecture on “Mentorship for young scientists: Developing scientific survival skills”. It presented an interesting forum to find out that, as scientists, we all seem to struggle with the same things. We want to do what we like (research), and invest a lot in this, unfortunately external forces (the struggle for funding/job security) complicate life. Frederico centered his lecture around three points of importance/goals a young scientist should always try to be aware of:

  • Know yourself
  • Plan ahead
  • Find a mentor

Although these are lofty goals, they tend to be quite none-trivial in the current-day scientific environment. Finding a mentor, i.e. a senior scientist with time on his/her hands not involved in your projects, is a bit like looking for a unicorn. Unlike the unicorn, they do exist, but there are very very few of them (How many professors do you know with spare time?). Planning ahead, and following your own plans are nice in theory, unfortunately a young scientists’ life (i.e. everyone below tenured professorship) tends to be ruled by funding in a kind of life or dead setup. I am not saying this is not the case for tenured professors, however, it is not their own life and death. For all other scientist : No funding=no job=end scientific career. As such, the pressure to publish (yes, funding agencies only count your papers, not how good/bad they are even though that is the official statement) is high, and will have a detrimental influence on the quality of science and of what is being published (if it isn’t already the case). I truly wished, the world could be as Prof. Rosei envisions it. Back to more happy subjects.

Friday was the last day of the conference. In the morning I again attended the virtual materials design session, but as with all other sessions several presentations were canceled, apparently snow is wreaking some havoc in New York airport, preventing several presenters not to be able to make it. Luckily Eva Zarkadoula made it to the conference to present her very nice modeling work: “Molecular Dynamic Simulations of Synergistic Effects in Ion Track Formation”. Using classical molecular dynamics simulations, she simulated how incoming high energy radiation traces a path trough a material allowing one to use this material as a detector. In contrast to what I would have imagined, perfectly crystalline material shows very little damage after the radiation has passed through. Even though initially a clear trace is visible, the system appear to relax back to a more or less perfect crystalline solid upon relaxation, making it a rather poor sensor material. However, if defects are present in the material, the track made by the radiation remains clearly visible. An extremely nice bonus to this work is the fact that direct comparison to experiments is possible.

The conference ended at noon, leaving some time to have a walk on the beach, find some souvenirs, and have a last dinner with colleagues from the conference.

 

39th ICACC: Day 2

Today I had to get to work myself, unlike the first day of the ICACC where I was only a spectator. At 8h30 I kicked of the GYIF as the first speaker in a session called “Theoretical Modeling and Applications” with my presentation: Computational Materials Science: Where Theory meets Experiment. Intended for a general audience of (expectedly mainly) experimental materials scientist this presentation was aimed to show the audience that computational materials science has in the last decades developed to the level where results relevant for real-life applications in materials science can be obtained. For this I used three examples of my own work: (1) Pt nanowires on Ge(001), where I used simulated STM to build an accurate atomistic model of these wires, (2) doped cerium oxides, where the influence of doping on lattice parameter and thermal expansion coefficients was studied for the purpose of matching them to those of other materials, and (3) Metal-Organic Frameworks, where I showed that the spin-configuration of the MIL-47(V) MOF is linked to the transition pressure inducing breathing.

Later on during the morning, I had the pleasure of chairing the session “Additive Manufacturing” together with Valerie Wiesner from NASA. During this session Mahref Vali and Lisa Rueschoff presented their most recent work on 3D printing of ceramic materials, a technique which will allow the printing of ceramic components in the future. The third speaker of this session, Rumi Kitazawa, delighted us with an inspiring talk on the “Engineering applications of Menger sponges”. A Menger sponge is a fractal related to the Cantor set, as such a fully developed Menger sponge can not be build with any real material, however, using 3D printing it is possible to build a structure with Menger sponge like features (i.e. with holes down to a certain size). By comparing experimental stress tests on such 3D-printed systems with calculations of the strain energy in such a structure, Rumi was able to show these Menger materials have a peculiar, albeit very organised, strain pattern along the main diagonal of the material. The combination of large and small pores present in these Menger sponge materials may make this behavior relevant for MOF (and other porous) materials, where the large pores reflect inter-grain pores, while the small pores are the pores of the MOF. So this is definitely a topic to remember.

Next to presentations, conferences also contain social events. Today, there were two social events; at noon there was a luncheon by the GYIF where the young investigators could mix with people from industry and senior group leaders. In the evening there was the first poster session and booth-stand where companies try to sell their services and lab equipment. Somehow, as a theoretician, I am always a bit at a loss at such events. To draw in more people, there was also a shot-glass contest. No, it was not the goal to drink as much as possible, but to build a protective structure around a shot-glass using only 15 drinking straws (no tape, wire, paper, staples,… allowed). To find the most protective structure, the shot-glass and straw constructions were dropped from various heights. Twenty four teams started at a drop height of 3 feet (~1 m), where already the first shot-glass didn’t survive the drop. Every round, the drop height was increased by 3 feet. For a drop of 20 feet (~6 m) there were only two teams remaining, including our team. Altough our glass survived its first bounce, the second bounce unfortunately broke our glass (darn). Then it was the turn of our remaining competitor, who’s shot-glass exploded into shards on first impact. Officially the result was a draw, although it is clear our construct had clearly the upper hand 😎 .

Our mixed Theoretical-Experimental "international-multi-university" research  team. Left to right: me (UGhent, Belgium), Bert Conings (UHasselt, Belgium), and Chenxin Jin  (Dalhousie University, Canada)  On the right hand side, you can see our construct of straws  around the shot glass.

Our mixed Theoretical-Experimental “international-multi-university” research team. Left to right: me (UGhent, Belgium), Bert Conings (UHasselt, Belgium), and Chenxin Jin (Dalhousie University, Canada)
On the right hand side, you can see our construct of straws around the shot glass.

39th ICACC: Day 1

Today the 39th International Conference  and exposition on Advanced Ceramics and Composites (ICACC) started in Daytona Beach, Florida. Here, scientist from all over the world will be discussing their latest work and findings in the field of ceramic materials during the coming week. Although my project on ceramic materials has been completed for over 2 years now, culminating in a second PhD, I was invited to present my work here. As with many of this kind of conferences, I am afraid that, as a computational materials scientist, I belong to a minority, strongly outnumbered by the experimental (materials) scientists present. This is an aspect that I will need to consider preparing my presentations.

The first morning session consisted of four plenary lectures (general overview presentations in which celebrated group leaders present the overall picture of the work done in their group and their hopes/views on the future). We started of an interesting lecture on “Thermal Barrier Coatings for Gas Turbines” by Prof. David Clarke, where we learned that, since a major part of the world wide energy production is gas based, improving gas-turbine efficiency by only 1% would produce more energy than all renewable production currently in play. This efficiency improvement can be obtained by operating the gas-turbines at higher temperatures. Unfortunately, the metal fans of such a turbine start to degrade if temperatures are too high. By coating them with materials that have a low thermal conductivity, it is possible to operate at the required high temperatures, while the metal fans experience an acceptable operating temperature which is a few hundred degrees lower. Next Prof. Sanjay Mathur from the university of Cologne presented his  groups work on the development of precursor libraries (these ideas are similar as those behind computational high-throughput projects). In precursor chemistry, where changing functional groups or doping leads to changes in the surface morphology, such libraries would then present an interesting tool for designing new materials for energy and health applications. As this apparently is a very hot topic, the fire-alarm of the conference center went into overdrive and everyone needed to be evacuated. When the conference resumed half an hour later Prof. Mathur reassured us he did not intend to fire things up like this. His presentation was followed by that of Prof. Cato Laurencin, who showed us how ceramic materials could be used in a new field he wishes to launch: “Regenerative engineering”. Here, combinations of micro- and nanostructured ceramics are used as matrices to grow and differentiate stem cells intended to heal fractured bones and cartilage. The final presentation by Prof. Kazushige Ohno discussed next generation filters for diesel particulates, which should provide us with a cleaner future.

In the afternoon the parallel symposia started, where I followed the 4th Global Young Investigator Forum (GYIF). Here, Prof. Ricardo Castro presented an interesting method for experimentally obtaining the surface energy of nanoparticles. His quest originated from the simple observation that existing phase diagrams for bulk materials no longer hold when one is working with nano-particles. In such systems, the energy contributions due to the surface of the particle become comparable to those of the inner bulk. Interestingly, one of the example systems Prof. Castro looked into was Mn3+ doped CeO2. In his work he found that the Mn was mainly located at the surface of the CeO2 particles, something I also expected from my own work on aliovalent doped CeO2, based on the defect formation energies of Cu and Co doping. Further presentations discussed Organometal trihalide perovskite solar cells. Although these solar cells still are rather unstable, they do show promise with regard to their efficiency.(The origin of this efficiency is unfortunately not really understood. Maybe other perovskite MOFs are more stable?)

Comment on ‘Europium doping induced symmetry deviation and its impact on the second harmonic generation of doped ZnO nanowires’

Authors: Danny E. P. Vanpoucke
Journal: Nanotechnology 25(45), 458001 (2014)
doi: 10.1088/0957-4484/25/45/458001
IF(2014): 3.821
export: bibtex
pdf: <Nanotechnology>

Abstract

In Dhara et al. 2014 Nanotechnology 25 225202, the authors reported on the synthesis of Eu-doped ZnO nanowires (NWs) and investigated the influence of Eu doping on the second harmonic generation (SHG). Maximum SHG was found to correlate strongly with the structural deformation attributed to Eu3+ doping. In this comment, we show the deformation of interest is due to the presence of Eu2+ dopants, based on both the experimental data presented by Dhara et al. and ab-initio density functional theory calculations.

Quasi-1D physics in metal-organic frameworks: MIL-47(V) from first principles

Authors: Danny E. P. Vanpoucke, Jan W. Jaeken, Stijn De Baerdemacker, Kurt Lejaeghere
and Veronique Van Speybroeck
Journal: Beilstein J. Nanotechnol. 5, 1738-1748 (2014)
doi: 10.3762/bjnano.5.184
IF(2014): 2.670
export: bibtex
pdf: <Beilstein> (open access)
Graphical Abstract: (left) Spin density of anti-ferromagnetic MIL-47(V) with ferromagnetic chains. (right) Electronic band structure and density of states.
Graphical Abstract: The MIL-47(V) MOF has one unpaired electron per V site. As a result, different spin configurations are possible, several of which lead to an anti-ferromagnetic state. The spin density of an antiferromagnetic state, containing only ferromagnetic chains is shown on the left. On the right, the electronic band structure of the same system is presented.

Abstract

The geometric and electronic structure of the MIL-47(V) metal-organic framework (MOF) is investigated by using ab initio density functional theory (DFT) calculations. Special focus is placed on the relation between the spin configuration and the properties of the MOF. The ground state is found to be antiferromagnetic, with an equilibrium volume of 1554.70 Å3. The transition pressure of the pressure-induced large-pore-to-narrow-pore phase transition is calculated to be 82 MPa and 124 MPa for systems with ferromagnetic and antiferromagnetic chains, respectively. For a mixed system, the transition pressure is found to be a weighted average of the ferromagnetic and antiferromagnetic transition pressures. Mapping DFT energies onto a simple-spin Hamiltonian shows both the intra- and inter-chain coupling to be antiferromagnetic, with the latter coupling constant being two orders of magnitude smaller than the former, suggesting the MIL-47(V) to present quasi-1D behavior. The electronic structure of the different spin configurations is investigated and it shows that the band gap position varies strongly with the spin configuration. The valence and conduction bands show a clear V d-character. In addition, these bands are flat in directions orthogonal to VO6 chains, while showing dispersion along the the direction of the VO6 chains, similar as for other quasi-1D materials.

Aliovalent Doping of CeO2: DFT study of oxidation state and vacancy effects

Authors: Danny E. P. Vanpoucke, Patrick Bultinck, Stefaan Cottenier, Veronique Van Speybroeck, and Isabel Van Driessche,
Journal: J. Mater. Chem. A 2(33), 13723-13737 (2014)
doi: 10.1039/C4TA02449D
IF(2014): 7.443
export: bibtex
pdf: <JMaterChemA> <arXiv>

Abstract

The modification of CeO2 properties by means of aliovalent doping is investigated within the ab initio density functional theory framework. Lattice parameters, dopant atomic radii, bulk moduli and thermal expansion coefficients of fluorite type Ce1-xMxO2-y (with M = Mg, V, Co, Cu, Zn, Nb, Ba, La, Sm, Gd, Yb, and Bi) are presented for 0.00 ≤ x ≤ 0.25. The relative stability of the doped systems is discussed, and the influence of oxygen vacancies is investigated. It is shown that oxygen vacancies tend to increase the lattice parameter, and strongly decrease the bulk modulus. Defect formation energies are correlated with calculated crystal radii and covalent radii of the dopants, and are shown to present no simple trend. The previously observed inverse relationship between the thermal expansion coefficient and the bulk modulus in group IV doped CeO2 [J. Am. Ceram. Soc. 97(1), 258 (2014)] is shown to persist independent of the inclusion of charge compensating vacancies.

Tetravalent Doping of CeO2: The impact of valence electron character on group IV dopant influence

Authors: Danny E. P. Vanpoucke, Stefaan Cottenier, Veronique Van Speybroeck, Isabel Van Driessche, and Patrick Bultinck
Journal: J. Am. Ceram. Soc. 97(1), 258-266 (2014)
doi: 10.1111/jace.12650
IF(2014): 2.610
export: bibtex
pdf: <J.Am.Ceram.Soc.> <arXiv>

Abstract

Fluorite CeO2 doped with group IV elements is studied within the density functional theory (DFT) and DFT + U framework. Concentration-dependent formation energies are calculated for Ce1−xZxO2 (Z = C, Si, Ge, Sn, Pb, Ti, Zr, Hf) with 0 ≤ x ≤ 0.25 and a roughly decreasing trend with ionic radius is observed. The influence of the valence and near valence electronic configuration is discussed, indicating the importance of filled d and f shells near the Fermi level for all properties investigated. A clearly different behavior of group IVa and IVb dopants is observed: the former are more suitable for surface modifications and the latter are more suitable for bulk modifications. For the entire set of group IV dopants, there exists an inverse relation between the change, due to doping, of the bulk modulus, and the thermal expansion coefficients. Hirshfeld-I atomic charges show that charge-transfer effects due to doping are limited to the nearest-neighbor oxygen atoms.