@article{912, abstract = {We consider a many-body system of fermionic atoms interacting via a local pair potential and subject to an external potential within the framework of Bardeen-Cooper-Schrieffer (BCS) theory. We measure the free energy of the whole sample with respect to the free energy of a reference state which allows us to define a BCS functional with boundary conditions at infinity. Our main result is a lower bound for this energy functional in terms of expressions that typically appear in Ginzburg-Landau functionals. }, author = {Deuchert, Andreas}, issn = {0022-2488}, journal = {Journal of Mathematical Physics}, number = {8}, publisher = {AIP Publishing}, title = {{A lower bound for the BCS functional with boundary conditions at infinity}}, doi = {10.1063/1.4996580}, volume = {58}, year = {2017}, } @article{744, abstract = {In evolutionary game theory interactions between individuals are often assumed obligatory. However, in many real-life situations, individuals can decide to opt out of an interaction depending on the information they have about the opponent. We consider a simple evolutionary game theoretic model to study such a scenario, where at each encounter between two individuals the type of the opponent (cooperator/defector) is known with some probability, and where each individual either accepts or opts out of the interaction. If the type of the opponent is unknown, a trustful individual accepts the interaction, whereas a suspicious individual opts out of the interaction. If either of the two individuals opt out both individuals remain without an interaction. We show that in the prisoners dilemma optional interactions along with suspicious behaviour facilitates the emergence of trustful cooperation.}, author = {Priklopil, Tadeas and Chatterjee, Krishnendu and Nowak, Martin}, issn = {0022-5193}, journal = {Journal of Theoretical Biology}, pages = {64 -- 72}, publisher = {Elsevier}, title = {{Optional interactions and suspicious behaviour facilitates trustful cooperation in prisoners dilemma}}, doi = {10.1016/j.jtbi.2017.08.025}, volume = {433}, year = {2017}, } @article{1076, abstract = {Signatures of the Coulomb corrections in the photoelectron momentum distribution during laser-induced ionization of atoms or ions in tunneling and multiphoton regimes are investigated analytically in the case of a one-dimensional problem. A high-order Coulomb-corrected strong-field approximation is applied, where the exact continuum state in the S matrix is approximated by the eikonal Coulomb-Volkov state including the second-order corrections to the eikonal. Although without high-order corrections our theory coincides with the known analytical R-matrix (ARM) theory, we propose a simplified procedure for the matrix element derivation. Rather than matching the eikonal Coulomb-Volkov wave function with the bound state as in the ARM theory to remove the Coulomb singularity, we calculate the matrix element via the saddle-point integration method by time as well as by coordinate, and in this way avoiding the Coulomb singularity. The momentum shift in the photoelectron momentum distribution with respect to the ARM theory due to high-order corrections is analyzed for tunneling and multiphoton regimes. The relation of the quantum corrections to the tunneling delay time is discussed.}, author = {Klaiber, Michael and Daněk, Jiří and Yakaboylu, Enderalp and Hatsagortsyan, Karen and Keitel, Christoph}, issn = {2469-9926}, journal = {Physical Review A}, number = {2}, publisher = {American Physical Society}, title = {{Strong-field ionization via a high-order Coulomb-corrected strong-field approximation}}, doi = {10.1103/PhysRevA.95.023403}, volume = {95}, year = {2017}, } @inproceedings{6526, abstract = {This paper studies the complexity of estimating Rényi divergences of discrete distributions: p observed from samples and the baseline distribution q known a priori. Extending the results of Acharya et al. (SODA'15) on estimating Rényi entropy, we present improved estimation techniques together with upper and lower bounds on the sample complexity. We show that, contrarily to estimating Rényi entropy where a sublinear (in the alphabet size) number of samples suffices, the sample complexity is heavily dependent on events occurring unlikely in q, and is unbounded in general (no matter what an estimation technique is used). For any divergence of integer order bigger than 1, we provide upper and lower bounds on the number of samples dependent on probabilities of p and q (the lower bounds hold for non-integer orders as well). We conclude that the worst-case sample complexity is polynomial in the alphabet size if and only if the probabilities of q are non-negligible. This gives theoretical insights into heuristics used in the applied literature to handle numerical instability, which occurs for small probabilities of q. Our result shows that they should be handled with care not only because of numerical issues, but also because of a blow up in the sample complexity.}, author = {Skórski, Maciej}, booktitle = {2017 IEEE International Symposium on Information Theory}, isbn = {9781509040964}, location = {Aachen, Germany}, publisher = {IEEE}, title = {{On the complexity of estimating Rènyi divergences}}, doi = {10.1109/isit.2017.8006529}, year = {2017}, } @article{467, abstract = {Recently there has been a significant effort to handle quantitative properties in formal verification and synthesis. While weighted automata over finite and infinite words provide a natural and flexible framework to express quantitative properties, perhaps surprisingly, some basic system properties such as average response time cannot be expressed using weighted automata or in any other known decidable formalism. In this work, we introduce nested weighted automata as a natural extension of weighted automata, which makes it possible to express important quantitative properties such as average response time. In nested weighted automata, a master automaton spins off and collects results from weighted slave automata, each of which computes a quantity along a finite portion of an infinite word. Nested weighted automata can be viewed as the quantitative analogue of monitor automata, which are used in runtime verification. We establish an almost-complete decidability picture for the basic decision problems about nested weighted automata and illustrate their applicability in several domains. In particular, nested weighted automata can be used to decide average response time properties.}, author = {Chatterjee, Krishnendu and Henzinger, Thomas A and Otop, Jan}, issn = {1529-3785}, journal = {ACM Transactions on Computational Logic}, number = {4}, publisher = {ACM}, title = {{Nested weighted automata}}, doi = {10.1145/3152769}, volume = {18}, year = {2017}, } @phdthesis{6287, abstract = {The main objects considered in the present work are simplicial and CW-complexes with vertices forming a random point cloud. In particular, we consider a Poisson point process in R^n and study Delaunay and Voronoi complexes of the first and higher orders and weighted Delaunay complexes obtained as sections of Delaunay complexes, as well as the Čech complex. Further, we examine theDelaunay complex of a Poisson point process on the sphere S^n, as well as of a uniform point cloud, which is equivalent to the convex hull, providing a connection to the theory of random polytopes. Each of the complexes in question can be endowed with a radius function, which maps its cells to the radii of appropriately chosen circumspheres, called the radius of the cell. Applying and developing discrete Morse theory for these functions, joining it together with probabilistic and sometimes analytic machinery, and developing several integral geometric tools, we aim at getting the distributions of circumradii of typical cells. For all considered complexes, we are able to generalize and obtain up to constants the distribution of radii of typical intervals of all types. In low dimensions the constants can be computed explicitly, thus providing the explicit expressions for the expected numbers of cells. In particular, it allows to find the expected density of simplices of every dimension for a Poisson point process in R^4, whereas the result for R^3 was known already in 1970's.}, author = {Nikitenko, Anton}, issn = {2663-337X}, pages = {86}, publisher = {Institute of Science and Technology Austria}, title = {{Discrete Morse theory for random complexes }}, doi = {10.15479/AT:ISTA:th_873}, year = {2017}, } @phdthesis{1127, abstract = {Plant hormone auxin and its transport between cells belong to the most important mechanisms controlling plant development. Auxin itself could change localization of PINs and thereby control direction of its own flow. We performed an expression profiling experiment in Arabidopsis roots to identify potential regulators of PIN polarity which are transcriptionally regulated by auxin signalling. We identified several novel regulators and performed a detailed characterization of the transcription factor WRKY23 (At2g47260) and its role in auxin feedback on PIN polarity. Gain-of-function and dominant-negative mutants revealed that WRKY23 plays a crucial role in mediating the auxin effect on PIN polarity. In concordance, typical polar auxin transport processes such as gravitropism and leaf vascular pattern formation were disturbed by interfering with WRKY23 function. In order to identify direct targets of WRKY23, we performed consequential expression profiling experiments using a WRKY23 inducible gain-of-function line and dominant-negative WRKY23 line that is defunct in PIN re-arrangement. Among several genes mostly related to the groups of cell wall and defense process regulators, we identified LYSINE-HISTIDINE TRANSPORTER 1 (LHT1; At5g40780), a small amino acid permease gene from the amino acid/auxin permease family (AAAP), we present its detailed characterisation in auxin feedback on PIN repolarization, identified its transcriptional regulation, we propose a potential mechanism of its action. Moreover, we identified also a member of receptor-like protein kinase LRR-RLK (LEUCINE-RICH REPEAT TRANSMEMBRANE PROTEIN KINASE PROTEIN 1; LRRK1; At1g05700), which also affects auxin-dependent PIN re-arrangement. We described its transcriptional behaviour, subcellular localization. Based on global expression data, we tried to identify ligand responsible for mechanism of signalling and suggest signalling partner and interactors. Additionally, we described role of novel phytohormone group, strigolactone, in auxin-dependent PIN re-arrangement, that could be a fundament for future studies in this field. Our results provide first insights into an auxin transcriptional network targeting PIN localization and thus regulating plant development. We highlighted WRKY23 transcriptional network and characterised its mediatory role in plant development. We identified direct effectors of this network, LHT1 and LRRK1, and describe their roles in PIN re-arrangement and PIN-dependent auxin transport processes.}, author = {Prat, Tomas}, issn = {2663-337X}, pages = {131}, publisher = {Institute of Science and Technology Austria}, title = {{Identification of novel regulators of PIN polarity and development of novel auxin sensor}}, year = {2017}, } @phdthesis{839, abstract = {This thesis describes a brittle fracture simulation method for visual effects applications. Building upon a symmetric Galerkin boundary element method, we first compute stress intensity factors following the theory of linear elastic fracture mechanics. We then use these stress intensities to simulate the motion of a propagating crack front at a significantly higher resolution than the overall deformation of the breaking object. Allowing for spatial variations of the material's toughness during crack propagation produces visually realistic, highly-detailed fracture surfaces. Furthermore, we introduce approximations for stress intensities and crack opening displacements, resulting in both practical speed-up and theoretically superior runtime complexity compared to previous methods. While we choose a quasi-static approach to fracture mechanics, ignoring dynamic deformations, we also couple our fracture simulation framework to a standard rigid-body dynamics solver, enabling visual effects artists to simulate both large scale motion, as well as fracturing due to collision forces in a combined system. As fractures inside of an object grow, their geometry must be represented both in the coarse boundary element mesh, as well as at the desired fine output resolution. Using a boundary element method, we avoid complicated volumetric meshing operations. Instead we describe a simple set of surface meshing operations that allow us to progressively add cracks to the mesh of an object and still re-use all previously computed entries of the linear boundary element system matrix. On the high resolution level, we opt for an implicit surface representation. We then describe how to capture fracture surfaces during crack propagation, as well as separate the individual fragments resulting from the fracture process, based on this implicit representation. We show results obtained with our method, either solving the full boundary element system in every time step, or alternatively using our fast approximations. These results demonstrate that both of these methods perform well in basic test cases and produce realistic fracture surfaces. Furthermore we show that our fast approximations substantially out-perform the standard approach in more demanding scenarios. Finally, these two methods naturally combine, using the full solution while the problem size is manageably small and switching to the fast approximations later on. The resulting hybrid method gives the user a direct way to choose between speed and accuracy of the simulation. }, author = {Hahn, David}, issn = {2663-337X}, pages = {124}, publisher = {Institute of Science and Technology Austria}, title = {{Brittle fracture simulation with boundary elements for computer graphics}}, doi = {10.15479/AT:ISTA:th_855}, year = {2017}, } @misc{5568, abstract = {Includes source codes, test cases, and example data used in the thesis Brittle Fracture Simulation with Boundary Elements for Computer Graphics. Also includes pre-built binaries of the HyENA library, but not sources - please contact the HyENA authors to obtain these sources if required (https://mech.tugraz.at/hyena)}, author = {Hahn, David}, keywords = {Boundary elements, brittle fracture, computer graphics, fracture simulation}, publisher = {Institute of Science and Technology Austria}, title = {{Source codes: Brittle fracture simulation with boundary elements for computer graphics}}, doi = {10.15479/AT:ISTA:73}, year = {2017}, } @inproceedings{559, abstract = {Proofs of space (PoS) were suggested as more ecological and economical alternative to proofs of work, which are currently used in blockchain designs like Bitcoin. The existing PoS are based on rather sophisticated graph pebbling lower bounds. Much simpler and in several aspects more efficient schemes based on inverting random functions have been suggested, but they don’t give meaningful security guarantees due to existing time-memory trade-offs. In particular, Hellman showed that any permutation over a domain of size N can be inverted in time T by an algorithm that is given S bits of auxiliary information whenever (Formula presented). For functions Hellman gives a weaker attack with S2· T≈ N2 (e.g., S= T≈ N2/3). To prove lower bounds, one considers an adversary who has access to an oracle f: [ N] → [N] and can make T oracle queries. The best known lower bound is S· T∈ Ω(N) and holds for random functions and permutations. We construct functions that provably require more time and/or space to invert. Specifically, for any constant k we construct a function [N] → [N] that cannot be inverted unless Sk· T∈ Ω(Nk) (in particular, S= T≈ (Formula presented). Our construction does not contradict Hellman’s time-memory trade-off, because it cannot be efficiently evaluated in forward direction. However, its entire function table can be computed in time quasilinear in N, which is sufficient for the PoS application. Our simplest construction is built from a random function oracle g: [N] × [N] → [ N] and a random permutation oracle f: [N] → N] and is defined as h(x) = g(x, x′) where f(x) = π(f(x′)) with π being any involution without a fixed point, e.g. flipping all the bits. For this function we prove that any adversary who gets S bits of auxiliary information, makes at most T oracle queries, and inverts h on an ϵ fraction of outputs must satisfy S2· T∈ Ω(ϵ2N2).}, author = {Abusalah, Hamza M and Alwen, Joel F and Cohen, Bram and Khilko, Danylo and Pietrzak, Krzysztof Z and Reyzin, Leonid}, isbn = {978-331970696-2}, location = {Hong Kong, China}, pages = {357 -- 379}, publisher = {Springer}, title = {{Beyond Hellman’s time-memory trade-offs with applications to proofs of space}}, doi = {10.1007/978-3-319-70697-9_13}, volume = {10625}, year = {2017}, } @article{1144, abstract = {We show that matrix elements of functions of N × N Wigner matrices fluctuate on a scale of order N−1/2 and we identify the limiting fluctuation. Our result holds for any function f of the matrix that has bounded variation thus considerably relaxing the regularity requirement imposed in [7, 11].}, author = {Erdös, László and Schröder, Dominik J}, journal = {Electronic Communications in Probability}, publisher = {Institute of Mathematical Statistics}, title = {{Fluctuations of functions of Wigner matrices}}, doi = {10.1214/16-ECP38}, volume = {21}, year = {2017}, } @article{795, abstract = {We introduce a common generalization of the strong Hanani–Tutte theorem and the weak Hanani–Tutte theorem: if a graph G has a drawing D in the plane where every pair of independent edges crosses an even number of times, then G has a planar drawing preserving the rotation of each vertex whose incident edges cross each other evenly in D. The theorem is implicit in the proof of the strong Hanani–Tutte theorem by Pelsmajer, Schaefer and Štefankovič. We give a new, somewhat simpler proof.}, author = {Fulek, Radoslav and Kynčl, Jan and Pálvölgyi, Dömötör}, issn = {1077-8926}, journal = {The Electronic Journal of Combinatorics}, number = {3}, publisher = {Electronic Journal of Combinatorics}, title = {{Unified Hanani Tutte theorem}}, doi = {10.37236/6663}, volume = {24}, year = {2017}, } @article{914, abstract = {Infections with potentially lethal pathogens may negatively affect an individual’s lifespan and decrease its reproductive value. The terminal investment hypothesis predicts that individuals faced with a reduced survival should invest more into reproduction instead of maintenance and growth. Several studies suggest that individuals are indeed able to estimate their body condition and to increase their reproductive effort with approaching death, while other studies gave ambiguous results. We investigate whether queens of a perennial social insect (ant) are able to boost their reproduction following infection with an obligate killing pathogen. Social insect queens are special with regard to reproduction and aging, as they outlive conspecific non-reproductive workers. Moreover, in the ant Cardiocondyla obscurior, fecundity increases with queen age. However, it remained unclear whether this reflects negative reproductive senescence or terminal investment in response to approaching death. Here, we test whether queens of C. obscurior react to infection with the entomopathogenic fungus Metarhizium brunneum by an increased egg-laying rate. We show that a fungal infection triggers a reinforced investment in reproduction in queens. This adjustment of the reproductive rate by ant queens is consistent with predictions of the terminal investment hypothesis and is reported for the first time in a social insect.}, author = {Giehr, Julia and Grasse, Anna V and Cremer, Sylvia and Heinze, Jürgen and Schrempf, Alexandra}, issn = {2054-5703}, journal = {Royal Society Open Science}, number = {7}, publisher = {Royal Society}, title = {{Ant queens increase their reproductive efforts after pathogen infection}}, doi = {10.1098/rsos.170547}, volume = {4}, year = {2017}, } @misc{9853, abstract = {Egg laying rates and infection loads of C. obscurior queens}, author = {Giehr, Julia and Grasse, Anna V and Cremer, Sylvia and Heinze, Jürgen and Schrempf, Alexandra}, publisher = {The Royal Society}, title = {{Raw data from ant queens increase their reproductive efforts after pathogen infection}}, doi = {10.6084/m9.figshare.5117788.v1}, year = {2017}, } @article{953, abstract = {The role of natural selection in the evolution of adaptive phenotypes has undergone constant probing by evolutionary biologists, employing both theoretical and empirical approaches. As Darwin noted, natural selection can act together with other processes, including random changes in the frequencies of phenotypic differences that are not under strong selection, and changes in the environment, which may reflect evolutionary changes in the organisms themselves. As understanding of genetics developed after 1900, the new genetic discoveries were incorporated into evolutionary biology. The resulting general principles were summarized by Julian Huxley in his 1942 book Evolution: the modern synthesis. Here, we examine how recent advances in genetics, developmental biology and molecular biology, including epigenetics, relate to today's understanding of the evolution of adaptations. We illustrate how careful genetic studies have repeatedly shown that apparently puzzling results in a wide diversity of organisms involve processes that are consistent with neo-Darwinism. They do not support important roles in adaptation for processes such as directed mutation or the inheritance of acquired characters, and therefore no radical revision of our understanding of the mechanism of adaptive evolution is needed.}, author = {Charlesworth, Deborah and Barton, Nicholas H and Charlesworth, Brian}, journal = {Proceedings of the Royal Society of London Series B Biological Sciences}, number = {1855}, publisher = {Royal Society}, title = {{The sources of adaptive evolution}}, doi = {10.1098/rspb.2016.2864}, volume = {284}, year = {2017}, } @article{557, abstract = {PURPOSE. Gene therapy of retinal ganglion cells (RGCs) has promise as a powerful therapeutic for the rescue and regeneration of these cells after optic nerve damage. However, early after damage, RGCs undergo atrophic changes, including gene silencing. It is not known if these changes will deleteriously affect transduction and transgene expression, or if the therapeutic protein can influence reactivation of the endogenous genome. METHODS. Double-transgenic mice carrying a Rosa26-(LoxP)-tdTomato reporter, and a mutant allele for the proapoptotic Bax gene were reared. The Bax mutant blocks apoptosis, but RGCs still exhibit nuclear atrophy and gene silencing. At times ranging from 1 hour to 4 weeks after optic nerve crush (ONC), eyes received an intravitreal injection of AAV2 virus carrying the Cre recombinase. Successful transduction was monitored by expression of the tdTomato reporter. Immunostaining was used to localize tdTomato expression in select cell types. RESULTS. Successful transduction of RGCs was achieved at all time points after ONC using AAV2 expressing Cre from the phosphoglycerate kinase (Pgk) promoter, but not the CMV promoter. ONC promoted an increase in the transduction of cell types in the inner nuclear layer, including Müller cells and rod bipolar neurons. There was minimal evidence of transduction of amacrine cells and astrocytes in the inner retina or optic nerve. CONCLUSIONS. Damaged RGCs can be transduced and at least some endogenous genes can be subsequently activated. Optic nerve damage may change retinal architecture to allow greater penetration of an AAV2 virus to transduce several additional cell types in the inner nuclear layer.}, author = {Nickells, Robert and Schmitt, Heather and Maes, Margaret E and Schlamp, Cassandra}, issn = {0146-0404}, journal = {Investigative Ophthalmology and Visual Science}, number = {14}, pages = {6091 -- 6104}, publisher = {Association for Research in Vision & Ophthalmology}, title = {{AAV2 mediated transduction of the mouse retina after optic nerve injury}}, doi = {10.1167/iovs.17-22634}, volume = {58}, year = {2017}, } @article{1078, abstract = {One of the key questions in understanding plant development is how single cells behave in a larger context of the tissue. Therefore, it requires the observation of the whole organ with a high spatial- as well as temporal resolution over prolonged periods of time, which may cause photo-toxic effects. This protocol shows a plant sample preparation method for light-sheet microscopy, which is characterized by mounting the plant vertically on the surface of a gel. The plant is mounted in such a way that the roots are submerged in a liquid medium while the leaves remain in the air. In order to ensure photosynthetic activity of the plant, a custom-made lighting system illuminates the leaves. To keep the roots in darkness the water surface is covered with sheets of black plastic foil. This method allows long-term imaging of plant organ development in standardized conditions. }, author = {Von Wangenheim, Daniel and Hauschild, Robert and Friml, Jirí}, journal = {Journal of Visualized Experiments}, number = {119}, publisher = {MyJove Corporation}, title = {{Light sheet fluorescence microscopy of plant roots growing on the surface of a gel}}, doi = {10.3791/55044}, volume = {2017}, year = {2017}, } @misc{5565, abstract = {One of the key questions in understanding plant development is how single cells behave in a larger context of the tissue. Therefore, it requires the observation of the whole organ with a high spatial- as well as temporal resolution over prolonged periods of time, which may cause photo-toxic effects. This protocol shows a plant sample preparation method for light-sheet microscopy, which is characterized by mounting the plant vertically on the surface of a gel. The plant is mounted in such a way that the roots are submerged in a liquid medium while the leaves remain in the air. In order to ensure photosynthetic activity of the plant, a custom-made lighting system illuminates the leaves. To keep the roots in darkness the water surface is covered with sheets of black plastic foil. This method allows long-term imaging of plant organ development in standardized conditions. The Video is licensed under a CC BY NC ND license. }, author = {Von Wangenheim, Daniel and Hauschild, Robert and Friml, Jirí}, publisher = {Institute of Science and Technology Austria}, title = {{Light Sheet Fluorescence microscopy of plant roots growing on the surface of a gel}}, doi = {10.15479/AT:ISTA:66}, year = {2017}, } @article{825, abstract = {What data is needed about data? Describing the process to answer this question for the institutional data repository IST DataRep.}, author = {Petritsch, Barbara}, issn = {1022-2588}, journal = {Mitteilungen der Vereinigung Österreichischer Bibliothekarinnen und Bibliothekare}, number = {2}, pages = {200 -- 207}, publisher = {Vereinigung Österreichischer Bibliothekarinnen und Bibliothekare}, title = {{Metadata for research data in practice}}, doi = {10.31263/voebm.v70i2.1678}, volume = {70}, year = {2017}, } @article{1162, abstract = {Selected universal experimental properties of high-temperature superconducting (HTS) cuprates have been singled out in the last decade. One of the pivotal challenges in this field is the designation of a consistent interpretation framework within which we can describe quantitatively the universal features of those systems. Here we analyze in a detailed manner the principal experimental data and compare them quantitatively with the approach based on a single-band model of strongly correlated electrons supplemented with strong antiferromagnetic (super)exchange interaction (the so-called t−J−U model). The model rationale is provided by estimating its microscopic parameters on the basis of the three-band approach for the Cu-O plane. We use our original full Gutzwiller wave-function solution by going beyond the renormalized mean-field theory (RMFT) in a systematic manner. Our approach reproduces very well the observed hole doping (δ) dependence of the kinetic-energy gain in the superconducting phase, one of the principal non-Bardeen-Cooper-Schrieffer features of the cuprates. The calculated Fermi velocity in the nodal direction is practically δ-independent and its universal value agrees very well with that determined experimentally. Also, a weak doping dependence of the Fermi wave vector leads to an almost constant value of the effective mass in a pure superconducting phase which is both observed in experiment and reproduced within our approach. An assessment of the currently used models (t−J, Hubbard) is carried out and the results of the canonical RMFT as a zeroth-order solution are provided for comparison to illustrate the necessity of the introduced higher-order contributions.}, author = {Spałek, Jozef and Zegrodnik, Michał and Kaczmarczyk, Jan}, issn = {2469-9950}, journal = {Physical Review B}, number = {2}, publisher = {American Physical Society}, title = {{Universal properties of high temperature superconductors from real space pairing t-J-U model and its quantitative comparison with experiment}}, doi = {10.1103/PhysRevB.95.024506}, volume = {95}, year = {2017}, }