The role of superconductivity in development of the Standard Model

In 1986, Steven Weinberg published an article, Superconductivity for Particular Theorists, in which he stated

"No one did more than Nambu to bring the idea of spontaneously broken symmetries to the attention of elementary particle physicists. And, as he acknowledged in his ground-breaking 1960 article  "Axial Current Conservation in Weak Interactions'', Nambu was guided in this work by an analogy with the theory of superconductivity,..."

In the 1960 PRL, referenced by Weinberg, Nambu states that in the BCS theory, as refined by Bogoliubov, [and Anderson]

"gauge invariance, the energy gap, and the collective excitations are logically related to each other as was shown by the author. [Y. Nambu, Phys. Rev. 117, 648 (1960)] In the present case we have only to replace them by (chiral) (gamma_5) invariance, baryon mass, and the mesons." 

This connection is worked out explicitly in two papers in 1961. The first is

Dynamical Model of Elementary Particles Based on an Analogy with Superconductivity. I

Y. Nambu and G. Jona-Lasinio

They acknowledge, 

"that the model treated here is not realistic enough to be compared with the actual nucleon problem. Our purpose was to show that a new possibility exists for field theory to be richer and more complex than has been hitherto envisaged,"

Hence, I consider this to be a toy model for an emergent phenomena.

The model consists of a massless fermion field with a quartic interaction that has chiral invariance, i.e., unchanged by global gauge transformations associated with the gamma_5 matrix. (The Lagrangian is given above.) At the mean-field level, this symmetry is broken. Excitations include massless bosons (associated with the symmetry breaking and similar to those found earlier by Goldstone) and bound fermion pairs. It was conjectured that these could be analogues of mesons and baryons, respectively. The model was proposed before quarks and QCD. Now, the fermion degrees of freedom would be identified with quarks, and the model illustrates the dynamical generation of quark masses. When generalised to include SU(2) or SU(3) symmetry the model is considered to be an effective field theory for QCD, such as chiral effective theory.

Multi-step spin-state transitions in organometallics and frustrated antiferromagnetic Ising models

In previous posts, I discussed how "spin-crossover" material is a misnomer because many of these materials do not undergo crossovers but phase transitions due to collective effects. Furthermore, they exhibit rich behaviours, including hysteresis, incomplete transitions, and multiple-step transitions. Ising models can capture some of these effects.

Here, I discuss how an antiferromagnetic Ising model with frustrated interactions can give multi-step transitions. This has been studied previously by Paez-Espejo, Sy and Boukheddaden, and my UQ colleagues Jace Cruddas and Ben Powell. In their case, they start with a lattice "balls and spring" model and derive Ising models with an infinite-range ferromagnetic interaction and short-range antiferromagnetic interactions. They show that when the range of these interactions (and thus the frustration) is increased, more and more steps are observed.

Here, I do something simpler to illustrate some key physics and some subtleties and cautions.

fcc lattice

Consider the antiferromagnetic Ising model on the face-centred-cubic lattice in a magnetic field. 

[Historical trivia: the model was studied by William Shockley back in 1938, in the context of understanding alloys of gold and copper.]

The picture below shows a tetrahedron of four nearest neighbours in the fcc lattice.

Even with just nearest-neighbour interactions, the lattice is frustrated. On a tetrahedron, you cannot satisfy all six AFM interactions. Four bonds are satisfied, and two are unsatisfied.

The phase diagram of the model was studied using Monte Carlo by Kammerer et al. in 1996. It is shown above as a function of temperature and field. All the transition lines are (weakly) first-order.

The AB phase has AFM order within the [100] planes. It has an equal number of up and down spins.

The A3B phase has alternating FM and AFM order between neighbouring planes. Thus, 3/4 of the spins have the same direction as the magnetic field.

The stability of these ordered states is subtle. At zero temperature, both the AB and A3B states are massively degenerate. For a system of 4 x L^3 spins, there are 3 x 2^2L AB states, and 6 x 2^L   A3B states. At finite temperature, the system exhibits “order by disorder”.

On the phase diagram, I have shown three straight lines (blue, red, and dashed-black) representing a temperature sweep for three different spin-crossover systems. The "field" is given by h=1/2(Delta H - T Delta S). In the lower panel, I have shown the temperature dependence of the High Spin (HS) population for the three different systems. For clarity, I have not shown the effects of the hysteresis associated with the first-order transitions.

