@inproceedings{18241,
  abstract     = {Multiple-input multiple-output (MIMO) radar is one of the leading depth sensing modalities. However, the usage of multiple receive channels lead to relative high costs and prevent the penetration of MIMOs in many areas such as the automotive industry. Over the last years, few studies concentrated on designing reduced measurement schemes and image reconstruction schemes for MIMO radars, however these problems have been so far addressed separately. On the other hand, recent works in optical computational imaging have demonstrated growing success of simultaneous learning-based design of the acquisition and reconstruction schemes, manifesting significant improvement in the reconstruction quality. Inspired by these successes, in this work, we propose to learn MIMO acquisition parameters in the form of receive (Rx) antenna elements locations jointly with an image neural-network based reconstruction. To this end, we propose an algorithm for training the combined acquisition-reconstruction pipeline end-to-end in a differentiable way. We demonstrate the significance of using our learned acquisition parameters with and without the neural-network reconstruction. Code and datasets will be released upon publication.},
  author       = {Weiss, Tomer and Peretz, Nissim and Vedula, Sanketh and Feuer, Arie and Bronstein, Alexander},
  booktitle    = {31st International Workshop on Machine Learning for Signal Processing},
  location     = {Gold Coast, Australia},
  publisher    = {Institute of Electrical and Electronics Engineers},
  title        = {{Joint optimization of system design and reconstruction in MIMO radar imaging}},
  doi          = {10.1109/mlsp52302.2021.9596168},
  volume       = {4},
  year         = {2021},
}

