@article{10178, abstract = {In dense biological tissues, cell types performing different roles remain segregated by maintaining sharp interfaces. To better understand the mechanisms for such sharp compartmentalization, we study the effect of an imposed heterotypic tension at the interface between two distinct cell types in a fully 3D Voronoi model for confluent tissues. We find that cells rapidly sort and self-organize to generate a tissue-scale interface between cell types, and cells adjacent to this interface exhibit signature geometric features including nematic-like ordering, bimodal facet areas, and registration, or alignment, of cell centers on either side of the two-tissue interface. The magnitude of these features scales directly with the magnitude of the imposed tension, suggesting that biologists can estimate the magnitude of tissue surface tension between two tissue types simply by segmenting a 3D tissue. To uncover the underlying physical mechanisms driving these geometric features, we develop two minimal, ordered models using two different underlying lattices that identify an energetic competition between bulk cell shapes and tissue interface area. When the interface area dominates, changes to neighbor topology are costly and occur less frequently, which generates the observed geometric features.}, author = {Sahu, Preeti and Schwarz, J. M. and Manning, M. Lisa}, issn = {1367-2630}, journal = {New Journal of Physics}, number = {9}, publisher = {IOP Publishing}, title = {{Geometric signatures of tissue surface tension in a three-dimensional model of confluent tissue}}, doi = {10.1088/1367-2630/ac23f1}, volume = {23}, year = {2021}, } @article{9909, abstract = {Roots are composed of different root types and, in the dicotyledonous Arabidopsis, typically consist of a primary root that branches into lateral roots. Adventitious roots emerge from non-root tissue and are formed upon wounding or other types of abiotic stress. Here, we investigated adventitious root (AR) formation in Arabidopsis hypocotyls under conditions of altered abscisic acid (ABA) signaling. Exogenously applied ABA suppressed AR formation at 0.25 µM or higher doses. AR formation was less sensitive to the synthetic ABA analog pyrabactin (PB). However, PB was a more potent inhibitor at concentrations above 1 µM, suggesting that it was more selective in triggering a root inhibition response. Analysis of a series of phosphonamide and phosphonate pyrabactin analogs suggested that adventitious root formation and lateral root branching are differentially regulated by ABA signaling. ABA biosynthesis and signaling mutants affirmed a general inhibitory role of ABA and point to PYL1 and PYL2 as candidate ABA receptors that regulate AR inhibition.}, author = {Zeng, Yinwei and Verstraeten, Inge and Trinh, Hoang Khai and Heugebaert, Thomas and Stevens, Christian V. and Garcia-Maquilon, Irene and Rodriguez, Pedro L. and Vanneste, Steffen and Geelen, Danny}, issn = {2073-4425}, journal = {Genes}, number = {8}, publisher = {MDPI}, title = {{Arabidopsis hypocotyl adventitious root formation is suppressed by ABA signaling}}, doi = {10.3390/genes12081141}, volume = {12}, year = {2021}, } @article{9333, abstract = {We revise a previous result about the Fröhlich dynamics in the strong coupling limit obtained in Griesemer (Rev Math Phys 29(10):1750030, 2017). In the latter it was shown that the Fröhlich time evolution applied to the initial state φ0⊗ξα, where φ0 is the electron ground state of the Pekar energy functional and ξα the associated coherent state of the phonons, can be approximated by a global phase for times small compared to α2. In the present note we prove that a similar approximation holds for t=O(α2) if one includes a nontrivial effective dynamics for the phonons that is generated by an operator proportional to α−2 and quadratic in creation and annihilation operators. Our result implies that the electron ground state remains close to its initial state for times of order α2, while the phonon fluctuations around the coherent state ξα can be described by a time-dependent Bogoliubov transformation.}, author = {Mitrouskas, David Johannes}, issn = {1573-0530}, journal = {Letters in Mathematical Physics}, publisher = {Springer Nature}, title = {{A note on the Fröhlich dynamics in the strong coupling limit}}, doi = {10.1007/s11005-021-01380-7}, volume = {111}, year = {2021}, } @article{9244, abstract = {Organ function depends on tissues adopting the correct architecture. However, insights into organ architecture are currently hampered by an absence of standardized quantitative 3D analysis. We aimed to develop a robust technology to visualize, digitalize, and segment the architecture of two tubular systems in 3D: double resin casting micro computed tomography (DUCT). As proof of principle, we applied DUCT to a mouse model for Alagille syndrome (Jag1Ndr/Ndr mice), characterized by intrahepatic bile duct paucity, that can spontaneously generate a biliary system in adulthood. DUCT identified increased central biliary branching and peripheral bile duct tortuosity as two compensatory processes occurring in distinct regions of Jag1Ndr/Ndr liver, leading to full reconstitution of wild-type biliary volume and phenotypic recovery. DUCT is thus a powerful new technology for 3D analysis, which can reveal novel phenotypes and provide a standardized method of defining liver architecture in mouse models.}, author = {Hankeova, Simona and Salplachta, Jakub and Zikmund, Tomas and Kavkova, Michaela and Van Hul, Noémi and Brinek, Adam and Smekalova, Veronika and Laznovsky, Jakub and Dawit, Feven and Jaros, Josef and Bryja, Vítězslav and Lendahl, Urban and Ellis, Ewa and Nemeth, Antal and Fischler, Björn and Hannezo, Edouard B and Kaiser, Jozef and Andersson, Emma Rachel}, issn = {2050-084X}, journal = {eLife}, publisher = {eLife Sciences Publications}, title = {{DUCT reveals architectural mechanisms contributing to bile duct recovery in a mouse model for alagille syndrome}}, doi = {10.7554/eLife.60916}, volume = {10}, year = {2021}, } @article{9306, abstract = {Assemblies of actin and its regulators underlie the dynamic morphology of all eukaryotic cells. To understand how actin regulatory proteins work together to generate actin-rich structures such as filopodia, we analyzed the localization of diverse actin regulators within filopodia in Drosophila embryos and in a complementary in vitro system of filopodia-like structures (FLSs). We found that the composition of the regulatory protein complex where actin is incorporated (the filopodial tip complex) is remarkably heterogeneous both in vivo and in vitro. Our data reveal that different pairs of proteins correlate with each other and with actin bundle length, suggesting the presence of functional subcomplexes. This is consistent with a theoretical framework where three or more redundant subcomplexes join the tip complex stochastically, with any two being sufficient to drive filopodia formation. We provide an explanation for the observed heterogeneity and suggest that a mechanism based on multiple components allows stereotypical filopodial dynamics to arise from diverse upstream signaling pathways.}, author = {Dobramysl, Ulrich and Jarsch, Iris Katharina and Inoue, Yoshiko and Shimo, Hanae and Richier, Benjamin and Gadsby, Jonathan R. and Mason, Julia and Szałapak, Alicja and Ioannou, Pantelis Savvas and Correia, Guilherme Pereira and Walrant, Astrid and Butler, Richard and Hannezo, Edouard B and Simons, Benjamin D. and Gallop, Jennifer L.}, issn = {1540-8140}, journal = {Journal of Cell Biology}, number = {4}, publisher = {Rockefeller University Press}, title = {{Stochastic combinations of actin regulatory proteins are sufficient to drive filopodia formation}}, doi = {10.1083/jcb.202003052}, volume = {220}, year = {2021}, } @article{9205, abstract = {Cryo-EM grid preparation is an important bottleneck in protein structure determination, especially for membrane proteins, typically requiring screening of a large number of conditions. We systematically investigated the effects of buffer components, blotting conditions and grid types on the outcome of grid preparation of five different membrane protein samples. Aggregation was the most common type of problem which was addressed by changing detergents, salt concentration or reconstitution of proteins into nanodiscs or amphipols. We show that the optimal concentration of detergent is between 0.05 and 0.4% and that the presence of a low concentration of detergent with a high critical micellar concentration protects the proteins from denaturation at the air-water interface. Furthermore, we discuss the strategies for achieving an adequate ice thickness, particle coverage and orientation distribution on free ice and on support films. Our findings provide a clear roadmap for comprehensive screening of conditions for cryo-EM grid preparation of membrane proteins.}, author = {Kampjut, Domen and Steiner, Julia and Sazanov, Leonid A}, issn = {2589-0042}, journal = {iScience}, number = {3}, publisher = {Elsevier}, title = {{Cryo-EM grid optimization for membrane proteins}}, doi = {10.1016/j.isci.2021.102139}, volume = {24}, year = {2021}, } @article{9334, abstract = {Polaritons with directional in-plane propagation and ultralow losses in van der Waals (vdW) crystals promise unprecedented manipulation of light at the nanoscale. However, these polaritons present a crucial limitation: their directional propagation is intrinsically determined by the crystal structure of the host material, imposing forbidden directions of propagation. Here, we demonstrate that directional polaritons (in-plane hyperbolic phonon polaritons) in a vdW crystal (α-phase molybdenum trioxide) can be directed along forbidden directions by inducing an optical topological transition, which emerges when the slab is placed on a substrate with a given negative permittivity (4H–silicon carbide). By visualizing the transition in real space, we observe exotic polaritonic states between mutually orthogonal hyperbolic regimes, which unveil the topological origin of the transition: a gap opening in the dispersion. This work provides insights into optical topological transitions in vdW crystals, which introduce a route to direct light at the nanoscale.}, author = {Duan, J. and Álvarez-Pérez, G. and Voronin, K. V. and Prieto Gonzalez, Ivan and Taboada-Gutiérrez, J. and Volkov, V. S. and Martín-Sánchez, J. and Nikitin, A. Y. and Alonso-González, P.}, issn = {2375-2548}, journal = {Science Advances}, number = {14}, publisher = {AAAS}, title = {{Enabling propagation of anisotropic polaritons along forbidden directions via a topological transition}}, doi = {10.1126/sciadv.abf2690}, volume = {7}, year = {2021}, } @article{9412, abstract = {We extend our recent result [22] on the central limit theorem for the linear eigenvalue statistics of non-Hermitian matrices X with independent, identically distributed complex entries to the real symmetry class. We find that the expectation and variance substantially differ from their complex counterparts, reflecting (i) the special spectral symmetry of real matrices onto the real axis; and (ii) the fact that real i.i.d. matrices have many real eigenvalues. Our result generalizes the previously known special cases where either the test function is analytic [49] or the first four moments of the matrix elements match the real Gaussian [59, 44]. The key element of the proof is the analysis of several weakly dependent Dyson Brownian motions (DBMs). The conceptual novelty of the real case compared with [22] is that the correlation structure of the stochastic differentials in each individual DBM is non-trivial, potentially even jeopardising its well-posedness.}, author = {Cipolloni, Giorgio and Erdös, László and Schröder, Dominik J}, issn = {1083-6489}, journal = {Electronic Journal of Probability}, publisher = {Institute of Mathematical Statistics}, title = {{Fluctuation around the circular law for random matrices with real entries}}, doi = {10.1214/21-EJP591}, volume = {26}, year = {2021}, } @article{9361, abstract = {The multimeric matrix (M) protein of clinically relevant paramyxoviruses orchestrates assembly and budding activity of viral particles at the plasma membrane (PM). We identified within the canine distemper virus (CDV) M protein two microdomains, potentially assuming α-helix structures, which are essential for membrane budding activity. Remarkably, while two rationally designed microdomain M mutants (E89R, microdomain 1 and L239D, microdomain 2) preserved proper folding, dimerization, interaction with the nucleocapsid protein, localization at and deformation of the PM, the virus-like particle formation, as well as production of infectious virions (as monitored using a membrane budding-complementation system), were, in sharp contrast, strongly impaired. Of major importance, raster image correlation spectroscopy (RICS) revealed that both microdomains contributed to finely tune M protein mobility specifically at the PM. Collectively, our data highlighted the cornerstone membrane budding-priming activity of two spatially discrete M microdomains, potentially by coordinating the assembly of productive higher oligomers at the PM.}, author = {Gast, Matthieu and Kadzioch, Nicole P. and Milius, Doreen and Origgi, Francesco and Plattet, Philippe}, issn = {2379-5042}, journal = {mSphere}, number = {2}, publisher = {American Society for Microbiology}, title = {{Oligomerization and cell egress controlled by two microdomains of canine distemper virus matrix protein}}, doi = {10.1128/mSphere.01024-20}, volume = {6}, year = {2021}, } @article{9362, abstract = {A central goal in systems neuroscience is to understand the functions performed by neural circuits. Previous top-down models addressed this question by comparing the behaviour of an ideal model circuit, optimised to perform a given function, with neural recordings. However, this requires guessing in advance what function is being performed, which may not be possible for many neural systems. To address this, we propose an inverse reinforcement learning (RL) framework for inferring the function performed by a neural network from data. We assume that the responses of each neuron in a network are optimised so as to drive the network towards ‘rewarded’ states, that are desirable for performing a given function. We then show how one can use inverse RL to infer the reward function optimised by the network from observing its responses. This inferred reward function can be used to predict how the neural network should adapt its dynamics to perform the same function when the external environment or network structure changes. This could lead to theoretical predictions about how neural network dynamics adapt to deal with cell death and/or varying sensory stimulus statistics.}, author = {Chalk, Matthew J and Tkačik, Gašper and Marre, Olivier}, issn = {1932-6203}, journal = {PLoS ONE}, number = {4}, publisher = {Public Library of Science}, title = {{Inferring the function performed by a recurrent neural network}}, doi = {10.1371/journal.pone.0248940}, volume = {16}, year = {2021}, } @article{9392, abstract = {Humans conceptualize the diversity of life by classifying individuals into types we call ‘species’1. The species we recognize influence political and financial decisions and guide our understanding of how units of diversity evolve and interact. Although the idea of species may seem intuitive, a debate about the best way to define them has raged even before Darwin2. So much energy has been devoted to the so-called ‘species problem’ that no amount of discourse will ever likely solve it2,3. Dozens of species concepts are currently recognized3, but we lack a concrete understanding of how much researchers actually disagree and the factors that cause them to think differently1,2. To address this, we used a survey to quantify the species problem for the first time. The results indicate that the disagreement is extensive: two randomly chosen respondents will most likely disagree on the nature of species. The probability of disagreement is not predicted by researcher experience or broad study system, but tended to be lower among researchers with similar focus, training and who study the same organism. Should we see this diversity of perspectives as a problem? We argue that we should not.}, author = {Stankowski, Sean and Ravinet, Mark}, issn = {1879-0445}, journal = {Current Biology}, number = {9}, pages = {R428--R429}, publisher = {Cell Press}, title = {{Quantifying the use of species concepts}}, doi = {10.1016/j.cub.2021.03.060}, volume = {31}, year = {2021}, } @article{9679, abstract = {The relative motion of three impenetrable particles on a ring, in our case two identical fermions and one impurity, is isomorphic to a triangular quantum billiard. Depending on the ratio κ of the impurity and fermion masses, the billiards can be integrable or non-integrable (also referred to in the main text as chaotic). To set the stage, we first investigate the energy level distributions of the billiards as a function of 1/κ ∈ [0, 1] and find no evidence of integrable cases beyond the limiting values 1/κ = 1 and 1/κ = 0. Then, we use machine learning tools to analyze properties of probability distributions of individual quantum states. We find that convolutional neural networks can correctly classify integrable and non-integrable states. The decisive features of the wave functions are the normalization and a large number of zero elements, corresponding to the existence of a nodal line. The network achieves typical accuracies of 97%, suggesting that machine learning tools can be used to analyze and classify the morphology of probability densities obtained in theory or experiment.}, author = {Huber, David and Marchukov, Oleksandr V. and Hammer, Hans Werner and Volosniev, Artem}, issn = {1367-2630}, journal = {New Journal of Physics}, number = {6}, publisher = {IOP Publishing}, title = {{Morphology of three-body quantum states from machine learning}}, doi = {10.1088/1367-2630/ac0576}, volume = {23}, year = {2021}, } @article{9410, abstract = {Antibiotic concentrations vary dramatically in the body and the environment. Hence, understanding the dynamics of resistance evolution along antibiotic concentration gradients is critical for predicting and slowing the emergence and spread of resistance. While it has been shown that increasing the concentration of an antibiotic slows resistance evolution, how adaptation to one antibiotic concentration correlates with fitness at other points along the gradient has not received much attention. Here, we selected populations of Escherichia coli at several points along a concentration gradient for three different antibiotics, asking how rapidly resistance evolved and whether populations became specialized to the antibiotic concentration they were selected on. Populations selected at higher concentrations evolved resistance more slowly but exhibited equal or higher fitness across the whole gradient. Populations selected at lower concentrations evolved resistance rapidly, but overall fitness in the presence of antibiotics was lower. However, these populations readily adapted to higher concentrations upon subsequent selection. Our results indicate that resistance management strategies must account not only for the rates of resistance evolution but also for the fitness of evolved strains.}, author = {Lagator, Mato and Uecker, Hildegard and Neve, Paul}, issn = {1744-957X}, journal = {Biology letters}, number = {5}, publisher = {Royal Society of London}, title = {{Adaptation at different points along antibiotic concentration gradients}}, doi = {10.1098/rsbl.2020.0913}, volume = {17}, year = {2021}, } @article{8601, abstract = {We consider large non-Hermitian real or complex random matrices X with independent, identically distributed centred entries. We prove that their local eigenvalue statistics near the spectral edge, the unit circle, coincide with those of the Ginibre ensemble, i.e. when the matrix elements of X are Gaussian. This result is the non-Hermitian counterpart of the universality of the Tracy–Widom distribution at the spectral edges of the Wigner ensemble.}, author = {Cipolloni, Giorgio and Erdös, László and Schröder, Dominik J}, issn = {1432-2064}, journal = {Probability Theory and Related Fields}, publisher = {Springer Nature}, title = {{Edge universality for non-Hermitian random matrices}}, doi = {10.1007/s00440-020-01003-7}, year = {2021}, } @article{9318, abstract = {We consider a system of N bosons in the mean-field scaling regime for a class of interactions including the repulsive Coulomb potential. We derive an asymptotic expansion of the low-energy eigenstates and the corresponding energies, which provides corrections to Bogoliubov theory to any order in 1/N.}, author = {Bossmann, Lea and Petrat, Sören P and Seiringer, Robert}, issn = {2050-5094}, journal = {Forum of Mathematics, Sigma}, publisher = {Cambridge University Press}, title = {{Asymptotic expansion of low-energy excitations for weakly interacting bosons}}, doi = {10.1017/fms.2021.22}, volume = {9}, year = {2021}, } @article{8198, abstract = {We investigate how the critical driving amplitude at the Floquet many-body localized (MBL) to ergodic phase transition differs between smooth and nonsmooth drives. To this end, we numerically study a disordered spin-1/2 chain which is periodically driven by a sine or square-wave drive over a wide range of driving frequencies. In both cases the critical driving amplitude increases monotonically with the frequency, and at large frequencies it is identical for the two drives. However, at low and intermediate frequencies the critical amplitude of the square-wave drive depends strongly on the frequency, while that of the sinusoidal drive is almost constant over a wide frequency range. By analyzing the density of drive-induced resonances we conclude that this difference is due to resonances induced by the higher harmonics which are present (absent) in the Fourier spectrum of the square-wave (sine) drive. Furthermore, we suggest a numerically efficient method for estimating the frequency dependence of the critical driving amplitudes for different drives which is based on calculating the density of drive-induced resonances. We conclude that delocalization occurs once the density of drive-induced resonances reaches a critical value determined only by the static system.}, author = {Diringer, Asaf A. and Gulden, Tobias}, issn = {2469-9969}, journal = {Physical Review B}, number = {21}, publisher = {American Physical Society}, title = {{Impact of drive harmonics on the stability of Floquet many-body localization}}, doi = {10.1103/PhysRevB.103.214204}, volume = {103}, year = {2021}, } @article{9603, abstract = {Mosaic analysis with double markers (MADM) offers one approach to visualize and concomitantly manipulate genetically defined cells in mice with single-cell resolution. MADM applications include the analysis of lineage, single-cell morphology and physiology, genomic imprinting phenotypes, and dissection of cell-autonomous gene functions in vivo in health and disease. Yet, MADM can only be applied to <25% of all mouse genes on select chromosomes to date. To overcome this limitation, we generate transgenic mice with knocked-in MADM cassettes near the centromeres of all 19 autosomes and validate their use across organs. With this resource, >96% of the entire mouse genome can now be subjected to single-cell genetic mosaic analysis. Beyond a proof of principle, we apply our MADM library to systematically trace sister chromatid segregation in distinct mitotic cell lineages. We find striking chromosome-specific biases in segregation patterns, reflecting a putative mechanism for the asymmetric segregation of genetic determinants in somatic stem cell division.}, author = {Contreras, Ximena and Amberg, Nicole and Davaatseren, Amarbayasgalan and Hansen, Andi H and Sonntag, Johanna and Andersen, Lill and Bernthaler, Tina and Streicher, Carmen and Heger, Anna-Magdalena and Johnson, Randy L. and Schwarz, Lindsay A. and Luo, Liqun and Rülicke, Thomas and Hippenmeyer, Simon}, issn = {2211-1247}, journal = {Cell Reports}, number = {12}, publisher = {Cell Press}, title = {{A genome-wide library of MADM mice for single-cell genetic mosaic analysis}}, doi = {10.1016/j.celrep.2021.109274}, volume = {35}, year = {2021}, } @article{9363, abstract = {Optogenetics has been harnessed to shed new mechanistic light on current and future therapeutic strategies. This has been to date achieved by the regulation of ion flow and electrical signals in neuronal cells and neural circuits that are known to be affected by disease. In contrast, the optogenetic delivery of trophic biochemical signals, which support cell survival and are implicated in degenerative disorders, has never been demonstrated in an animal model of disease. Here, we reengineered the human and Drosophila melanogaster REarranged during Transfection (hRET and dRET) receptors to be activated by light, creating one-component optogenetic tools termed Opto-hRET and Opto-dRET. Upon blue light stimulation, these receptors robustly induced the MAPK/ERK proliferative signaling pathway in cultured cells. In PINK1B9 flies that exhibit loss of PTEN-induced putative kinase 1 (PINK1), a kinase associated with familial Parkinson’s disease (PD), light activation of Opto-dRET suppressed mitochondrial defects, tissue degeneration and behavioral deficits. In human cells with PINK1 loss-of-function, mitochondrial fragmentation was rescued using Opto-dRET via the PI3K/NF-кB pathway. Our results demonstrate that a light-activated receptor can ameliorate disease hallmarks in a genetic model of PD. The optogenetic delivery of trophic signals is cell type-specific and reversible and thus has the potential to inspire novel strategies towards a spatio-temporal regulation of tissue repair.}, author = {Inglés Prieto, Álvaro and Furthmann, Nikolas and Crossman, Samuel H. and Tichy, Alexandra Madelaine and Hoyer, Nina and Petersen, Meike and Zheden, Vanessa and Bicher, Julia and Gschaider-Reichhart, Eva and György, Attila and Siekhaus, Daria E and Soba, Peter and Winklhofer, Konstanze F. and Janovjak, Harald L}, issn = {1553-7404}, journal = {PLoS genetics}, number = {4}, pages = {e1009479}, publisher = {Public Library of Science}, title = {{Optogenetic delivery of trophic signals in a genetic model of Parkinson's disease}}, doi = {10.1371/journal.pgen.1009479}, volume = {17}, year = {2021}, } @article{9640, abstract = {Selection and random drift determine the probability that novel mutations fixate in a population. Population structure is known to affect the dynamics of the evolutionary process. Amplifiers of selection are population structures that increase the fixation probability of beneficial mutants compared to well-mixed populations. Over the past 15 years, extensive research has produced remarkable structures called strong amplifiers which guarantee that every beneficial mutation fixates with high probability. But strong amplification has come at the cost of considerably delaying the fixation event, which can slow down the overall rate of evolution. However, the precise relationship between fixation probability and time has remained elusive. Here we characterize the slowdown effect of strong amplification. First, we prove that all strong amplifiers must delay the fixation event at least to some extent. Second, we construct strong amplifiers that delay the fixation event only marginally as compared to the well-mixed populations. Our results thus establish a tight relationship between fixation probability and time: Strong amplification always comes at a cost of a slowdown, but more than a marginal slowdown is not needed.}, author = {Tkadlec, Josef and Pavlogiannis, Andreas and Chatterjee, Krishnendu and Nowak, Martin A.}, issn = {2041-1723}, journal = {Nature Communications}, number = {1}, publisher = {Springer Nature}, title = {{Fast and strong amplifiers of natural selection}}, doi = {10.1038/s41467-021-24271-w}, volume = {12}, year = {2021}, } @article{10401, abstract = {Theoretical and experimental studies of the interaction between spins and temperature are vital for the development of spin caloritronics, as they dictate the design of future devices. In this work, we propose a two-terminal cold-atom simulator to study that interaction. The proposed quantum simulator consists of strongly interacting atoms that occupy two temperature reservoirs connected by a one-dimensional link. First, we argue that the dynamics in the link can be described using an inhomogeneous Heisenberg spin chain whose couplings are defined by the local temperature. Second, we show the existence of a spin current in a system with a temperature difference by studying the dynamics that follows the spin-flip of an atom in the link. A temperature gradient accelerates the impurity in one direction more than in the other, leading to an overall spin current similar to the spin Seebeck effect.}, author = {Barfknecht, Rafael E. and Foerster, Angela and Zinner, Nikolaj T. and Volosniev, Artem}, issn = {2399-3650}, journal = {Communications Physics}, number = {1}, publisher = {Springer Nature}, title = {{Generation of spin currents by a temperature gradient in a two-terminal device}}, doi = {10.1038/s42005-021-00753-7}, volume = {4}, year = {2021}, }