If Delta H is smaller than the values shown in the figure, then at low temperatures, the spin-crossover system will never reach the complete low-spin state.

Main points.

Multiple steps are possible even in a simple model. This is because frustration stabilises new phases in a magnetic field. Similar phenomena occur in other frustrated models, such as the triangular lattice, the J1-J2 model on a chain or a square lattice.

The number of steps may change depending on Delta S. This is because a temperature sweep traverses the field-temperature phase diagram asymmetrically.

Caution.

Fluctuations matter.

The mean-field theory phase diagram was studied by Beath and Ryan. Their phase diagram is below. Clearly, there are significant qualitative differences, particularly in the stability of the A3B phase.

The transition temperature at zero field is 3.5 J, compared to the value of 1.4J from Monte Carlo.

Monte Carlo simulations may be fraught.

Because of the many competing ordered states associated with frustration, Kammerer et al. note that “in a Monte Carlo simulation one needs unusually large systems in order observe the correct asymptotic behaviour, and that the effect gets worse with decreasing temperature because of the proximity of the phase transition to the less ordered phase at T=0”. 

Open questions.

The example above hints at what the essential physics may be how frustrated Ising models may capture it. However, to definitively establish the connection with real materials, several issues need to be resolved.

1. Show definitively how elastic interactions can produce the necessary Ising interactions. In particular, derive a formula for the interactions in terms of elastic properties of the high-spin and low-spin states. How do their structural differences, and the associated bond stretches or compressions, affect the elastic energy? What is the magnitude, range, and direction of the interactions?

[n.b. Different authors have different expressions for the Ising interactions for a range of toy models, using a range of approximations. It also needs to be done for a general atomic "force field".]

2. For specific materials, calculate the Ising interactions from a DFT-based method. Then show that the relevant Ising model does produce the steps and hysteresis observed experimentally.

"Ferromagnetic" Ising models for spin-state transitions in organometallics

In recent posts, I discussed how "spin crossover" is a misnomer for the plethora of organometallic compounds that undergo spin-state phase transitions (abrupt, first-order, hysteretic, multi-step,...)

In theory development, it is best to start with the simplest possible model and then gradually add new features to the model until (hopefully) arriving at a minimal model that can describe (almost) everything. Hence, I described how the two-state model can describe spin crossover. An Ising "spin" has values of +1 or -1, corresponding to high spin (HS) and low spin (LS) states. The "magnetic" field is half of the difference in Gibbs free energy between the two states. 

The model predicts equal numbers of HS and LS at a temperature

The two-state model is modified by adding Ising-type interactions between the “spins” (molecules). The Hamiltonian is then of the form

 The temperature dependence in the field arises because this is an effective Hamiltonian.

The Ising-type interactions are due to elastic effects. The spin-state transition in the iron atom leads to changes in the Fe-N bond lengths (an increase of about 10 per cent in going from LS to HS), changing the size of the metal-ligand (ML6 ) complex. This affects the interactions (ionic, pi-pi, H-bond, van der Waals) between the complexes. The volume of the ML6 complex changes by about 30 per cent, but typically the volume of the crystal unit cell changes by only a few per cent. The associated relaxation energies are related to the J’s. Calculating them is non-trivial and will be discussed elsewhere. There are many competing and contradictory models for the elastic origin of the J’s.

In this post, I only consider nearest-neighbour ferromagnetic interactions. Later, I will consider antiferromagnetic interactions and further-neighbour interactions that lead to frustration. 

Slichter-Drickamer model

This model was introduced in 1972 is beloved by experimentalists, especially chemists, because it provides a simple analytic formula that can be fit to experimental data.

The system is assumed to be a thermodynamic mixture of HS and LS. x=n_HS(T) is the fraction of HS. The Gibbs free energy is given by

This is minimised as a function of x to give the temperature dependence of the HS population.

The model is a natural extension of the two-state model, by adding a single parameter, Gamma, which is sometimes referred to as the cooperativity parameter.

The model is equivalent to the mean-field treatment of a ferromagnetic Ising model, with Gamma=2zJ, where z is the number of nearest neighbours. Some chemists do not seem to be aware of this connection to Ising. The model is also identical to the theory of binary mixtures, such as discussed in Thermal Physics by Schroeder, Section 5.4.

Successes of the model.

good quantitative agreement with experiments on many materials.

