Search

Switchability of materials properties by applying controlled stimuli such as voltage pulses is an emerging field of study with applicability in adaptive and programmable devices like neuromorphic transistors or non-emissive smart displays. One of the most exciting approaches to modulate materials performance is mobile ion/vacancy insertion for inducing changes in relevant electrical, optical, or magnetic properties, among others. Unveiling the interplay between changes in the concentration of mobile defects (like oxygen vacancies) and functional properties in relevant materials represents a step forward for underpinning the emerging oxide iontronics discipline. In this work, electrochemical oxide-ion pumping cells were fabricated for an analog control of the oxygen stoichiometry in thin films of mixed ionic-electronic conductor La0.5Sr0.5FeO3-{\delta}. We demonstrate over more than 4 orders of magnitude electronic conductivity control within the same crystallographic phase through the precise and continuous voltage control of the oxygen stoichiometry. We show that behind the modification of the transport properties of the material lays a paramagnetic-to-antiferromagnetic transition. We exploit such magnetoelectric coupling to show control over the exchange interaction between La0.5Sr0.5FeO3-{\delta} and a ferromagnetic Co layer deposited on top.

Neuromorphic hardware facilitates rapid and energy-efficient training and operation of neural network models for artificial intelligence. However, existing analog in-memory computing devices, like memristors, continue to face significant challenges that impede their commercialization. These challenges include high variability due to their stochastic nature. Microfabricated electrochemical synapses offer a promising approach by functioning as an analog programmable resistor based on deterministic ion-insertion mechanisms. Here, we developed an all-solid-state oxide-ion synaptic transistor employing $\text{Bi}_2\text{V}_{0.9}\text{Cu}_{0.1}\text{O}_{5.35}$ as a superior oxide-ion conductor electrolyte and $\text{La}_\text{0.5}\text{Sr}_\text{0.5}\text{F}\text{O}_\text{3-$\delta$}$ as a variable resistance channel able to efficiently operate at temperatures compatible with conventional electronics. Our transistor exhibits essential synaptic behaviors such as long- and short-term potentiation, paired-pulse facilitation, and post-tetanic potentiation, mimicking fundamental properties of biological neural networks. Key criteria for efficient neuromorphic computing are satisfied, including excellent linear and symmetric synaptic plasticity, low energy consumption per programming pulse, and high endurance with minimal cycle-to-cycle variation. Integrated into an artificial neural network (ANN) simulation for handwritten digit recognition, the presented synaptic transistor achieved a 96% accuracy on the MNIST dataset, illustrating the effective implementation of our device in ANNs. These findings demonstrate the potential of oxide-ion based synaptic transistors for effective implementation in analog neuromorphic computing based on iontronics.

Functional properties of mixed ionic electronic conductors (MIECs) can be radically modified by (de)insertion of mobile charged defects. A complete control of this dynamic behaviour has multiple applications in a myriad of fields including advanced computing, data processing, sensing or energy conversion. However, the effect of different MIECs state-of-charge is not fully understood yet and there is a lack of strategies for fully controlling the defect content in a material. In this work we present a model-less technique to characterize ionic defect concentration and ionic insertion kinetics in MIEC materials: Optoionic Impedance Spectroscopy (OIS). The proof of concept and advantages of OIS are demonstrated by studying the oxygen (de)insertion in thin films of hole-doped perovskite oxides. Ion migration into/out of the studied materials is achieved by the application of an electrochemical potential, achieving stable and reversible modification of its optical properties. By tracking the dynamic variation of optical properties depending on the gating conditions, OIS enables to extract electrochemical parameters involved in the electrochromic process. The results demonstrate the capability of the technique to effectively characterize the kinetics of single- and even multi- layer systems. The technique can be employed for studying underlying mechanisms of the response characteristics of MIEC-based devices.

The exploitation of advanced materials for novel energy, health and computing applications requires fundamental understanding of enabling physicochemical mechanisms, including ionic and electronic conductivity, defect formation processes and reaction kinetics. Therefore, access to underlying constants of functional materials via advanced but straightforward experimental techniques is key. We present a novel, cheap and widely applicable approach to analyze oxygen-tracer-diffusion in thin films with unprecedented time resolution based on the novel in situ isotope-exchange Raman spectroscopy (IERS) methodology. Raman spectroscopy is sensitive to changes in the local isotopic composition, manifested by a frequency shift of the oxygen Raman modes. Employing a Raman transparent capping layer allows to establish an in-plane tracer gradient and follow the isotope exchange and diffusion processes via consecutive spatial and time resolved in situ Raman line scans. Mass-transport coefficients are calculated from these isotopic gradients, similar to conventional techniques, but with an additional time-component, not accessible by other techniques. Here, we study gadolinium doped ceria (CGO) thin films, capped with Si3N4 or Al2O3. We report diffusion coefficients within the temperature range of interest for intermediate temperature emerging applications (300-500{\deg}C) and confirm the validity of the measurement procedure and extracted parameters by comparison with FEM simulations and literature results.

