• Precise segmentation of the left ventricle (LV) within cardiac MRI images is a prerequisite for the quantitative measurement of heart function. However, this task is challenging due to the limited availability of labeled data and motion artifacts from cardiac imaging. In this work, we present an iterative segmentation algorithm for LV delineation. By coupling deep learning with a novel dynamic-based labeling scheme, we present a new methodology where a policy model is learned to guide an agent to travel over the the image, tracing out a boundary of the ROI -- using the magnitude difference of the Poincar\'e map as a stopping criterion. Our method is evaluated on two datasets, namely the Sunnybrook Cardiac Dataset (SCD) and data from the STACOM 2011 LV segmentation challenge. Our method outperforms the previous research over many metrics. In order to demonstrate the transferability of our method we present encouraging results over the STACOM 2011 data, when using a model trained on the SCD dataset.
  • On the road of searching for Majorana zero modes (MZMs) in topological insulator-based Josephson junctions, a highly-sought signature is the protected full transparency of electron transport through the junctions due to the existence of the MZMs, associated with complete gap closing between the electron-like and hole-like Andreev bound states (ABSs). Here, we present direct experimental evidence of gap closing and full transparency in single Josephson junctions constructed on the surface of three-dimensional topological insulator (3D TI) Bi$_2$Te$_3$. Our results demonstrate that the 2D surface of 3D TIs provides a promising platform for hosting and manipulating MZMs.
  • Fine control of the dynamics of a quantum system is the key element to perform quantum information processing and coherent manipulations for atomic and molecular systems. In this paper we propose a control protocol using a tangent-pulse driven model and demonstrate that it indicates a desirable design, i.e., of being both fast and accurate for population transfer. As opposed to other existing strategies, a remarkable character of the present scheme is that high velocity of the nonadiabatic evolution itself not only will not lead to unwanted transitions but also can suppress the error caused by the truncation of the driving pulse.
  • For millimeter-wave (mm-wave) buffer-aided tandem networks consisting of relay nodes and multiple channels per hop, we consider two traffic allocation schemes, namely local allocation and global allocation, and investigate the end-to-end latency of a file transfer. We formulate the problem for generic multi-hop queuing systems and subsequently derive closed-form expressions of the end-to-end latency. We quantify the advantages of the global allocation scheme relative to its local allocation counterpart, and we conduct an asymptotic analysis on the performance gain when the number of channels in each hop increases to infinity. The traffic allocations and the analytical delay performance are validated through simulations. Furthermore, taking a specific two-hop mm-wave network as an example, we derive lower bounds on the average end-to-end latency, where Nakagami-$m$ fading is considered. Numerical results demonstrate that, compared to the local allocation scheme, the advantage of global allocation grows as the number of relay nodes increases, at the expense of higher complexity that linearly increases with the number of relay nodes. It is also demonstrated that a proper deployment of relay nodes in a linear mm-wave network plays an important role in reducing the average end-to-end latency, and the average latency decays as the mm-wave channels become more deterministic. These findings provide insights for designing multi-hop mm-wave networks with low end-to-end latency.
  • Deep convolutional neural networks (CNNs) have been intensively used for multi-class segmentation of data from different modalities and achieved state-of-the-art performances. However, a common problem when dealing with large, high resolution 3D data is that the volumes input into the deep CNNs has to be either cropped or downsampled due to limited memory capacity of computing devices. These operations lead to loss of resolution and increment of class imbalance in the input data batches, which can downgrade the performances of segmentation algorithms. Inspired by the architecture of image super-resolution CNN (SRCNN) and self-normalization network (SNN), we developed a two-stage modified Unet framework that simultaneously learns to detect a ROI within the full volume and to classify voxels without losing the original resolution. Experiments on a variety of multi-modal volumes demonstrated that, when trained with a simply weighted dice coefficients and our customized learning procedure, this framework shows better segmentation performances than state-of-the-art Deep CNNs with advanced similarity metrics.
  • We revisit skip-gram negative sampling (SGNS), a popular neural-network based approach to learning distributed word representation. We first point out the ambiguity issue undermining the SGNS model, in the sense that the word vectors can be entirely distorted without changing the objective value. To resolve this issue, we rectify the SGNS model with quadratic regularization. A theoretical justification, which provides a novel insight into quadratic regularization, is presented. Preliminary experiments are also conducted on Google's analytical reasoning task to support the modified SGNS model.
  • Drones are driving numerous and evolving use cases, and creating transformative socio-economic benefits. Drone operation needs wireless connectivity for communication between drones and ground control systems, among drones, and between drones and air traffic management systems. Mobile networks are well positioned to identify, track, and control the growing fleet of drones. The wide-area, quality, and secure connectivity provided by mobile networks can enhance the efficiency and effectiveness of drone operations beyond visual line-of-sight range. In this article, we elaborate how the drone ecosystem can benefit from mobile technologies, summarize key capabilities required by drone applications, and analyze the service requirements on mobile networks. We present field trial results collected in LTE-Advanced networks to gain insights into the capabilities of the current 4G+ networks for connected drones and share our vision on how 5G networks can further support diversified drone applications.
  • Most recently, a brand new phosphorus allotrope called green phosphorus has been predicted, which has a direct band gap up to 2.4 eV, and its single-layer form termed green phosphorene shows high stability. Here the mechanical properties and the uniaxial strain effect on the electronic band structure of green phosphorene along two perpendicular in-plane directions are investigated. Remarkably, we find that this material can sustain a tensile strain in the armchair direction up to a threshold of 35\% which is larger than that of black phosphorene, suggesting that green phosphorene is more puckered. Our calculations also show the Young's modulus and Poisson's ratio in the zigzag direction are four times larger than those in the armchair direction, which confirm the anisotropy of the material. Furthermore, the uniaxial strain can trigger the direct-indirect band gap transition for green phosphorene and the critical strains for the band gap transition are revealed.
  • This paper investigates two strategies to reduce the communication delay in future wireless networks: traffic dispersion and network densification. A hybrid scheme that combines these two strategies is also considered. The probabilistic delay and effective capacity are used to evaluate performance. For probabilistic delay, the violation probability of delay, i.e., the probability that the delay exceeds a given tolerance level, is characterized in terms of upper bounds, which are derived by applying stochastic network calculus theory. In addition, to characterize the maximum affordable arrival traffic for mmWave systems, the effective capacity, i.e., the service capability with a given quality-of-service (QoS) requirement, is studied. The derived bounds on the probabilistic delay and effective capacity are validated through simulations. These numerical results show that, for a given average system gain, traffic dispersion, network densification, and the hybrid scheme exhibit different potentials to reduce the end-to-end communication delay. For instance, traffic dispersion outperforms network densification, given high average system gain and arrival rate, while it could be the worst option, otherwise. Furthermore, it is revealed that, increasing the number of independent paths and/or relay density is always beneficial, while the performance gain is related to the arrival rate and average system gain, jointly. Therefore, a proper transmission scheme should be selected to optimize the delay performance, according to the given conditions on arrival traffic and system service capability.
  • Heterogeneous network (HetNet) is a key enabler to largely boost network coverage and capacity in the forthcoming fifth-generation (5G) and beyond. To support the explosively growing mobile data volumes, wireless communications with millimeter-wave (mm-wave) radios have attracted massive attention, which is widely considered as a promising candidate in 5G HetNets. In this article, we give an overview on the end-to-end latency of HetNets with mm-wave communications. In general, it is rather challenging for formulating and optimizing the delay problem with buffers in mm-wave communications, since conventional graph-based network optimization techniques are not applicable when queues are considered. Toward this end, we develop an adaptive low-latency strategy, which uses cooperative networking to reduce the end-to-end latency. Then, we evaluate the performance of the introduced strategy. Results reveal the importance of proper cooperative networking in reducing the end-to-end latency. In addition, we have identified several challenges in future research for low-latency mm-wave HetNets.
  • The weight decision problem, which requires to determine the Hamming weight of a given binary string, is a natural and important problem, with applications in cryptanalysis, coding theory, fault-tolerant circuit design and so on. In particular, both Deutsch-Jozsa problem and Grover search problem can be interpreted as special cases of weight decision problems. In this work, we investigate the exact quantum query complexity of weight decision problems, where the quantum algorithm must always output the correct answer. More specifically we consider a partial Boolean function which distinguishes whether the Hamming weight of the length-$n$ input is $k$ or it is $l$. Our contribution includes both upper bounds and lower bounds for the precise number of queries. Furthermore, for most choices of $(\frac{k}{n},\frac{l}{n})$ and sufficiently large $n$, the gap between our upper and lower bounds is no more than one. To get the results, we first build the connection between Chebyshev polynomials and our problem, then determine all the boundary cases of $(\frac{k}{n},\frac{l}{n})$ with matching upper and lower bounds, and finally we generalize to other cases via a new \emph{quantum padding} technique. This quantum padding technique can be of independent interest in designing other quantum algorithms.
  • Using the determinant quantum Monte-Carlo method, we elucidate the strain tuning of edge magnetism in zigzag graphene nanoribbons. Our intensive numerical results show that a relatively weak Coulomb interaction may induce a ferromagnetic-like behaviour with a proper strain, and the edge magnetism can be enhanced greatly as the strain along the zigzag edge increases, which provides another way to control graphene magnetism even at room temperature.
  • We experimentally observe soliton self-organization and pulsation in a passively mode-locked fiber laser. The optomechanical interaction in the optical fiber is key to the formation of equidistant soliton bunches. These solitons simultaneously undergo a pulsation process with a period corresponding to tens of the cavity round trip time. Using the dispersive Fourier transformation technique, we find that the Kelly sidebands in the shot-to-shot spectra appear periodically, synchronizing with the pulsation.
  • We present a detailed spectral analysis of the brightest Active Galactic Nuclei (AGN) identified in the 7Ms Chandra Deep Field South (CDF-S) survey over a time span of 16 years. Using a model of an intrinsically absorbed power-law plus reflection, with possible soft excess and narrow Fe K$\alpha$ line, we perform a systematic X-ray spectral analysis, both on the total 7Ms exposure and in four different periods with lengths of 2-21 months. With this approach, we not only present the power-law slopes, column densities $N_H$, observed fluxes, and absorption-corrected 2-10~keV luminosities $L_X$ for our sample of AGNs, but also identify significant spectral variabilities among them on time scales of years. We find that the $N_H$ variabilities can be ascribed to two different types of mechanisms, either flux-driven or flux-independent. We also find that the correlation between the narrow Fe line EW and $N_H$ can be well explained by the continuum suppression with increasing $N_H$. Accounting for the sample incompleteness and bias, we measure the intrinsic distribution of $N_H$ for the CDF-S AGN population and present re-selected subsamples which are complete with respect to $N_H$. The $N_H$-complete subsamples enable us to decouple the dependences of $N_H$ on $L_X$ and on redshift. Combining our data with that from C-COSMOS, we confirm the anti-correlation between the average $N_H$ and $L_X$ of AGN, and find a significant increase of the AGN obscured fraction with redshift at any luminosity. The obscured fraction can be described as $f_{obscured}\thickapprox 0.42\ (1+z)^{0.60}$.
  • A major challenge in brain tumor treatment planning and quantitative evaluation is determination of the tumor extent. The noninvasive magnetic resonance imaging (MRI) technique has emerged as a front-line diagnostic tool for brain tumors without ionizing radiation. Manual segmentation of brain tumor extent from 3D MRI volumes is a very time-consuming task and the performance is highly relied on operator's experience. In this context, a reliable fully automatic segmentation method for the brain tumor segmentation is necessary for an efficient measurement of the tumor extent. In this study, we propose a fully automatic method for brain tumor segmentation, which is developed using U-Net based deep convolutional networks. Our method was evaluated on Multimodal Brain Tumor Image Segmentation (BRATS 2015) datasets, which contain 220 high-grade brain tumor and 54 low-grade tumor cases. Cross-validation has shown that our method can obtain promising segmentation efficiently.
  • Purpose: Atrial fibrillation (AF) is the most common cardiac arrhythmia and is correlated with increased morbidity and mortality. It is associated with atrial fibrosis, which may be assessed non-invasively using late gadolinium-enhanced (LGE) magnetic resonance imaging (MRI) where scar tissue is visualised as a region of signal enhancement. In this study, we proposed a novel fully automatic pipeline to achieve an accurate and objective atrial scarring segmentation and assessment of LGE MRI scans for the AF patients. Methods: Our fully automatic pipeline uniquely combined: (1) a multi-atlas based whole heart segmentation (MA-WHS) to determine the cardiac anatomy from an MRI Roadmap acquisition which is then mapped to LGE MRI, and (2) a super-pixel and supervised learning based approach to delineate the distribution and extent of atrial scarring in LGE MRI. Results: Both our MA-WHS and atrial scarring segmentation showed accurate delineations of cardiac anatomy (mean Dice = 89%) and atrial scarring (mean Dice =79%) respectively compared to the established ground truth from manual segmentation. Compared with previously studied methods with manual interventions, our innovative pipeline demonstrated comparable results, but was computed fully automatically. Conclusion: The proposed segmentation methods allow LGE MRI to be used as an objective assessment tool for localisation, visualisation and quantification of atrial scarring.
  • Fast Magnetic Resonance Imaging (MRI) is highly in demand for many clinical applications in order to reduce the scanning cost and improve the patient experience. This can also potentially increase the image quality by reducing the motion artefacts and contrast washout. However, once an image field of view and the desired resolution are chosen, the minimum scanning time is normally determined by the requirement of acquiring sufficient raw data to meet the Nyquist-Shannon sampling criteria. Compressive Sensing (CS) theory has been perfectly matched to the MRI scanning sequence design with much less required raw data for the image reconstruction. Inspired by recent advances in deep learning for solving various inverse problems, we propose a conditional Generative Adversarial Networks-based deep learning framework for de-aliasing and reconstructing MRI images from highly undersampled data with great promise to accelerate the data acquisition process. By coupling an innovative content loss with the adversarial loss our de-aliasing results are more realistic. Furthermore, we propose a refinement learning procedure for training the generator network, which can stabilise the training with fast convergence and less parameter tuning. We demonstrate that the proposed framework outperforms state-of-the-art CS-MRI methods, in terms of reconstruction error and perceptual image quality. In addition, our method can reconstruct each image in 0.22ms--0.37ms, which is promising for real-time applications.
  • High speed wireless access on 60 GHz spectrum relies on high-gain directional antennas to overcome the severe signal attenuation. However, perfect alignment between transmitting and receiving antenna beams is rare in practice and overheard signals from concurrent transmissions may cause significant interference. In this paper we analyze the impact of antenna beam misalignment on the system performance of 60 GHz wireless access. We quantify the signal power loss caused by beam misalignment and the interference power accumulated from neighboring concurrent transmissions whose signals are leaked either via the main-beam pointing in the similar direction or via side-lobe emission, and derive the probability distribution of the signal to interference plus noise power ratio (SINR). For scenarios where interfering transmitters are distributed uniformly at random, we derive upper and lower bounds on the cumulative distribution function (abbreviated as CDF or c.d.f.) of SINR, which can be easily applied to evaluate system performance. We validate our analytical results by simulations where random nodes are uniformly distributed within a circular hall, and evaluate the sensitivity of average throughput and outage probability against two parameters: the half-power (3 dB) beamwidth to main-lobe beamwidth ratio and the beam misalignment deviation to main-lobe beamwidth ratio. Our results indicate that the derived lower bound performs well when the half-power beamwidth to main-lobe beamwidth ratio or the number of concurrent transmission links is small. When the number of active links is high, it is desirable in antenna design to balance the degradation caused by beam misalignment (wider beam is better) and the interference from concurrent transmission (narrower beam is better).
  • There have been continuous efforts in searching for unconventional superconductivity over the past five decades. Compared to the well-established d-wave superconductivity in cuprates, the existence of superconductivity with other high-angular-momentum pairing symmetries is less conclusive. Bi/Ni epitaxial bilayer is a potential unconventional superconductor with broken time reversal symmetry (TRS), for that it demonstrates superconductivity and ferromagnetism simultaneously at low temperatures. We employ a specially designed superconducting quantum interference device (SQUID) to detect, on the Bi/Ni bilayer, the orbital magnetic moment which is expected if the TRS is broken. An anomalous hysteretic magnetic response has been observed in the superconducting state, providing the evidence for the existence of chiral superconducting domains in the material.
  • Artificial neural networks have gone through a recent rise in popularity, achieving state-of-the-art results in various fields, including image classification, speech recognition, and automated control. Both the performance and computational complexity of such models are heavily dependant on the design of characteristic hyper-parameters (e.g., number of hidden layers, nodes per layer, or choice of activation functions), which have traditionally been optimized manually. With machine learning penetrating low-power mobile and embedded areas, the need to optimize not only for performance (accuracy), but also for implementation complexity, becomes paramount. In this work, we present a multi-objective design space exploration method that reduces the number of solution networks trained and evaluated through response surface modelling. Given spaces which can easily exceed 1020 solutions, manually designing a near-optimal architecture is unlikely as opportunities to reduce network complexity, while maintaining performance, may be overlooked. This problem is exacerbated by the fact that hyper-parameters which perform well on specific datasets may yield sub-par results on others, and must therefore be designed on a per-application basis. In our work, machine learning is leveraged by training an artificial neural network to predict the performance of future candidate networks. The method is evaluated on the MNIST and CIFAR-10 image datasets, optimizing for both recognition accuracy and computational complexity. Experimental results demonstrate that the proposed method can closely approximate the Pareto-optimal front, while only exploring a small fraction of the design space.
  • Understanding phonon scattering by topological defects in graphene is of particular interest for thermal management in graphene-based devices. We present a study that quantifies the roles of the different mechanisms governing defect phonon scattering by comparing the effects of ten different defect structures using molecular dynamics. Our results show that phonon scattering is mainly influenced by mass density difference, with general trends governed by the defect formation energy and typical softening behaviors in the phonon density of state. The phonon scattering cross-section is found to be far larger than that geometrically occupied by the defects. We also show that the lattice thermal conductivity can be reduced by a factor of up to ~30 in the presence of the grain boundaries formed by these defects.
  • The Astropy Problem (1610.03159)

    Demitri Muna, Michael Alexander, Alice Allen, Richard Ashley, Daniel Asmus, Ruyman Azzollini, Michele Bannister, Rachael Beaton, Andrew Benson, G. Bruce Berriman, Maciej Bilicki, Peter Boyce, Joanna Bridge, Jan Cami, Eryn Cangi, Xian Chen, Nicholas Christiny, Christopher Clark, Michelle Collins, Johan Comparat, Neil Cook, Darren Croton, Isak Delberth Davids, Éric Depagne, John Donor, Leonardo A. dos Santos, Stephanie Douglas, Alan Du, Meredith Durbin, Dawn Erb, Daniel Faes, J. G. Fernández-Trincado, Anthony Foley, Sotiria Fotopoulou, Søren Frimann, Peter Frinchaboy, Rafael Garcia-Dias, Artur Gawryszczak, Elizabeth George, Sebastian Gonzalez, Karl Gordon, Nicholas Gorgone, Catherine Gosmeyer, Katie Grasha, Perry Greenfield, Rebekka Grellmann, James Guillochon, Mark Gurwell, Marcel Haas, Alex Hagen, Daryl Haggard, Tim Haines, Patrick Hall, Wojciech Hellwing, Edmund Christian Herenz, Samuel Hinton, Renee Hlozek, John Hoffman, Derek Holman, Benne Willem Holwerda, Anthony Horton, Cameron Hummels, Daniel Jacobs, Jens Juel Jensen, David Jones, Arna Karick, Luke Kelley, Matthew Kenworthy, Ben Kitchener, Dominik Klaes, Saul Kohn, Piotr Konorski, Coleman Krawczyk, Kyler Kuehn, Teet Kuutma, Michael T. Lam, Richard Lane, Jochen Liske, Diego Lopez-Camara, Katherine Mack, Sam Mangham, Qingqing Mao, David J. E. Marsh, Cecilia Mateu, Loïc Maurin, James McCormac, Ivelina Momcheva, Hektor Monteiro, Michael Mueller, Roberto Munoz, Rohan Naidu, Nicholas Nelson, Christian Nitschelm, Chris North, Juan Nunez-Iglesias, Sara Ogaz, Russell Owen, John Parejko, Vera Patrício, Joshua Pepper, Marshall Perrin, Timothy Pickering, Jennifer Piscionere, Richard Pogge, Radek Poleski, Alkistis Pourtsidou, Adrian M. Price-Whelan, Meredith L. Rawls, Shaun Read, Glen Rees, Hanno Rein, Thomas Rice, Signe Riemer-Sørensen, Naum Rusomarov, Sebastian F. Sanchez, Miguel Santander-García, Gal Sarid, William Schoenell, Aleks Scholz, Robert L. Schuhmann, William Schuster, Peter Scicluna, Marja Seidel, Lijing Shao, Pranav Sharma, Aleksandar Shulevski, David Shupe, Cristóbal Sifón, Brooke Simmons, Manodeep Sinha, Ian Skillen, Bjoern Soergel, Thomas Spriggs, Sundar Srinivasan, Abigail Stevens, Ole Streicher, Eric Suchyta, Joshua Tan, O. Grace Telford, Romain Thomas, Chiara Tonini, Grant Tremblay, Sarah Tuttle, Tanya Urrutia, Sam Vaughan, Miguel Verdugo, Alexander Wagner, Josh Walawender, Andrew Wetzel, Kyle Willett, Peter K. G. Williams, Guang Yang, Guangtun Zhu, Andrea Zonca
    Oct. 11, 2016 physics.soc-ph, astro-ph.IM
    The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.
  • We exploit the 7 Ms \textit{Chandra} observations in the \chandra\,Deep Field-South (\mbox{CDF-S}), the deepest X-ray survey to date, coupled with CANDELS/GOODS-S data, to measure the total X-ray emission arising from 2076 galaxies at $3.5\leq z < 6.5$. This aim is achieved by stacking the \textit{Chandra} data at the positions of optically selected galaxies, reaching effective exposure times of $\geq10^9\mathrm{s}$. We detect significant ($>3.7\sigma$) X-ray emission from massive galaxies at $z\approx4$. We also report the detection of massive galaxies at $z\approx5$ at a $99.7\%$ confidence level ($2.7\sigma$), the highest significance ever obtained for X-ray emission from galaxies at such high redshifts. No significant signal is detected from galaxies at even higher redshifts. The stacking results place constraints on the BHAD associated with the known high-redshift galaxy samples, as well as on the SFRD at high redshift, assuming a range of prescriptions for X-ray emission due to X- ray binaries. We find that the X-ray emission from our sample is likely dominated by processes related to star formation. Our results show that low-rate mass accretion onto SMBHs in individually X-ray-undetected galaxies is negligible, compared with the BHAD measured for samples of X-ray detected AGN, for cosmic SMBH mass assembly at high redshift. We also place, for the first time, constraints on the faint-end of the AGN X-ray luminosity function ($\mathrm{logL_X\sim42}$) at $z>4$, with evidence for fairly flat slopes. The implications of all of these findings are discussed in the context of the evolution of the AGN population at high redshift.
  • To meet the ever-increasing demands on higher throughput and better network delay performance, 60 GHZ networking is proposed as a promising solution for the next generation of wireless communications. To successfully deploy such networks, its important to understand their performance first. However, due to the unique fading characteristic of the 60 GHz channel, the characterization of the corresponding service process, offered by the channel, using the conventional methodologies may not be tractable. In this work, we provide an alternative approach to derive a closed-form expression that characterizes the cumulative service process of the 60 GHz channel in terms of the moment generating function (MGF) of its instantaneous channel capacity. We then use this expression to derive probabilistic upper bounds on the backlog and delay that are experienced by a flow traversing this network, using results from the MGF-based network calculus. The computed bounds are validated using simulation. We provide numerical results for different networking scenarios and for different traffic and channel parameters and we show that the 60 GHz wireless network is capable of satisfying stringent quality-of-Service (QoS) requirements, in terms of network delay and reliability. With this analysis approach at hand, a larger scale 60 GHz network design and optimization is possible.
  • Searching for room temperature ferromagnetic semiconductors has evolved into a broad field of material science and spintronics for decades, nevertheless, these novel states remain rare. Phosphorene, a monolayer black phosphorus with a puckered honeycomb lattice structure possessing a finite band gap and high carrier mobility, has been synthesized recently. Here we show, by means of two different large scale quantum Monte-Carlo methods, that relatively weak interactions can lead to remarkable edge magnetism in the phosphorene nanoribbons. The ground state constrained path quantum Monte-Carlo simulations reveal strong ferromagnetic correlations along the zigzag edges, and the finite temperature determinant quantum Monte-Carlo calculations show a high Curie temperature up to room temperature.