a first-order transition with hysteresis for T_1/2 < Tc =z J.

a steep and continuous (abrupt) transition for T_1/2 slightly larger than Tc.

Values of Gamma are in the range 1-10 kJ/mol. Corresponding vaules of J are in the range 10-200 K, depending on what value of z is assumed.

Weaknesses of the model.

It cannot explain multi-step transitions.

Mean-field theory is quantitatively, and sometimes qualitatively, wrong, especially in one and two dimensions.

The description of hysteresis is an artefact of the mean-field theory, as discussed below.

Figure. Phase diagram of a ferromagnetic Ising model in a magnetic field. (Fig. 8.7.1, Chaikin and Lubensky). Vertical axis is the magnetic field, and the horizontal axis is temperature. Tc denotes the critical temperature, and the double-line denotes a first-order phase transition between paramagnetic phases where the magnetisation is parallel to the direction of the applied field.

Curves show the free energy as a function of the order parameter (magnetisation) in mean-field theory. The dashed lines are the lines of metastability deduced from these free-energy curves. Inside these lines, the free energy has two minima: the equilibrium one and a metastable one. The lines are sometimes referred to as spinodal curves.

The consequences of the metastability for a field sweep at constant temperature are shown in the Figure below, taken from Banerjee and Bar.

How does this relate to thermally induced spin-state transitions?

Consider the phase diagram shown above of a ferromagnetic Ising model in a magnetic field. The red and blue lines correspond to temperature scans for two SCO materials that have different values of the parameters Delta H and DeltaS.

The occurrence of qualitatively different behaviour is determined by where the lines intercept the temperature and field axes, i.e. the values of T_1/2 /J and Delta H/J. If the former is larger than Tc/J, as it is for the blue line, then no phase transition is observed. 

The parameter Delta H/J determines whether at low temperatures, the complete HS state is formed.

The figure below is a sketch of the temperature dependence of the population of HS for the red and blue cases.

Note that because of the non-zero slope of the red line, the temperature  T_1/2 is not the average of the temperatures at which the transition occurs on the up and down temperature sweeps.

Deconstructing hysteresis.

The physical picture above of metastability is an artefact (oversimplification) of mean-field theory. It predicts that an infinite system would take an infinite time to reach the equilibrium state from the metastable state.

(Aside: In the context of the corresponding discrete-choice models in economics, this has important and amusing consequences, as discussed by Bouchaud.)

In reality, the transition to the equilibrium state can occur via nucleation of finite domains or in some regimes via a perturbation with a non-zero wavevector. This is discussed in detail by Chaikin and Lubensky, chapter 4.

The consequence of this “metastability” for a first-order transition in an SCO system is that the width of the hysteresis region (in temperature) may depend on the rate at which the temperature is swept and whether the system is allowed to relax before the magnetisation (fraction of HS) is measured at any temperature. Emprically, this is observed and has been highlighted by Brooker, albeit without reference to the theoretical subtleties I am highlighting here. She points out that up to 2014, chemists seemed to have been oblivious to these issues and reported results without testing whether their observations depended on the sweep rate or whether they waited for relaxation.

(Aside. The dynamics are different for conserved and non-conserved order parameters. In a binary liquid mixture, the order parameter is conserved, i.e., the number of A and B atoms is fixed. In an SCO material, the number of HS and LS is not conserved.)

In the next post, I will discuss how an antiferromagnetic Ising model can give a two-step transition and models with frustrated interactions can give multi-step transitions.

The two-state model for spin crossover in organometallics

Previously, I discussed how spin-crossover is a misnomer for organometallic compounds and proposed that an effective Hamiltonian to describe the rich states and phase transitions is an Ising model in "magnetic field".

I introduce the two-state model that defines the model without the Ising interactions. To save me time on formatting in HTML, here is a pdf file that describes the model and what comparisons with experimental data (such as that below) tells us.

Future posts will consider how elastic interactions produce the Ising interaction and how frustrated interactions can produce multi-step transitions.

My review article on emergence

I just posted on the arXiv a long review article on emergence

Emergence: from physics to biology, sociology, and computer science

The abstract is below.

I welcome feedback. 

------