Holography relates gravitational theories in five dimensions to four-dimensional quantum field theories in flat space. Under this map, the equation of state of the field theory is encoded in the black hole solutions of the gravitational theory. Solving the five-dimensional Einstein's equations to determine the equation of state is an algorithmic, direct problem. Determining the gravitational theory that gives rise to a prescribed equation of state is a much more challenging, inverse problem. We present a novel approach to solve this problem based on physics-informed neural networks. The resulting algorithm is not only data-driven but also informed by the physics of the Einstein's equations. We successfully apply it to theories with crossovers, first- and second-order phase transitions.

LiMn2O4 (LMO), cathodes present large stability when cycled in aqueous electrolytes, contrasting its behavior in conventional organic electrolytes in Lithium-ion batteries (LIBs). To elucidate the mechanisms underlying this distinctive behavior, we employ unconventional characterization techniques, including Variable Energy Positron Annihilation Lifetime Spectroscopy (VEPALS), Tip-Enhanced Raman Spectroscopy (TERS) and macro-Raman Spectroscopy (with mm-size laser spot). These still rather unexplored techniques in the battery field provide complementary information across different length scales, revealing previously hidden features. VEPALS offers atomic-scale insights, uncovering cationic defects and sub-nanometer pores that tend to collapse with cycling. TERS, operating at the nanometric range at the surface, captured the presence of Mn3O4 and its dissolution with cycling, elucidating dynamic changes during operation. Additionally, TERS highlights SO42- accumulation at grain boundaries. Macro-Raman Spectroscopy focuses on the micrometer scale, depicting small changes in the cathode's long-range order, suggesting a slow but progressive loss of crystalline quality under operation. Integrating these techniques provides a comprehensive assessment of LMO cathode stability in aqueous electrolytes, offering multifaceted insights into phase and defect evolution that can help to rationalize the origin of such stability when compared to conventional organic electrolytes. Our findings advance the understanding of LMO behavior in aqueous environments and provide guidelines for its development for next-generation LIBs.

Rejuvenation and memory, long considered the distinguishing features of spin glasses, have recently been proven to result from the growth of multiple length scales. This insight, enabled by simulations on the Janus~II supercomputer, has opened the door to a quantitative analysis. We combine numerical simulations with comparable experiments to introduce two coefficients that quantify memory. A third coefficient has been recently presented by Freedberg et al. We show that these coefficients are physically equivalent by studying their temperature and waiting-time dependence., Comment: 10 pages, 8 figures

Ion intercalation of perovskite oxides in liquid electrolytes is a very promising method for controlling their functional properties while storing charge, which opens the potential application in different energy and information technologies. Although the role of defect chemistry in the oxygen intercalation in a gaseous environment is well established, the mechanism of ion intercalation in liquid electrolytes at room temperature is poorly understood. In this study, the defect chemistry during ion intercalation of La0.5Sr0.5FeO3-{\delta} thin films in alkaline electrolytes is studied. Oxygen and proton intercalation into the LSF perovskite structure is observed at moderate electrochemical potentials (0.5 V to -0.4 V), giving rise to a change in the oxidation state of Fe (as a charge compensation mechanism). The variation of the concentration of holes as a function of the intercalation potential was characterized by in-situ ellipsometry and the concentration of electron holes was indirectly quantified for different electrochemical potentials. Finally, a dilute defect chemistry model that describes the variation of defect species during ionic intercalation was developed.

We unveil the multifractal behavior of Ising spin glasses in their low-temperature phase. Using the Janus II custom-built supercomputer, the spin-glass correlation function is studied locally. Dramatic fluctuations are found when pairs of sites at the same distance are compared. The scaling of these fluctuations, as the spin-glass coherence length grows with time, is characterized through the computation of the singularity spectrum and its corresponding Legendre transform. A comparatively small number of site pairs controls the average correlation that governs the response to a magnetic field. We explain how this scenario of dramatic fluctuations (at length scales smaller than the coherence length) can be reconciled with the smooth, self-averaging behavior that has long been considered to describe spin-glass dynamics., Comment: 15 pages, 15 figures

We have developed a novel in situ methodology for the direct study of mass transport properties in oxides with spatial and unprecedented time resolution, based on Raman spectroscopy coupled to isothermal isotope exchanges. Changes in the isotope concentration, resulting in a Raman frequency shift, can be followed in real time, not accessible by conventional methods, enabling complementary insights for the study of ion transport properties of electrode and electrolyte materials for advanced solid-state electrochemical devices. The proof of concept and strengths of isotope exchange Raman spectroscopy (IERS) are demonstrated by studying the oxygen isotope back-exchange in gadolinium-doped ceria (CGO) thin films. Resulting oxygen self-diffusion and surface exchange coefficients are compared to conventional time-of-flight secondary ion mass spectrometry (ToF-SIMS) characterisation and literature values, showing good agreement, while at the same time providing additional insight, challenging established assumptions. IERS captivates through its rapidity, simple setup, non-destructive nature, cost effectiveness and versatile fields of application and thus can readily be integrated as new standard tool for in situ and operando characterization in many laboratories worldwide. The applicability of this method is expected to consolidate our understanding of elementary physicochemical processes and impact various emerging fields including solid oxide cells, battery research and beyond.