@inbook{18242,
  abstract     = {Fiber tractography is an important tool of computational neuroscience that enables reconstructing the spatial connectivity and organization of white matter of the brain. Fiber tractography takes advantage of diffusion Magnetic Resonance Imaging (dMRI) which allows measuring the apparent diffusivity of cerebral water along different spatial directions. Unfortunately, collecting such data comes at the price of reduced spatial resolution and substantially elevated acquisition times, which limits the clinical applicability of dMRI. This problem has been thus far addressed using two principal strategies. Most of the efforts have been extended towards improving the quality of signal estimation for any, yet fixed sampling scheme (defined through the choice of diffusion-encoding gradients). On the other hand, optimization over the sampling scheme has also proven to be effective. Inspired by the previous results, the present work consolidates the above strategies into a unified estimation framework, in which the optimization is carried out with respect to both estimation model and sampling design concurrently. The proposed solution offers substantial improvements in the quality of signal estimation as well as the accuracy of ensuing analysis by means of fiber tractography. While proving the optimality of the learned estimation models would probably need more extensive evaluation, we nevertheless claim that the learned sampling schemes can be of immediate use, offering a way to improve the dMRI analysis without the necessity of deploying the neural network used for their estimation. We present a comprehensive comparative analysis based on the Human Connectome Project data. Code and learned sampling designs available at https://github.com/tomer196/Learned_dMRI.},
  author       = {Weiss, Tomer and Vedula, Sanketh and Senouf, Ortal and Michailovich, Oleg and Bronstein, Alexander},
  booktitle    = {Computational Diffusion MRI},
  editor       = {Gyori, Noemi and Hutter, Jana and Nath, Vishwesh and Palombo, Marco and Pizzolato, Marco and Zhang, Fan},
  isbn         = {9783030730178},
  issn         = {1612-3786},
  location     = {Lima, Peru/Virtual},
  pages        = {13--28},
  publisher    = {Springer Nature},
  title        = {{Towards learned optimal q-space sampling in diffusion MRI}},
  doi          = {10.1007/978-3-030-73018-5_2},
  year         = {2021},
}

@inproceedings{18243,
  abstract     = {What is the best way to match the nodes of two graphs? This graph alignment problem generalizes graph isomorphism and arises in applications from social network analysis to bioinformatics. Existing solutions either require auxiliary information such as node attributes, or provide a single-scale view of the graph by translating the problem into aligning node embeddings.

In this paper, we transfer the shape-analysis concept of functional maps from the continuous to the discrete case, and treat the graph alignment problem as a special case of the problem of finding a mapping between functions on graphs. We present GRASP, a method that captures multiscale structural characteristics from the eigenvectors of the graph’s Laplacian and uses this information to align two graphs.Our experimental study, featuring noise levels higher than anything used in previous studies, shows that GRASP outperforms state-of-the-art methods for graph alignment across noise levels and graph types.},
  author       = {Hermanns, Judith and Tsitsulin, Anton and Munkhoeva, Marina and Bronstein, Alexander and Mottin, Davide and Karras, Panagiotis},
  booktitle    = {International Joint Conference on Asia-Paciﬁc Web and Web-Age Information Management},
  isbn         = {9783030858957},
  issn         = {0302-9743},
  location     = {Guangzhou, China},
  number       = {Part I},
  pages        = {44 -- 52},
  publisher    = {Springer Nature},
  title        = {{GRASP: Graph alignment through spectral signatures}},
  doi          = {10.1007/978-3-030-85896-4_4},
  volume       = {12858},
  year         = {2021},
}

@inproceedings{18244,
  abstract     = {Some face recognition methods are designed to utilize geometric information extracted from depth sensors to overcome the weaknesses of single-image based recognition technologies. However, the accurate acquisition of the depth profile is an expensive and challenging process. Here, we introduce a novel method that learns to recognize faces from stereo camera systems without the need to explicitly compute the facial surface or depth map. The raw face stereo images along with the location in the image from which the face is extracted allow the proposed CNN to improve the recognition task while avoiding the need to explicitly handle the geometric structure of the face. This way, we keep the simplicity and cost efficiency of identity authentication from a single image, while enjoying the benefits of geometric data without explicitly reconstructing it. We demonstrate that the suggested method outperforms both existing single-image and explicit depth based methods on largescale benchmarks, and even capable of recognize spoofing attacks. We also provide an ablation study that shows that the suggested method uses the face locations in the left and right images to encode informative features that improve the overall performance.},
  author       = {Livne, Amir and Aviv, Ziv and Grofit, Shahaf and Bronstein, Alexander and Kimmel, Ron},
  booktitle    = {2020 International Conference on 3D Vision (3DV)},
  isbn         = {9781728181295},
  issn         = {2475-7888},
  location     = {Fukuoka, Japan},
  publisher    = {IEEE},
  title        = {{Do we need depth in state-uf-the-art face authentication?}},
  doi          = {10.1109/3dv50981.2020.00099},
  year         = {2021},
}

@article{7685,
  abstract     = {We consider a gas of interacting bosons trapped in a box of side length one in the Gross–Pitaevskii limit. We review the proof of the validity of Bogoliubov’s prediction for the ground state energy and the low-energy excitation spectrum. This note is based on joint work with C. Brennecke, S. Cenatiempo and B. Schlein.},
  author       = {Boccato, Chiara},
  issn         = {0129-055X},
  journal      = {Reviews in Mathematical Physics},
  number       = {1},
  publisher    = {World Scientific Publishing},
  title        = {{The excitation spectrum of the Bose gas in the Gross-Pitaevskii regime}},
  doi          = {10.1142/S0129055X20600065},
  volume       = {33},
  year         = {2021},
}

@article{7900,
  abstract     = {Hartree–Fock theory has been justified as a mean-field approximation for fermionic systems. However, it suffers from some defects in predicting physical properties, making necessary a theory of quantum correlations. Recently, bosonization of many-body correlations has been rigorously justified as an upper bound on the correlation energy at high density with weak interactions. We review the bosonic approximation, deriving an effective Hamiltonian. We then show that for systems with Coulomb interaction this effective theory predicts collective excitations (plasmons) in accordance with the random phase approximation of Bohm and Pines, and with experimental observation.},
  author       = {Benedikter, Niels P},
  issn         = {1793-6659},
  journal      = {Reviews in Mathematical Physics},
  number       = {1},
  publisher    = {World Scientific Publishing},
  title        = {{Bosonic collective excitations in Fermi gases}},
  doi          = {10.1142/s0129055x20600090},
  volume       = {33},
  year         = {2021},
}

@article{7901,
  abstract     = {We derive rigorously the leading order of the correlation energy of a Fermi gas in a scaling regime of high density and weak interaction. The result verifies the prediction of the random-phase approximation. Our proof refines the method of collective bosonization in three dimensions. We approximately diagonalize an effective Hamiltonian describing approximately bosonic collective excitations around the Hartree–Fock state, while showing that gapless and non-collective excitations have only a negligible effect on the ground state energy.},
  author       = {Benedikter, Niels P and Nam, Phan Thành and Porta, Marcello and Schlein, Benjamin and Seiringer, Robert},
  issn         = {1432-1297},
  journal      = {Inventiones Mathematicae},
  pages        = {885--979},
  publisher    = {Springer},
  title        = {{Correlation energy of a weakly interacting Fermi gas}},
  doi          = {10.1007/s00222-021-01041-5},
  volume       = {225},
  year         = {2021},
}

@article{7905,
  abstract     = {We investigate a sheaf-theoretic interpretation of stratification learning from geometric and topological perspectives. Our main result is the construction of stratification learning algorithms framed in terms of a sheaf on a partially ordered set with the Alexandroff topology. We prove that the resulting decomposition is the unique minimal stratification for which the strata are homogeneous and the given sheaf is constructible. In particular, when we choose to work with the local homology sheaf, our algorithm gives an alternative to the local homology transfer algorithm given in Bendich et al. (Proceedings of the 23rd Annual ACM-SIAM Symposium on Discrete Algorithms, pp. 1355–1370, ACM, New York, 2012), and the cohomology stratification algorithm given in Nanda (Found. Comput. Math. 20(2), 195–222, 2020). Additionally, we give examples of stratifications based on the geometric techniques of Breiding et al. (Rev. Mat. Complut. 31(3), 545–593, 2018), illustrating how the sheaf-theoretic approach can be used to study stratifications from both topological and geometric perspectives. This approach also points toward future applications of sheaf theory in the study of topological data analysis by illustrating the utility of the language of sheaf theory in generalizing existing algorithms.},
  author       = {Brown, Adam and Wang, Bei},
  issn         = {1432-0444},
  journal      = {Discrete and Computational Geometry},
  pages        = {1166--1198},
  publisher    = {Springer Nature},
  title        = {{Sheaf-theoretic stratification learning from geometric and topological perspectives}},
  doi          = {10.1007/s00454-020-00206-y},
  volume       = {65},
  year         = {2021},
}

@article{7925,
  abstract     = {In this paper, we introduce a relaxed CQ method with alternated inertial step for solving split feasibility problems. We give convergence of the sequence generated by our method under some suitable assumptions. Some numerical implementations from sparse signal and image deblurring are reported to show the efficiency of our method.},
  author       = {Shehu, Yekini and Gibali, Aviv},
  issn         = {1862-4480},
  journal      = {Optimization Letters},
  pages        = {2109--2126},
  publisher    = {Springer Nature},
  title        = {{New inertial relaxed method for solving split feasibilities}},
  doi          = {10.1007/s11590-020-01603-1},
  volume       = {15},
  year         = {2021},
}

@article{7939,
  abstract     = {We design fast deterministic algorithms for distance computation in the Congested Clique model. Our key contributions include:
    A (2+ϵ)-approximation for all-pairs shortest paths in O(log2n/ϵ) rounds on unweighted undirected graphs. With a small additional additive factor, this also applies for weighted graphs. This is the first sub-polynomial constant-factor approximation for APSP in this model.
    A (1+ϵ)-approximation for multi-source shortest paths from O(n−−√) sources in O(log2n/ϵ) rounds on weighted undirected graphs. This is the first sub-polynomial algorithm obtaining this approximation for a set of sources of polynomial size.

Our main techniques are new distance tools that are obtained via improved algorithms for sparse matrix multiplication, which we leverage to construct efficient hopsets and shortest paths. Furthermore, our techniques extend to additional distance problems for which we improve upon the state-of-the-art, including diameter approximation, and an exact single-source shortest paths algorithm for weighted undirected graphs in O~(n1/6) rounds. },
  author       = {Censor-Hillel, Keren and Dory, Michal and Korhonen, Janne and Leitersdorf, Dean},
  issn         = {1432-0452},
  journal      = {Distributed Computing},
  pages        = {463--487},
  publisher    = {Springer Nature},
  title        = {{Fast approximate shortest paths in the congested clique}},
  doi          = {10.1007/s00446-020-00380-5},
  volume       = {34},
  year         = {2021},
}

@inbook{7941,
  abstract     = {Expansion microscopy is a recently developed super-resolution imaging technique, which provides an alternative to optics-based methods such as deterministic approaches (e.g. STED) or stochastic approaches (e.g. PALM/STORM). The idea behind expansion microscopy is to embed the biological sample in a swellable gel, and then to expand it isotropically, thereby increasing the distance between the fluorophores. This approach breaks the diffraction barrier by simply separating the emission point-spread-functions of the fluorophores. The resolution attainable in expansion microscopy is thus directly dependent on the separation that can be achieved, i.e. on the expansion factor. The original implementation of the technique achieved an expansion factor of fourfold, for a resolution of 70–80 nm. The subsequently developed X10 method achieves an expansion factor of 10-fold, for a resolution of 25–30 nm. This technique can be implemented with minimal technical requirements on any standard fluorescence microscope, and is more easily applied for multi-color imaging than either deterministic or stochastic super-resolution approaches. This renders X10 expansion microscopy a highly promising tool for new biological discoveries, as discussed here, and as demonstrated by several recent applications.},
  author       = {Truckenbrodt, Sven M and Rizzoli, Silvio O.},
  booktitle    = {Methods in Cell Biology},
  isbn         = {978012820807-6},
  issn         = {0091-679X},
  pages        = {33--56},
  publisher    = {Elsevier},
  title        = {{Simple multi-color super-resolution by X10 microscopy}},
  doi          = {10.1016/bs.mcb.2020.04.016},
  volume       = {161},
  year         = {2021},
}

@article{8196,
  abstract     = {This paper aims to obtain a strong convergence result for a Douglas–Rachford splitting method with inertial extrapolation step for finding a zero of the sum of two set-valued maximal monotone operators without any further assumption of uniform monotonicity on any of the involved maximal monotone operators. Furthermore, our proposed method is easy to implement and the inertial factor in our proposed method is a natural choice. Our method of proof is of independent interest. Finally, some numerical implementations are given to confirm the theoretical analysis.},
  author       = {Shehu, Yekini and Dong, Qiao-Li and Liu, Lu-Lu and Yao, Jen-Chih},
  issn         = {1573-2924},
  journal      = {Optimization and Engineering},
  pages        = {2627--2653},
  publisher    = {Springer Nature},
  title        = {{New strong convergence method for the sum of two maximal monotone operators}},
  doi          = {10.1007/s11081-020-09544-5},
  volume       = {22},
  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{8248,
  abstract     = {We consider the following setting: suppose that we are given a manifold M in Rd with positive reach. Moreover assume that we have an embedded simplical complex A without boundary, whose vertex set lies on the manifold, is sufficiently dense and such that all simplices in A have sufficient quality. We prove that if, locally, interiors of the projection of the simplices onto the tangent space do not intersect, then A is a triangulation of the manifold, that is, they are homeomorphic.},
  author       = {Boissonnat, Jean-Daniel and Dyer, Ramsay and Ghosh, Arijit and Lieutier, Andre and Wintraecken, Mathijs},
  issn         = {1432-0444},
  journal      = {Discrete and Computational Geometry},
  pages        = {666--686},
  publisher    = {Springer Nature},
  title        = {{Local conditions for triangulating submanifolds of Euclidean space}},
  doi          = {10.1007/s00454-020-00233-9},
  volume       = {66},
  year         = {2021},
}

@article{8253,
  abstract     = {Brains process information in spiking neural networks. Their intricate connections shape the diverse functions these networks perform. In comparison, the functional capabilities of models of spiking networks are still rudimentary. This shortcoming is mainly due to the lack of insight and practical algorithms to construct the necessary connectivity. Any such algorithm typically attempts to build networks by iteratively reducing the error compared to a desired output. But assigning credit to hidden units in multi-layered spiking networks has remained challenging due to the non-differentiable nonlinearity of spikes. To avoid this issue, one can employ surrogate gradients to discover the required connectivity in spiking network models. However, the choice of a surrogate is not unique, raising the question of how its implementation influences the effectiveness of the method. Here, we use numerical simulations to systematically study how essential design parameters of surrogate gradients impact learning performance on a range of classification problems. We show that surrogate gradient learning is robust to different shapes of underlying surrogate derivatives, but the choice of the derivative’s scale can substantially affect learning performance. When we combine surrogate gradients with a suitable activity regularization technique, robust information processing can be achieved in spiking networks even at the sparse activity limit. Our study provides a systematic account of the remarkable robustness of surrogate gradient learning and serves as a practical guide to model functional spiking neural networks.},
  author       = {Zenke, Friedemann and Vogels, Tim P},
  issn         = {1530-888X},
  journal      = {Neural Computation},
  number       = {4},
  pages        = {899--925},
  publisher    = {MIT Press},
  title        = {{The remarkable robustness of surrogate gradient learning for instilling complex function in spiking neural networks}},
  doi          = {10.1162/neco_a_01367},
  volume       = {33},
  year         = {2021},
}

@article{8317,
  abstract     = {When can a polyomino piece of paper be folded into a unit cube? Prior work studied tree-like polyominoes, but polyominoes with holes remain an intriguing open problem. We present sufficient conditions for a polyomino with one or several holes to fold into a cube, and conditions under which cube folding is impossible. In particular, we show that all but five special “basic” holes guarantee foldability.},
  author       = {Aichholzer, Oswin and Akitaya, Hugo A. and Cheung, Kenneth C. and Demaine, Erik D. and Demaine, Martin L. and Fekete, Sándor P. and Kleist, Linda and Kostitsyna, Irina and Löffler, Maarten and Masárová, Zuzana and Mundilova, Klara and Schmidt, Christiane},
  issn         = {1879-081X},
  journal      = {Computational Geometry: Theory and Applications},
  publisher    = {Elsevier},
  title        = {{Folding polyominoes with holes into a cube}},
  doi          = {10.1016/j.comgeo.2020.101700},
  volume       = {93},
  year         = {2021},
}

@article{8338,
  abstract     = {Canonical parametrisations of classical confocal coordinate systems are introduced and exploited to construct non-planar analogues of incircular (IC) nets on individual quadrics and systems of confocal quadrics. Intimate connections with classical deformations of quadrics that are isometric along asymptotic lines and circular cross-sections of quadrics are revealed. The existence of octahedral webs of surfaces of Blaschke type generated by asymptotic and characteristic lines that are diagonally related to lines of curvature is proved theoretically and established constructively. Appropriate samplings (grids) of these webs lead to three-dimensional extensions of non-planar IC nets. Three-dimensional octahedral grids composed of planes and spatially extending (checkerboard) IC-nets are shown to arise in connection with systems of confocal quadrics in Minkowski space. In this context, the Laguerre geometric notion of conical octahedral grids of planes is introduced. The latter generalise the octahedral grids derived from systems of confocal quadrics in Minkowski space. An explicit construction of conical octahedral grids is presented. The results are accompanied by various illustrations which are based on the explicit formulae provided by the theory.},
  author       = {Akopyan, Arseniy and Bobenko, Alexander I. and Schief, Wolfgang K. and Techter, Jan},
  issn         = {1432-0444},
  journal      = {Discrete and Computational Geometry},
  pages        = {938--976},
  publisher    = {Springer Nature},
  title        = {{On mutually diagonal nets on (confocal) quadrics and 3-dimensional webs}},
  doi          = {10.1007/s00454-020-00240-w},
  volume       = {66},
  year         = {2021},
}

@article{8373,
  abstract     = {It is well known that special Kubo-Ando operator means admit divergence center interpretations, moreover, they are also mean squared error estimators for certain metrics on positive definite operators. In this paper we give a divergence center interpretation for every symmetric Kubo-Ando mean. This characterization of the symmetric means naturally leads to a definition of weighted and multivariate versions of a large class of symmetric Kubo-Ando means. We study elementary properties of these weighted multivariate means, and note in particular that in the special case of the geometric mean we recover the weighted A#H-mean introduced by Kim, Lawson, and Lim.},
  author       = {Pitrik, József and Virosztek, Daniel},
  issn         = {0024-3795},
  journal      = {Linear Algebra and its Applications},
  keywords     = {Kubo-Ando mean, weighted multivariate mean, barycenter},
  pages        = {203--217},
  publisher    = {Elsevier},
  title        = {{A divergence center interpretation of general symmetric Kubo-Ando means, and related weighted multivariate operator means}},
  doi          = {10.1016/j.laa.2020.09.007},
  volume       = {609},
  year         = {2021},
}

@article{8429,
  abstract     = {We develop a Bayesian model (BayesRR-RC) that provides robust SNP-heritability estimation, an alternative to marker discovery, and accurate genomic prediction, taking 22 seconds per iteration to estimate 8.4 million SNP-effects and 78 SNP-heritability parameters in the UK Biobank. We find that only ≤10% of the genetic variation captured for height, body mass index, cardiovascular disease, and type 2 diabetes is attributable to proximal regulatory regions within 10kb upstream of genes, while 12-25% is attributed to coding regions, 32–44% to introns, and 22-28% to distal 10-500kb upstream regions. Up to 24% of all cis and coding regions of each chromosome are associated with each trait, with over 3,100 independent exonic and intronic regions and over 5,400 independent regulatory regions having ≥95% probability of contributing ≥0.001% to the genetic variance of these four traits. Our open-source software (GMRM) provides a scalable alternative to current approaches for biobank data.},
  author       = {Patxot, Marion and Trejo Banos, Daniel and Kousathanas, Athanasios and Orliac, Etienne J and Ojavee, Sven E and Moser, Gerhard and Sidorenko, Julia and Kutalik, Zoltan and Magi, Reedik and Visscher, Peter M and Ronnegard, Lars and Robinson, Matthew Richard},
  issn         = {2041-1723},
  journal      = {Nature Communications},
  number       = {1},
  publisher    = {Springer Nature},
  title        = {{Probabilistic inference of the genetic architecture underlying functional enrichment of complex traits}},
  doi          = {10.1038/s41467-021-27258-9},
  volume       = {12},
  year         = {2021},
}

@article{8430,
  abstract     = {While recent advancements in computation and modelling have improved the analysis of complex traits, our understanding of the genetic basis of the time at symptom onset remains limited. Here, we develop a Bayesian approach (BayesW) that provides probabilistic inference of the genetic architecture of age-at-onset phenotypes in a sampling scheme that facilitates biobank-scale time-to-event analyses. We show in extensive simulation work the benefits BayesW provides in terms of number of discoveries, model performance and genomic prediction. In the UK Biobank, we find many thousands of common genomic regions underlying the age-at-onset of high blood pressure (HBP), cardiac disease (CAD), and type-2 diabetes (T2D), and for the genetic basis of onset reflecting the underlying genetic liability to disease. Age-at-menopause and age-at-menarche are also highly polygenic, but with higher variance contributed by low frequency variants. Genomic prediction into the Estonian Biobank data shows that BayesW gives higher prediction accuracy than other approaches.},
  author       = {Ojavee, Sven E and Kousathanas, Athanasios and Trejo Banos, Daniel and Orliac, Etienne J and Patxot, Marion and Lall, Kristi and Magi, Reedik and Fischer, Krista and Kutalik, Zoltan and Robinson, Matthew Richard},
  issn         = {2041-1723},
  journal      = {Nature Communications},
  number       = {1},
  publisher    = {Nature Research},
  title        = {{Genomic architecture and prediction of censored time-to-event phenotypes with a Bayesian genome-wide analysis}},
  doi          = {10.1038/s41467-021-22538-w},
  volume       = {12},
  year         = {2021},
}