Many systems of interest to scientists involve a large number of interacting parts and the whole system can have properties that the individual parts do not. The system is qualitatively different to its parts. More is different. I take this novelty as the defining characteristic of an emergent property. Many other characteristics have been associated with emergence are reviewed, including universality, order, complexity, unpredictability, irreducibility, diversity, self-organisation, discontinuities, and singularities. However, it has not been established whether these characteristics are necessary or sufficient for novelty. A wide range of examples are given to show how emergent phenomena are ubiquitous across most sub-fields of physics and many areas of biology and social sciences. Emergence is central to many of the biggest scientific and societal challenges today. Emergence can be understood in terms of scales (energy, time, length, complexity) and the associated stratification of reality. At each stratum (level) there is a distinct ontology (properties, phenomena, processes, entities, and effective interactions) and epistemology (theories, concepts, models, and methods). This stratification of reality leads to semi-autonomous scientific disciplines and sub-disciplines. A common challenge is understanding the relationship between emergent properties observed at the macroscopic scale (the whole system) and what is known about the microscopic scale: the components and their interactions. A key and profound insight is to identify a relevant emergent mesoscopic scale (i.e., a scale intermediate between the macro- and micro- scales) at which new entities emerge and interact with one another weakly. In different words, modular structures may emerge at the mesoscale. Key theoretical methods are the development and study of effective theories and toy models. Effective theories describe phenomena at a particular scale and sometimes can be derived from more microscopic descriptions. Toy models involve minimal degrees of freedom, interactions, and parameters. Toy models are amenable to analytical and computational analysis and may reveal the minimal requirements for an emergent property to occur. The Ising model is an emblematic toy model that elucidates not just critical phenomena but also key characteristics of emergence. Many examples are given from condensed matter physics to illustrate the characteristics of emergence. A wide range of areas of physics are discussed, including chaotic dynamical systems, fluid dynamics, nuclear physics, and quantum gravity. The ubiquity of emergence in other fields is illustrated by neural networks, protein folding, and social segregation. An emergent perspective matters for scientific strategy, as it shapes questions, choice of research methodologies, priorities, and allocation of resources. Finally, the elusive goal of the design and control of emergent properties is considered.

Spin crossover is a misnomer

There are hundreds of organometallic compounds that are classified as spin-crossover compounds. As the temperature is varied the average spin per molecule can undergo a transition between low-spin and high-spin states.

The figure below shows several classes of transitions that have been observed. The vertical axis represents the fraction of molecules in the high-spin state, and the horizontal axis represents temperature.

Taken from Gutlich, Garcia, and Spiering and reproduced by Nicolazzi and Bousseksou

a) A smooth crossover. At the temperature T_{1/2} there are equal numbers of high and low spins.

b) There is sharp transition with the curve having a very large slope at T_{1/2}.

c) There is a discontinuous change in the spin fraction at the transition temperature, the value of which depends on whether the temperature is increasing or decreasing, i.e., there is hysteresis. The discontinuity and hysteresis are characteristic of a first-order phase transition.

d) There is a step in the curve when the high-spin fraction is close to 0.5. This is known as a two-step transition.

e) Although a crossover occurs, the system never contains only low- or high-spins.

But, there is more. Over the past decade, multiple-step transitions have been observed. An example of a four-step transition is below.

The figure is taken from a paper by Clements, Price, Neville, and Kepert.

Hysteresis is present and is larger at lower temperatures.

In a few cases of multiple-step transitions on the down-temperature sweep, the first step is missing compared to the up-temperature step.

Given the diverse behaviour described above, including sharp transitions and first-order phase transitions, spin "crossover" is a misnomer.

More importantly, given the chemical and structural complexity materials involved, is there a simple model effective Hamiltonian that can capture all this diverse behaviour?

Yes. An Ising model in a field. A preliminary discussion is here. I hope to discuss this in future posts. But first I need to introduce the simple two-state model and show what it can and cannot explain.

Science job openings in sunny Brisbane, Australia

Bribie Island, just north of Brisbane.

The University of Queensland has just advertised several jobs that may be of interest to readers of this blog, particularly those seeking to flee the USA.

There is a junior faculty position for a theorist working at the interface of condensed matter, quantum chemistry, and quantum computing.

There is also a postdoc to work on the theory of strongly correlated electron systems with my colleagues Ben Powell and Carla Verdi.

There is a postdoc in experimental condensed matter, to work on scanning probe methods, such as STM, with my colleague Peter Jacobson.

Glasshouse Mountains. Just north of Brisbane.