The extended principle of superposition has been a touchstone of spin glass dynamics for almost thirty years. The Uppsala group has demonstrated its validity for the metallic spin glass, CuMn, for magnetic fields $H$ up to 10 Oe at the reduced temperature $T_\mathrm{r}=T/T_\mathrm{g} = 0.95$, where $T_\mathrm{g}$ is the spin glass condensation temperature. For $H > 10$ Oe, they observe a departure from linear response which they ascribe to the development of non-linear dynamics. The thrust of this paper is to develop a microscopic origin for this behavior by focusing on the time development of the spin glass correlation length, $\xi(t,t_\mathrm{w};H)$. Here, $t$ is the time after $H$ changes, and $t_\mathrm{w}$ is the time from the quench for $T>T_\mathrm{g}$ to the working temperature $T$ until $H$ changes. We connect the growth of $\xi(t,t_\mathrm{w};H)$ to the barrier heights $\Delta(t_\mathrm{w})$ that set the dynamics. The effect of $H$ on the magnitude of $\Delta(t_\mathrm{w})$ is responsible for affecting differently the two dynamical protocols associated with turning $H$ off (TRM, or thermoremanent magnetization) or on (ZFC, or zero field-cooled magnetization). In this paper, we display the difference between the zero-field cooled $\xi_{\text {ZFC}}(t,t_\mathrm{w};H)$ and the thermoremanent magnetization $\xi_{\text {TRM}}(t,t_\mathrm{w};H)$ correlation lengths as $H$ increases, both experimentally and through numerical simulations, corresponding to the violation of the extended principle of superposition in line with the finding of the Uppsala Group., Comment: 23 pages and 18 figures

Memory and rejuvenation effects in the magnetic response of off-equilibrium spin glasses have been widely regarded as the doorway into the experimental exploration of ultrametricity and temperature chaos (maybe the most exotic features in glassy free-energy landscapes). Unfortunately, despite more than twenty years of theoretical efforts following the experimental discovery of memory and rejuvenation, these effects have thus far been impossible to simulate reliably. Yet, three recent developments convinced us to accept this challenge: first, the custom-built Janus II supercomputer makes it possible to carry out "numerical experiments" in which the very same quantities that can be measured in single crystals of CuMn are computed from the simulation, allowing for parallel analysis of the simulation/experiment data. Second, Janus II simulations have taught us how numerical and experimental length scales should be compared. Third, we have recently understood how temperature chaos materializes in aging dynamics. All three aspects have proved crucial for reliably reproducing rejuvenation and memory effects on the computer. Our analysis shows that (at least) three different length scales play a key role in aging dynamics, while essentially all theoretical analyses of the aging dynamics emphasize the presence and the crucial role of a single glassy correlation length., Comment: Main text:13 pages, 4 figures SM: 7 pages, 8 figures

A novel non-destructive methodology for operando observation of ion intercalation and the state of charge on battery electrodes is presented based on spectroscopic ellipsometry (SE). The potentiality of this technique for performing time-resolved measurements of (de-)lithiation processes on electrode materials has been demonstrated using thin film spinel LiMn$_2$O$_4$ as a cathode for Li-ion batteries. The chemical diffusivity of Li+ ions in this material has been determined by the time evolution of the Li insertion into the film, and in the light of these results, it has been possible to relate the controversial pseudo-capacitive behavior observed in this material to the nanostructure of the layer. The results presented here establish a new way for operando characterization of battery materials and devices providing a powerful tool for the understanding of ion diffusion mechanisms in a collection of electrochemical devices., Comment: Showcasing research journal cover

Entering a new era of sustainable energy generation and consumption, new solutions for powering consumer electronics are required to tackle the limited capacity provided by the portable power sources employed nowadays. Hydrocarbon-fed micro-fuel cells represent a promising technology for this purpose, and micro-reactor technology can indeed enable their integration for portable applications. In this work, we present the design and fully scalable wafer-level fabrication of a MEMS-based catalytic micro-reactor, paving the way towards on-board hydrogen production for portable power generators. The device consists of an array of thousands of vertically-aligned micro-channels, 500 um in length and 50 um in diameter, for an overall superficial area per unit volume of 120 cm2 cm-3 and it embeds a thin-film heater for efficient reaction start-up. Functionalization of the active area was achieved by atomic layer deposition, resulting in the uniform coating of a Pt/Al2O3 heterogeneous catalyst. The temperature-dependent dimethyl ether (DME)-to-syngas conversion is tested through steam reforming (SR) and partial oxidation (POX) reactions. Here, conversion rates up to 74% and hydrogen selectivity of 60% are obtained by steam reforming at 650dC, while a specific volumetric hydrogen production of 4.5 mLH2 mL-1DME cm-3REACTOR at 600dC is obtained from DME POX in a standalone device tested by means of a 3D printed ceramic housing.

Micro-solid oxide fuel cells based on thin films have strong potential for use in portable power devices. However, devices based on silicon substrates typically involve thin-film metallic electrodes which are unstable at high temperatures. Devices based on bulk metal substrates overcome these limitations, though performance is hindered by the challenge of growing state-of-the-art epitaxial materials on metals. Here, we demonstrate for the first time the growth of epitaxial cathode materials on metal substrates (stainless steel) commercially supplied with epitaxial electrolyte layers (1.5 {um (Y2O3)0.15(ZrO2)0.85 (YSZ) + 50 nm CeO2). We create epitaxial mesoporous cathodes of (La0.60Sr0.40)0.95Co0.20Fe0.80O3 (LSCF) on the substrate by growing LSCF/MgO vertically aligned nanocomposite films by pulsed laser deposition, followed by selectively etching out the MgO. To enable valid comparison with the literature, the cathodes are also grown on single-crystal substrates, confirming state-of-the-art performance with an area specific resistance of 100ohmegacm2 at 500dC and activation energy down to 0.97 eV. The work marks an important step toward the commercialization of high-performance micro-solid oxide fuel cells for portable power applications.

Unveiling point defects concentration in transition metal oxide thin films is essential to understand and eventually control their functional properties, employed in an increasing number of applications and devices. Despite this unquestionable interest, there is a lack of available experimental techniques able to estimate the defect chemistry and equilibrium constants in such oxides at intermediate-to-low temperatures. In this study, the defect chemistry of a relevant material such as La1-xSrxFeO3-d (LSF) with (x = 0.2, 0.4 and 0.5 (LSF20, LSF40 and LSF50 respectively) is obtained by using a novel in situ spectroscopic ellipsometry approach applied to thin films. Through this technique, the concentration of holes in LSF is correlated to measured optical properties and its evolution with temperature and oxygen partial pressure is determined. In this way, a systematic description of defect chemistry in LSF thin films in the temperature range from 350dC to 500dC is obtained for the first time, which represents a step forward in the understanding of LSF20, LSF40 and LSF50 for emerging low temperature applications., Comment: Front cover of the journal

This work reports the design, manufacturing and catalytic activity characterization of a micro-reformer for hydrogen-rich gas generation integrated in portable-solid oxide fuel cells (u-SOFCs). The reformer has been designed as a silicon micro monolithic substrate compatible with the mainstream microelectronics fabrication technologies ensuring a cost-effective high reproducibility and reliability. Design and geometry of the system have been optimized comparing with the previous design, consisting in an array of more than 7x103 vertical through-silicon micro channels perfectly aligned (50 {um diameter) and a 5 W integrated serpentine heater consisting of three stacked metallic layers (TiW, W and Au) for perfect adhesion and passivation. Traditional fuels for SOFCs, such as ethanol or methanol, have been replaced by dimethyl ether (DME) and the chosen catalyst for DME conversion consists of Pd nanoparticles grafted on an alumina active support. The micro-channels have been coated by atomic layer deposition (ALD) with amorphous Al2O3 and the influence of rapid thermal processing (RTP) on such film has been studied. A customized ceramic 3D-printed holder has been designed to measure the specific hydrogen production rates, DME conversion and selectivity profiles of such catalyst at different temperatures.

The synergy between experiment, theory, and simulations enables a microscopic analysis of spin-glass dynamics in a magnetic field in the vicinity of and below the spin-glass transition temperature $T_\mathrm{g}$. The spin-glass correlation length, $\xi(t,t_\mathrm{w};T)$, is analysed both in experiments and in simulations in terms of the waiting time $t_\mathrm{w}$ after the spin glass has been cooled down to a stabilised measuring temperature $T

We find a dynamic effect in the non-equilibrium dynamics of a spin glass that closely parallels equilibrium temperature chaos. This effect, that we name dynamic temperature chaos, is spatially heterogeneous to a large degree. The key controlling quantity is the time-growing spin-glass coherence length. Our detailed characterization of dynamic temperature chaos paves the way for the analysis of recent and forthcoming experiments. This work has been made possible thanks to the most massive simulation to date of non-equilibrium dynamics, carried out on the Janus~II custom-built supercomputer., Comment: Version accepted for publication in Communication Physics 10 pages, 9 figures

The correlation length $\xi$, a key quantity in glassy dynamics, can now be precisely measured for spin glasses both in experiments and in simulations. However, known analysis methods lead to discrepancies either for large external fields or close to the glass temperature. We solve this problem by introducing a scaling law that takes into account both the magnetic field and the time-dependent spin-glass correlation length. The scaling law is successfully tested against experimental measurements in a CuMn single crystal and against large-scale simulations on the Janus II dedicated computer., Comment: Revised version, including supplemental material

Grain boundaries in Sr-doped LaMnO3 thin films have been shown to strongly influence the electronic and oxygen mass transport properties, being able to profoundly modify the nature of the material. The unique behaviour of the grain boundaries can be correlated with substantial modifications of the cation concentration at the interfaces, which can be tuned by changing the overall cationic ratio in the films. In this work, we study the electronic properties of La0.8Sr0.2Mn1-yO3 thin films with variable Mn content. The influence of the cationic composition on the grain boundary and grain bulk electronic properties is elucidated by studying the manganese valence state evolution using spectroscopy techniques and by confronting the electronic properties of epitaxial and polycrystalline films. Substantial differences in the electronic conduction mechanism are found in the presence of grain boundaries and depending on the manganese content. Moreover, the unique defect chemistry of the nanomaterial is elucidated by measuring the electrical resistance of the thin films as a function of oxygen partial pressure, disclosing the importance of the cationic local non-stoichiometry on the thin films behavior.

Interface-dominated materials such as nanocrystalline thin films have emerged as an enthralling class of materials able to engineer functional properties of transition metal oxides widely used in energy and information technologies. In particular, it has been proved that strain-induced defects in grain boundaries of manganites deeply impact their functional properties by boosting their oxygen mass transport while abating their electronic and magnetic order. In this work, the origin of these dramatic changes is correlated for the first time with strong modifications of the anionic and cationic composition in the vicinity of strained grain boundary regions. We are also able to alter the grain boundary composition by tuning the overall cationic content in the films, which represents a new and powerful tool, beyond the classical space charge layer effect, for engineering electronic and mass transport properties of metal oxide thin films useful for a collection of relevant solid state devices.

Strontium doped lanthanum manganite thin films were deposited by pulsed laser deposition on yttria stabilized zirconia single crystals for a comprehensive electrochemical characterization of the material acting as a cathode. A physically meaningful electrical model was employed to fit the electrochemical impedance spectroscopy results in order to extract the main oxygen mass transport parameters as a function of the temperature and oxygen partial pressure. The oxygen diffusion and surface exchange coefficients extracted from the analysis showed several orders of magnitude of enhancement with respect to the bulk values reported in the literature and an unexpectedly low dependence with the oxygen partial pressure. Different observations were combined to propose a mechanism for the enhanced incorporation of oxygen in interface dominated thin films mainly based on the high concentration of oxygen vacancies expected in the grain boundaries.

New emerging disciplines such as Nanoionics and Iontronics are dealing with the exploitation of mesoscopic size effects in materials, which become visible (if not predominant) when downsizing the system to the nanoscale. Driven by the worldwide standardisation of thin film deposition techniques, the access to radically different properties than those found in the bulk macroscopic systems can be accomplished. This opens up promising approaches for the development of advanced microdevices, by taking advantage of the nanostructural deviations found in nanometre sized, interface dominated materials compared to the ideal relaxed structure of the bulk. A completely new set of functionalities can be explored, with implications in many different fields such as energy conversion and storage, or information technologies. This manuscript reviews the strategies, employed and foreseen, for engineering mass transport properties in thin film ceramics, with the focus in oxide ionic and mixed ionic electronic conductors and their application in micro power sources.

There is a growing interest in the development of functional metal oxides with mixed ionic electronic conduction for their application in different strategic fields. In particular, ionic transport-related phenomena are of primary importance in energy transformation and storage devices, such as solid oxides fuel cells (SOFC) or batteries. Traditionally, the main issue that hindered a wide spread of solid state electrochemical devices is the slow ionic conduction and the reduced number of families of materials with pure ionic conductivity able to play the role of the electrolyte. This is due to the low mobility of the ionic charge carriers and their mechanical and coulombic interactions with their host crystal structures. For relatively big and charged ions, e.g. O2- or Na+, the low ionic conductivity imposes operating at high temperatures, which is incompatible with some applications such as portable electronics.

Ion transport in solid-state devices is of great interest for current and future energy and information technologies. A superior enhancement of several orders of magnitude of the oxygen diffusivity has been recently reported for grain boundaries in lanthanum strontium manganites. However, the significance and extent of this unique phenomenon is not yet established. Here, we fabricate a thin film continuous composition map of the La0.8Sr0.2(Mn1-xCox)0.85O3-d family revealing a substantial enhancement of the grain boundary oxygen mass transport properties for the entire range of compositions. Through isotope-exchange depth profiling coupled to secondary ion mass spectroscopy, we show that this excellent performance is not directly linked to the bulk of the material but to the intrinsic nature of the grain boundary. In particular, the great increase of the oxygen diffusion in Mn-rich compositions unveils an unprecedented catalytic performance in the field of Mixed Ionic Electronic Conductors. These results present grain boundaries engineering as a novel strategy for designing highly performing materials for solid state ionics based devices.

The Mpemba effect occurs when a hot system cools faster than an initially colder one, when both are refrigerated in the same thermal reservoir. Using the custom built supercomputer Janus II, we study the Mpemba effect in spin glasses and show that it is a non-equilibrium process, governed by the coherence length \xi of the system. The effect occurs when the bath temperature lies in the glassy phase, but it is not necessary for the thermal protocol to cross the critical temperature. In fact, the Mpemba effect follows from a strong relationship between the internal energy and \xi that turns out to be a sure-tell sign of being in the glassy phase. Thus, the Mpemba effect presents itself as an intriguing new avenue for the experimental study of the coherence length in supercooled liquids and other glass formers., Comment: Version accepted for publication in PNAS. 6 pages, 7 figures

Experiments on spin glasses can now make precise measurements of the exponent $z(T)$ governing the growth of glassy domains, while our computational capabilities allow us to make quantitative predictions for experimental scales. However, experimental and numerical values for $z(T)$ have differed. We use new simulations on the Janus II computer to resolve this discrepancy, finding a time-dependent $z(T, t_w)$, which leads to the experimental value through mild extrapolations. Furthermore, theoretical insight is gained by studying a crossover between the $T = T_c$ and $T = 0$ fixed points., Comment: Version accepted for publication in PRL. 12 pages, 9 figures

Portable electronic devices are already an indispensable part of our daily life; and their increasing number and demand for higher performance is becoming a challenge for the research community. In particular, a major concern is the way to efficiently power these energy-demanding devices, assuring long grid independency with high efficiency, sustainability and cheap production. In this context, technologies beyond Li-ion are receiving increasing attention, among which the development of micro solid oxide fuel cells ({\mu}SOFC) stands out. In particular, {\mu}SOFC provides a high energy density, high efficiency and opens the possibility to the use of different fuels, such as hydrocarbons. Yet, its high operating temperature has typically hindered its application as miniaturized portable device. Recent advances have however set a completely new range of lower operating temperatures, i.e. 350-450C, as compared to the typical >900C needed for classical bulk SOFC systems. In this work, a comprehensive review of the status of the technology is presented. The main achievements, as well as the most important challenges still pending are discussed, regarding (i.) the cell design and microfabrication, and (ii.) the integration of functional electrolyte and electrode materials. To conclude, the different strategies foreseen for a wide deployment of the technology as new portable power source are underlined.

Additive manufacturing represents a revolution due to its unique capabilities for freeform fabrication of near net shapes with strong reduction of waste material and capital cost. These unfair advantages are especially relevant for expensive and energy-demanding manufacturing processes of advanced ceramics such as Yttria-stabilized Zirconia, the state-of-the-art electrolyte in Solid Oxide Fuel Cell applications. In this study, self-supported electrolytes of yttria-stabilized zirconia have been printed by using a stereolithography three-dimensional printer. Printed electrolytes and complete cells fabricated with cathode and anode layers of lanthanum strontium manganite- and nickel oxide-yttria-stabilized zirconia composites, respectively, were electrochemical characterized showing full functionality. In addition, more complex configurations of the electrolyte have been printed yielding an increase of the performance entirely based on geometrical aspects. Complementary, a numerical model has been developed and validated as predictive tool for designing more advanced configurations that will enable highly performing and fully customized devices in the next future

We first reproduce on the Janus and Janus II computers a milestone experiment that measures the spin-glass coherence length through the lowering of free-energy barriers induced by the Zeeman effect. Secondly we determine the scaling behavior that allows a quantitative analysis of a new experiment reported in the companion Letter [S. Guchhait and R. Orbach, Phys. Rev. Lett. 118, 157203 (2017)]. The value of the coherence length estimated through the analysis of microscopic correlation functions turns out to be quantitatively consistent with its measurement through macroscopic response functions. Further, non-linear susceptibilities, recently measured in glass-forming liquids, scale as powers of the same microscopic length., Comment: 6 pages, 4 figures

The unifying feature of glass formers (such as polymers, supercooled liquids, colloids, granulars, spin glasses, superconductors, ...) is a sluggish dynamics at low temperatures. Indeed, their dynamics is so slow that thermal equilibrium is never reached in macroscopic samples: in analogy with living beings, glasses are said to age. Here, we show how to relate experimentally relevant quantities with the experimentally unreachable low-temperature equilibrium phase. We have performed a very accurate computation of the non-equilibrium fluctuation-dissipation ratio for the three-dimensional Edwards-Anderson Ising spin glass, by means of large-scale simulations on the special-purpose computers Janus and Janus II. This ratio (computed for finite times on very large, effectively infinite, systems) is compared with the equilibrium probability distribution of the spin overlap for finite sizes. The resulting quantitative statics-dynamics dictionary, based on observables that can be measured with current experimental methods, could allow the experimental exploration of important features of the spin-glass phase without uncontrollable extrapolations to infinite times or system sizes., Comment: Version accepted for publication in PNAS. Reuploaded to fix a typo in the list of authors

We perform equilibrium parallel-tempering simulations of the 3D Ising Edwards-Anderson spin glass in a field. A traditional analysis shows no signs of a phase transition. Yet, we encounter dramatic fluctuations in the behaviour of the model: Averages over all the data only describe the behaviour of a small fraction of it. Therefore we develop a new approach to study the equilibrium behaviour of the system, by classifying the measurements as a function of a conditioning variate. We propose a finite-size scaling analysis based on the probability distribution function of the conditioning variate, which may accelerate the convergence to the thermodynamic limit. In this way, we find a non-trivial spectrum of behaviours, where a part of the measurements behaves as the average, while the majority of them shows signs of scale invariance. As a result, we can estimate the temperature interval where the phase transition in a field ought to lie, if it exists. Although this would-be critical regime is unreachable with present resources, the numerical challenge is finally well posed., Comment: 42 pages, 19 figures. Minor changes and added figure (results unchanged)

We report a high-precision finite-size scaling study of the critical behavior of the three-dimensional Ising Edwards-Anderson model (the Ising spin glass). We have thermalized lattices up to L=40 using the Janus dedicated computer. Our analysis takes into account leading-order corrections to scaling. We obtain Tc = 1.1019(29) for the critical temperature, \nu = 2.562(42) for the thermal exponent, \eta = -0.3900(36) for the anomalous dimension and \omega = 1.12(10) for the exponent of the leading corrections to scaling. Standard (hyper)scaling relations yield \alpha = -5.69(13), \beta = 0.782(10) and \gamma = 6.13(11). We also compute several universal quantities at Tc., Comment: 9 pages, 5 figures

This paper describes the architecture, the development and the implementation of Janus II, a new generation application-driven number cruncher optimized for Monte Carlo simulations of spin systems (mainly spin glasses). This domain of computational physics is a recognized grand challenge of high-performance computing: the resources necessary to study in detail theoretical models that can make contact with experimental data are by far beyond those available using commodity computer systems. On the other hand, several specific features of the associated algorithms suggest that unconventional computer architectures, which can be implemented with available electronics technologies, may lead to order of magnitude increases in performance, reducing to acceptable values on human scales the time needed to carry out simulation campaigns that would take centuries on commercially available machines. Janus II is one such machine, recently developed and commissioned, that builds upon and improves on the successful JANUS machine, which has been used for physics since 2008 and is still in operation today. This paper describes in detail the motivations behind the project, the computational requirements, the architecture and the implementation of this new machine and compares its expected performances with those of currently available commercial systems., Comment: 28 pages, 6 figures

We study the off-equilibrium dynamics of the three-dimensional Ising spin glass in the presence of an external magnetic field. We have performed simulations both at fixed temperature and with an annealing protocol. Thanks to the Janus special-purpose computer, based on FPGAs, we have been able to reach times equivalent to 0.01 seconds in experiments. We have studied the system relaxation both for high and for low temperatures, clearly identifying a dynamical transition point. This dynamical temperature is strictly positive and depends on the external applied magnetic field. We discuss different possibilities for the underlying physics, which include a thermodynamical spin-glass transition, a mode-coupling crossover or an interpretation reminiscent of the random first-order picture of structural glasses., Comment: Version accepted in PRE. 23 pages, 23 figures

We describe Janus, a massively parallel FPGA-based computer optimized for the simulation of spin glasses, theoretical models for the behavior of glassy materials. FPGAs (as compared to GPUs or many-core processors) provide a complementary approach to massively parallel computing. In particular, our model problem is formulated in terms of binary variables, and floating-point operations can be (almost) completely avoided. The FPGA architecture allows us to run many independent threads with almost no latencies in memory access, thus updating up to 1024 spins per cycle. We describe Janus in detail and we summarize the physics results obtained in four years of operation of this machine; we discuss two types of physics applications: long simulations on very large systems (which try to mimic and provide understanding about the experimental non-equilibrium dynamics), and low-temperature equilibrium simulations using an artificial parallel tempering dynamics. The time scale of our non-equilibrium simulations spans eleven orders of magnitude (from picoseconds to a tenth of a second). On the other hand, our equilibrium simulations are unprecedented both because of the low temperatures reached and for the large systems that we have brought to equilibrium. A finite-time scaling ansatz emerges from the detailed comparison of the two sets of simulations. Janus has made it possible to perform spin-glass simulations that would take several decades on more conventional architectures. The paper ends with an assessment of the potential of possible future versions of the Janus architecture, based on state-of-the-art technology., Comment: 19 pages, 3 figures

Spin glasses are a longstanding model for the sluggish dynamics that appears at the glass transition. However, spin glasses differ from structural glasses for a crucial feature: they enjoy a time reversal symmetry. This symmetry can be broken by applying an external magnetic field, but embarrassingly little is known about the critical behaviour of a spin glass in a field. In this context, the space dimension is crucial. Simulations are easier to interpret in a large number of dimensions, but one must work below the upper critical dimension (i.e., in d<6) in order for results to have relevance for experiments. Here we show conclusive evidence for the presence of a phase transition in a four-dimensional spin glass in a field. Two ingredients were crucial for this achievement: massive numerical simulations were carried out on the Janus special-purpose computer, and a new and powerful finite-size scaling method., Comment: 10 pages, 6 figures

We study the sample-to-sample fluctuations of the overlap probability densities from large-scale equilibrium simulations of the three-dimensional Edwards-Anderson spin glass below the critical temperature. Ultrametricity, Stochastic Stability and Overlap Equivalence impose constraints on the moments of the overlap probability densities that can be tested against numerical data. We found small deviations from the Ghirlanda-Guerra predictions, which get smaller as system size increases. We also focus on the shape of the overlap distribution, comparing the numerical data to a mean-field-like prediction in which finite-size effects are taken into account by substituting delta functions with broad peaks, Comment: 11 pages, 8 figures; version accepted for publication on PRB and updated Refs

The number of people using online social networks in their everyday life is continuously growing at a pace never saw before. This new kind of communication has an enormous impact on opinions, cultural trends, information spreading and even in the commercial success of new products. More importantly, social online networks have revealed as a fundamental organizing mechanism in recent country-wide social movements. In this paper, we provide a quantitative analysis of the structural and dynamical patterns emerging from the activity of an online social network around the ongoing May 15th (15M) movement in Spain. Our network is made up by users that exchanged tweets in a time period of one month, which includes the birth and stabilization of the 15M movement. We characterize in depth the growth of such dynamical network and find that it is scale-free with communities at the mesoscale. We also find that its dynamics exhibits typical features of critical systems such as robustness and power-law distributions for several quantities. Remarkably, we report that the patterns characterizing the spreading dynamics are asymmetric, giving rise to a clear distinction between information sources and sinks. Our study represent a first step towards the use of data from online social media to comprehend modern societal dynamics., Comment: 16 pages, 7 figures

We numerically study the aging properties of the dynamical heterogeneities in the Ising spin glass. We find that a phase transition takes place during the aging process. Statics-dynamics correspondence implies that systems of finite size in equilibrium have static heterogeneities that obey Finite-Size Scaling, thus signaling an analogous phase transition in the thermodynamical limit. We compute the critical exponents and the transition point in the equilibrium setting, and use them to show that aging in dynamic heterogeneities can be described by a Finite-Time Scaling Ansatz, with potential implications for experimental work., Comment: Version accepted for publication in Phys. Rev. Lett. 4 pages, 3 postscript figures

We present a massive equilibrium simulation of the three-dimensional Ising spin glass at low temperatures. The Janus special-purpose computer has allowed us to equilibrate, using parallel tempering, L=32 lattices down to T=0.64 Tc. We demonstrate the relevance of equilibrium finite-size simulations to understand experimental non-equilibrium spin glasses in the thermodynamical limit by establishing a time-length dictionary. We conclude that non-equilibrium experiments performed on a time scale of one hour can be matched with equilibrium results on L=110 lattices. A detailed investigation of the probability distribution functions of the spin and link overlap, as well as of their correlation functions, shows that Replica Symmetry Breaking is the appropriate theoretical framework for the physically relevant length scales. Besides, we improve over existing methodologies to ensure equilibration in parallel tempering simulations., Comment: 48 pages, 19 postscript figures, 9 tables. Version accepted for publication in the Journal of Statistical Mechanics

We study the 3D Disordered Potts Model with p=5 and p=6. Our numerical simulations (that severely slow down for increasing p) detect a very clear spin glass phase transition. We evaluate the critical exponents and the critical value of the temperature, and we use known results at lower $p$ values to discuss how they evolve for increasing p. We do not find any sign of the presence of a transition to a ferromagnetic regime., Comment: 9 pages and 9 Postscript figures. Final version published in J. Stat. Mech

We have investigated the phase transition in the Heisenberg spin glass using massive numerical simulations to study larger sizes, 48x48x48, than have been attempted before at a spin glass phase transition. A finite-size scaling analysis indicates that the data is compatible with the most economical scenario: a common transition temperature for spins and chiralities., Comment: Version to appear in Phys. Rev. B

We perform numerical simulations, including parallel tempering, on the Potts glass model with binary random quenched couplings using the JANUS application-oriented computer. We find and characterize a glassy transition, estimating the location of the transition and the value of the critical exponents. We show that there is no ferromagnetic transition in a large temperature range around the glassy critical temperature. We also compare our results with those obtained recently on the "random permutation" Potts glass., Comment: 7 pages and 3 figures. Corrected minor typos

Using the dedicated computer Janus, we follow the nonequilibrium dynamics of the Ising spin glass in three dimensions for eleven orders of magnitude. The use of integral estimators for the coherence and correlation lengths allows us to study dynamic heterogeneities and the presence of a replicon mode and to obtain safe bounds on the Edwards-Anderson order parameter below the critical temperature. We obtain good agreement with experimental determinations of the temperature-dependent decay exponents for the thermoremanent magnetization. This magnitude is observed to scale with the much harder to measure coherence length, a potentially useful result for experimentalists. The exponents for energy relaxation display a linear dependence on temperature and reasonable extrapolations to the critical point. We conclude examining the time growth of the coherence length, with a comparison of critical and activated dynamics., Comment: 38 pages, 26 figures

We study numerically the nonequilibrium dynamics of the Ising Spin Glass, for a time that spans eleven orders of magnitude, thus approaching the experimentally relevant scale (i.e. {\em seconds}). We introduce novel analysis techniques that allow to compute the coherence length in a model-independent way. Besides, we present strong evidence for a replicon correlator and for overlap equivalence. The emerging picture is compatible with non-coarsening behavior., Comment: 4 pages, 4 eps color figures. Version accepted for publication in Phys. Rev. Lett

The author states an exact expression of the distribution of primes., Comment: 2 pages

We investigate the critical properties of the four-state commutative random permutation glassy Potts model in three and four dimensions by means of Monte Carlo simulation and of a finite size scaling analysis. Thanks to the use of a field programmable gate array we have been able to thermalize a large number of samples of systems with large volume. This has allowed us to observe a spin-glass ordered phase in d=4 and to study the critical properties of the transition. In d=3, our results are consistent with the presence of a Kosterlitz-Thouless transition, but we cannot exclude transient effects due to a value of the lower critical dimension slightly below 3., Comment: 9 pages, 8 Postscript figures

This paper describes JANUS, a modular massively parallel and reconfigurable FPGA-based computing system. Each JANUS module has a computational core and a host. The computational core is a 4x4 array of FPGA-based processing elements with nearest-neighbor data links. Processors are also directly connected to an I/O node attached to the JANUS host, a conventional PC. JANUS is tailored for, but not limited to, the requirements of a class of hard scientific applications characterized by regular code structure, unconventional data manipulation instructions and not too large data-base size. We discuss the architecture of this configurable machine, and focus on its use on Monte Carlo simulations of statistical mechanics. On this class of application JANUS achieves impressive performances: in some cases one JANUS processing element outperfoms high-end PCs by a factor ~ 1000. We also discuss the role of JANUS on other classes of scientific applications., Comment: 11 pages, 6 figures. Improved version, largely rewritten, submitted to Computing in Science & Engineering