• ### On the dimension effect of regularized linear discriminant analysis(1710.03136)

Sept. 21, 2018 math.ST, stat.TH
This paper studies the dimension effect of the linear discriminant analysis (LDA) and the regularized linear discriminant analysis (RLDA) classifiers for large dimensional data where the observation dimension $p$ is of the same order as the sample size $n$. More specifically, built on properties of the Wishart distribution and recent results in random matrix theory, we derive explicit expressions for the asymptotic misclassification errors of LDA and RLDA respectively, from which we gain insights of how dimension affects the performance of classification and in what sense. Motivated by these results, we propose adjusted classifiers by correcting the bias brought by the unequal sample sizes. The bias-corrected LDA and RLDA classifiers are shown to have smaller misclassification rates than LDA and RLDA respectively. Several interesting examples are discussed in detail and the theoretical results on dimension effect are illustrated via extensive simulation studies.
• ### Numerical Complete Solution for Random Genetic Drift by Energetic Variational Approach(1803.09436)

March 26, 2018 math.NA
In this paper, we focus on numerical solutions for random genetic drift problem, which is governed by a degenerated convection-dominated parabolic equation. Due to the fixation phenomenon of genes, Dirac delta singularities will develop at boundary points as time evolves. Based on an energetic variational approach (EnVarA), a balance between the maximal dissipation principle (MDP) and least action principle (LAP), we obtain the trajectory equation. In turn, a numerical scheme is proposed using a convex splitting technique, with the unique solvability (on a convex set) and the energy decay property (in time) justified at a theoretical level. Numerical examples are presented for cases of pure drift and drift with semi-selection. The remarkable advantage of this method is its ability to catch the Dirac delta singularity close to machine precision over any equidistant grid.
• ### Dynamic Output Feedback Guaranteed-Cost Synchronization for Multiagent Networks with Given Cost Budgets(1802.07923)

Feb. 22, 2018 cs.DC, cs.SY
The current paper addresses the distributed guaranteed-cost synchronization problems for general high-order linear multiagent networks. Existing works on the guaranteed-cost synchronization usually require all state information of neighboring agents and cannot give the cost budget previously. For both leaderless and leader-following interaction topologies, the current paper firstly proposes a dynamic output feedback synchronization protocol with guaranteed-cost constraints, which can realize the tradeoff design between the energy consumption and the synchronization regulation performance with the given cost budget. Then, according to different structure features of interaction topologies, leaderless and leader-following guaranteed-cost synchronization analysis and design criteria are presented, respectively, and an algorithm is proposed to deal with the impacts of nonlinear terms by using both synchronization analysis and design criteria. Especially, an explicit expression of the synchronization function is shown for leaderless cases, which is independent of protocol states and the given cost budget. Finally, numerical examples are presented to demonstrate theoretical results.
• ### Distinct differences in the nanoscale behaviors of the twist-bend liquid crystal phase of a flexible linear trimer and homologous dimer(1710.00922)

Feb. 19, 2018 cond-mat.soft
We synthesized the liquid crystal dimer and trimer members of a series of flexible linear oligomers and characterized their microscopic and nanoscopic properties using resonant soft x-ray scattering and a number of other experimental techniques. On the microscopic scale, the twist-bend phases of the dimer and trimer appear essentially identical. However, while the liquid crystal dimer exhibits a temperature-dependent variation of its twist-bend helical pitch varying from 100 - 170 {\AA} on heating, the trimer exhibits an essentially temperature-independent pitch of 66 {\AA}, significantly shorter than those reported for other twist-bend forming materials in the literature. We attribute this to a specific combination of intrinsic conformational bend of the trimer molecules and a sterically favorable intercalation of the trimers over a commensurate fraction (two-thirds) of the molecular length. We develop a geometric model of the twist-bend phase for these materials with the molecules arranging into helical chain structures, and we fully determine their respective geometric parameters.
• ### Laser guided ionic wind(1801.06778)

We report on a method to experimentally generate ionic wind by coupling an external high voltage electric field with an intense femtosecond laser induced air plasma filament. The measured ionic wind velocity could be as strong as >4 m/s. It could be optimized by changing the applied electric field and the laser induced plasma channel. The experimental observation was qualitatively confirmed by a numerical simulation of spatial distribution of the electric field. This technique is robust and free from sharp metallic electrodes for coronas; it opens a way to optically generate ionic wind at a distance.
• ### Composite Behavioral Modeling for Identity Theft Detection in Online Social Networks(1801.06825)

Jan. 21, 2018 cs.CR, cs.SI
In this work, we aim at building a bridge from poor behavioral data to an effective, quick-response, and robust behavior model for online identity theft detection. We concentrate on this issue in online social networks (OSNs) where users usually have composite behavioral records, consisting of multi-dimensional low-quality data, e.g., offline check-ins and online user generated content (UGC). As an insightful result, we find that there is a complementary effect among different dimensions of records for modeling users' behavioral patterns. To deeply exploit such a complementary effect, we propose a joint model to capture both online and offline features of a user's composite behavior. We evaluate the proposed joint model by comparing with some typical models on two real-world datasets: Foursquare and Yelp. In the widely-used setting of theft simulation (simulating thefts via behavioral replacement), the experimental results show that our model outperforms the existing ones, with the AUC values $0.956$ in Foursquare and $0.947$ in Yelp, respectively. Particularly, the recall (True Positive Rate) can reach up to $65.3\%$ in Foursquare and $72.2\%$ in Yelp with the corresponding disturbance rate (False Positive Rate) below $1\%$. It is worth mentioning that these performances can be achieved by examining only one composite behavior (visiting a place and posting a tip online simultaneously) per authentication, which guarantees the low response latency of our method. This study would give the cybersecurity community new insights into whether and how a real-time online identity authentication can be improved via modeling users' composite behavioral patterns.
• ### Artificial Impostors for Location Privacy Preservation(1801.06827)

Jan. 21, 2018 cs.CR, cs.SI
The progress of location-based services has led to serious concerns on location privacy leakage. For effective and efficient location privacy preservation (LPP), existing methods are still not fully competent. They are often vulnerable under the identification attack with side information, or hard to be implemented due to the high computational complexity. In this paper, we pursue the high protection efficacy and low computational complexity simultaneously. We propose a \emph{scalable} LPP method based on the paradigm of counterfeiting locations. To make fake locations extremely plausible, we forge them through synthesizing \emph{artificial impostors} (AIs). The AIs refer to the synthesized traces which have similar semantic features to the actual traces, and do \emph{not} contain any target location. Two dedicated techniques are devised: the \emph{sampling-based synthesis method} and \emph{population-level semantic model}. They play significant roles in two critical steps of synthesizing AIs. We conduct experiments on real datasets in two cities (Shanghai, China and Asturias, Spain) to validate the high efficacy and scalability of the proposed method. In these two datasets, the experimental results show that our method achieves the preservation efficacy of $97.65\%$ and $96.12\%$, and its run time of building the generators is only $230.47$ and $215.92$ seconds, respectively. This study would give the research community new insights into improving the practicality of the state-of-the-art LPP paradigm via counterfeiting locations.
• ### Traffic Sign Timely Visual Recognizability Evaluation Based on 3D Measurable Point Clouds(1710.03553)

Oct. 10, 2017 cs.CV
The timely provision of traffic sign information to drivers is essential for the drivers to respond, to ensure safe driving, and to avoid traffic accidents in a timely manner. We proposed a timely visual recognizability quantitative evaluation method for traffic signs in large-scale transportation environments. To achieve this goal, we first address the concept of a visibility field to reflect the visible distribution of three-dimensional (3D) space and construct a traffic sign Visibility Evaluation Model (VEM) to measure the traffic sign visibility for a given viewpoint. Then, based on the VEM, we proposed the concept of the Visual Recognizability Field (VRF) to reflect the visual recognizability distribution in 3D space and established a Visual Recognizability Evaluation Model (VREM) to measure a traffic sign visual recognizability for a given viewpoint. Next, we proposed a Traffic Sign Timely Visual Recognizability Evaluation Model (TSTVREM) by combining VREM, the actual maximum continuous visual recognizable distance, and traffic big data to measure a traffic sign visual recognizability in different lanes. Finally, we presented an automatic algorithm to implement the TSTVREM model through traffic sign and road marking detection and classification, traffic sign environment point cloud segmentation, viewpoints calculation, and TSTVREM model realization. The performance of our method for traffic sign timely visual recognizability evaluation is tested on three road point clouds acquired by a mobile laser scanning system (RIEGL VMX-450) according to Road Traffic Signs and Markings (GB 5768-1999 in China), showing that our method is feasible and efficient.
• ### Exploring the dark matter inelastic frontier with 79.6 days of PandaX-II data(1708.05825)

Oct. 3, 2017 hep-ph, hep-ex, astro-ph.CO
We report here the results of searching for inelastic scattering of dark matter (initial and final state dark matter particles differ by a small mass splitting) with nucleon with the first 79.6-day of PandaX-II data (Run 9). We set the upper limits for the spin independent WIMP-nucleon scattering cross section up to a mass splitting of 300 keV/c$^2$ at two benchmark dark matter masses of 1 and 10 TeV/c$^2$.
• ### Dark Matter Results From 54-Ton-Day Exposure of PandaX-II Experiment(1708.06917)

We report a new search of weakly interacting massive particles (WIMPs) using the combined low background data sets in 2016 and 2017 from the PandaX-II experiment in China. The latest data set contains a new exposure of 77.1 live day, with the background reduced to a level of 0.8$\times10^{-3}$ evt/kg/day, improved by a factor of 2.5 in comparison to the previous run in 2016. No excess events were found above the expected background. With a total exposure of 5.4$\times10^4$ kg day, the most stringent upper limit on spin-independent WIMP-nucleon cross section was set for a WIMP with mass larger than 100 GeV/c$^2$, with the lowest exclusion at 8.6$\times10^{-47}$ cm$^2$ at 40 GeV/c$^2$.
• ### RRA: Recurrent Residual Attention for Sequence Learning(1709.03714)

Sept. 12, 2017 cs.AI, cs.LG
In this paper, we propose a recurrent neural network (RNN) with residual attention (RRA) to learn long-range dependencies from sequential data. We propose to add residual connections across timesteps to RNN, which explicitly enhances the interaction between current state and hidden states that are several timesteps apart. This also allows training errors to be directly back-propagated through residual connections and effectively alleviates gradient vanishing problem. We further reformulate an attention mechanism over residual connections. An attention gate is defined to summarize the individual contribution from multiple previous hidden states in computing the current state. We evaluate RRA on three tasks: the adding problem, pixel-by-pixel MNIST classification and sentiment analysis on the IMDB dataset. Our experiments demonstrate that RRA yields better performance, faster convergence and more stable training compared to a standard LSTM network. Furthermore, RRA shows highly competitive performance to the state-of-the-art methods.
• ### Metasurface-assisted phase-matching-free second harmonic generation in lithium niobate waveguides(1708.02889)

Aug. 9, 2017 physics.optics
The phase-matching condition is a key aspect in nonlinear wavelength conversion processes, which requires the momenta of the photons involved in the processes to be conserved. Conventionally, nonlinear phase matching is achieved using either birefringent or periodically poled nonlinear crystals, which requires careful dispersion engineering and are usually narrowband. In recent years, metasurfaces consisting of densely packed arrays of optical antennas have been demonstrated to provide an effective optical momentum to bend light in arbitrary ways. Here, we demonstrate that gradient metasurface structures consisting of phased array antennas are able to circumvent the phase-matching requirement in on-chip nonlinear wavelength conversion. We experimentally demonstrate phase-matching-free second harmonic generation over many coherent lengths in thin film lithium niobate waveguides patterned with the gradient metasurfaces. Efficient second-harmonic generation (1660% W-1cm-2) in the metasurface-based devices was observed over a wide range of pump wavelengths (1580-1650 nm).
• ### A Second-Order Energy Stable Backward Differentiation Formula Method for the Epitaxial Thin Film Equation with Slope Selection(1706.01943)

June 27, 2017 math.NA
In this paper, we study a novel second-order energy stable Backward Differentiation Formula (BDF) finite difference scheme for the epitaxial thin film equation with slope selection (SS). One major challenge for the higher oder in time temporal discretization is how to ensure an unconditional energy stability and an efficient numerical implementation. We propose a general framework for designing the higher order in time numerical scheme with unconditional energy stability by using the BDF method with constant coefficient stabilized terms. Based on the unconditional energy stability property, we derive an $L^\infty_h (0,T; H_{h}^2)$ stability for the numerical solution and provide an optimal the convergence analysis. To deal with the 4-Laplacian solver in an $L^{2}$ gradient flow at each time step, we apply an efficient preconditioned steepest descent algorithm and preconditioned nonlinear conjugate gradient algorithm to solve the corresponding nonlinear system. Various numerical simulations are present to demonstrate the stability and efficiency of the proposed schemes and slovers.
• ### Enhanced betatron radiation by steering a low-energy-spread electron beam in a deflected laser-driven plasma wiggler(1706.01033)

June 4, 2017 physics.plasm-ph
Laser wakefield accelerators (LWFA) hold great potential to produce high-quality high-energy electron beams (e beams) and simultaneously bright x-ray sources via betatron radiation, which are very promising for pump-probe study in ultrafast science. However, in order to obtain a high-quality e beam, electron injection and acceleration should be carefully manipulated, where a large oscillation amplitude has to be avoided and thus the emitted x-ray yield is limited. Here, we report a new scheme to experimentally enhance betatron radiation significantly both in photon yield and photon energy by separating electron injection and acceleration from manipulation of the e-beam transverse oscillation in the wake via introducing a slanted thin plasma refraction slab. Particle-in-cell simulations indicate that the e-beam transverse oscillation amplitude can be increased by more than 10 folds, after being steered into the deflected laser-driven wakefield due to refraction at the slab's boundaries. Spectral broadening of the x-rays can be suppressed owing to the small variation in the peak energy of the low-energy-spread e beam in a plasma wiggler regime. We demonstrate that the high-quality e-beam generation, refracting and wiggling can act as a whole to realize the concurrence of monoenergetic e beam and bright x-rays in a compact LWFA.
• ### Structure of nanoscale-pitch helical phases: blue phase and twist-bend nematic phase resolved by resonant soft X-ray scattering(1703.08007)

April 18, 2017 cond-mat.soft
Periodic structures of phases with orientational order of molecules, but homogenous electron density distribution: a short pitch cholesteric, blue phase and twist-bend nematic phase, were probed by a resonant soft x-ray scattering (RSoXS) at the carbon K-edge. The theoretical model shows that in case of a simple heliconical nematic structure two resonant signals corresponding to the full and half pitch band should be present, while only the full pitch band is observed in experiment. This suggests that the twist-bend nematic phase has complex structure with a double-helix, built of two interlocked, shifted helices. We confirm that the helical pitch in the twist-bend nematic phase is in a 10 nm range, for both, the chiral and achiral materials. We also show that the symmetry of a blue phase can unambiguously be determined through a resonant enhancement of x-ray diffraction signals, by including polarization effects, which are found to be an important indicator in phase structure determination.
• ### A novel energy-based phase field model for ferrodroplet deformation and breakup in a uniform magnetic field(1704.00645)

March 30, 2017 physics.flu-dyn
In this paper, we propose a novel, thermodynamically consistent phase field model to simulate the deformation and breakup of a ferrodroplet that is immersed in a viscous medium and subject to an applied uniform magnetic field. Instead of using the magnetic body force in the traditional Rosensweig model, the key idea of this model is to propose a new magnetic energy that enables direct effects of the magnetic field on the interface evolution. The model can thereby be derived from the variational principle via minimizing the free energy of the total system. This energy based modeling idea can be easily extended to include more external fields in a convenient and consistent way for more complicated applications involving multiple external fields. We validate the model by performing a series of numerical simulations, including the comparison with analytic solutions, the investigation of the effect of different types of magnetic fields, the dynamical behaviors of the ferrodroplet breakup under the strong magnetic field, magnetic/velocity/pressure field distributions, the magnetic energy density, and the inertial phenomenon, etc.
• ### Nanophotonic Lithium Niobate Electro-optic Modulators(1701.06470)

Modern communication networks require high performance and scalable electro-optic modulators that convert electrical signals to optical signals at high speed. Existing lithium niobate modulators have excellent performance but are bulky and prohibitively expensive to scale up. Here we demonstrate scalable and high-performance nanophotonic electro-optic modulators made of single-crystalline lithium niobate microring resonators and micro-Mach-Zehnder interferometers. We show a half-wave electro-optic modulation efficiency of 1.8V$\cdot$cm and data rates up to 40 Gbps.
• ### Second harmonic generation in nano-structured thin-film lithium niobate waveguides(1610.04197)

Jan. 3, 2017 physics.optics
Integrated thin-film lithium niobate platform has recently emerged as a promising candidate for next-generation, high-efficiency wavelength conversion systems that allow dense packaging and mass-production. Here we demonstrate efficient, phase-matched second harmonic generation in lithographically-defined thin-film lithium niobate waveguides with sub-micron dimensions. Both modal phase matching in fixed-width waveguides and quasi-phase matching in periodically grooved waveguides are theoretically proposed and experimentally demonstrated. Our low-loss (~ 3.0 dB/cm) nanowaveguides possess normalized conversion efficiencies as high as 41% W-1cm-2.
• ### Preconditioned Steepest Descent Methods for some Nonlinear Elliptic Equations Involving p-Laplacian Terms(1607.01475)

Nov. 24, 2016 math.NA
We describe and analyze preconditioned steepest descent (PSD) solvers for fourth and sixth-order nonlinear elliptic equations that include p-Laplacian terms on periodic domains in 2 and 3 dimensions. The highest and lowest order terms of the equations are constant-coefficient, positive linear operators, which suggests a natural preconditioning strategy. Such nonlinear elliptic equations often arise from time discretization of parabolic equations that model various biological and physical phenomena, in particular, liquid crystals, thin film epitaxial growth and phase transformations. The analyses of the schemes involve the characterization of the strictly convex energies associated with the equations. We first give a general framework for PSD in generic Hilbert spaces. Based on certain reasonable assumptions of the linear pre-conditioner, a geometric convergence rate is shown for the nonlinear PSD iteration. We then apply the general the theory to the fourth and sixth-order problems of interest, making use of Sobolev embedding and regularity results to confirm the appropriateness of our pre-conditioners for the regularized p-Lapacian problems. Our results include a sharper theoretical convergence result for p-Laplacian systems compared to what may be found in existing works. We demonstrate rigorously how to apply the theory in the finite dimensional setting using finite difference discretization methods. Numerical simulations for some important physical application problems -- including thin film epitaxy with slope selection and the square phase field crystal model -- are carried out to verify the efficiency of the scheme.
• ### Convergence Analysis and Numerical Implementation of a Second Order Numerical Scheme for the Three-Dimensional Phase Field Crystal Equation(1611.06288)

Nov. 19, 2016 math.NA
In this paper we analyze and implement a second-order-in-time numerical scheme for the three-dimensional phase field crystal (PFC) equation. The numerical scheme was proposed in [46], with the unique solvability and unconditional energy stability established. However, its convergence analysis remains open. We present a detailed convergence analysis in this article, in which the maximum norm estimate of the numerical solution over grid points plays an essential role. Moreover, we outline the detailed multigrid method to solve the highly nonlinear numerical scheme over a cubic domain, and various three-dimensional numerical results are presented, including the numerical convergence test, complexity test of the multigrid solver and the polycrystal growth simulation.
• ### A Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations(1611.02967)

Nov. 9, 2016 math.NA
We present a second-order-in-time finite difference scheme for the Cahn-Hilliard-Hele-Shaw equations. This numerical method is uniquely solvable and unconditionally energy stable. At each time step, this scheme leads to a system of nonlinear equations that can be efficiently solved by a nonlinear multigrid solver. Owing to the energy stability, we derive an $\ell^2 (0,T; H_h^3)$ stability of the numerical scheme. To overcome the difficulty associated with the convection term $\nabla \cdot (\phi \boldsymbol{u})$, we perform an $\ell^\infty (0,T; H_h^1)$ error estimate instead of the classical $\ell^\infty (0,T; \ell^2)$ one to obtain the optimal rate convergence analysis. In addition, various numerical simulations are carried out, which demonstrate the accuracy and efficiency of the proposed numerical scheme.
• ### Privately Information Sharing with Delusive Paths for Data Forwarding in Vehicular Networks(1502.05551)

Oct. 24, 2016 cs.OH
We discuss how to efficiently forward data in vehicular networks. Existing solutions do not make full use of trajectory planning of nearby vehicles, or social attributes. The development of onboard navigation system provides drivers some traveling route information. The main novelty of our approach is to envision sharing partial traveling information to the encountered vehicles for better service. Our data forwarding algorithm utilizes this lightweight information under the delusive paths privacy preservation together with the social community structure in vehicular networks. We assume that data transmission is carried by vehicles and road side units (RSUs), while cellular network manages and coordinates relevant global information. The approximate destination set is the set of RSUs that are often passed by the destination vehicle. RSU importance is raised by summing encounter ratios of RSUs in the same connected component. We first define a concept of space-time approachability which is derived from shared partial traveling route and encounter information. It describes the capability of a vehicle to advance messages toward destination. Then, we design a novel data forwarding algorithm, called approachability based algorithm, which combines the space-time approachability with the social community attribute in vehicular networks. We evaluate our approachability based algorithm on data sets from San Francisco Cabspotting and Shanghai Taxi Movement. Results show that the partially shared traveling information plays a positive role in data forwarding in vehicular networks. Approachability based data forwarding algorithm achieves a better performance than existing social based algorithms in vehicular networks.
• ### Convergence Analysis for Second Order Accurate Convex Splitting Schemes for the Periodic Nonlocal Allen-Cahn and Cahn-Hilliard Equations(1610.06973)

Oct. 21, 2016 math.NA
In this paper we provide a detailed convergence analysis for fully discrete second order (in both time and space) numerical schemes for nonlocal Allen-Cahn (nAC) and nonlocal Cahn-Hilliard (nCH) equations. The unconditional unique solvability and energy stability ensures $\ell^4$ stability. The convergence analysis for the nAC equation follows the standard procedure of consistency and stability estimate for the numerical error function. For the nCH equation, due to the complicated form of the nonlinear term, a careful expansion of its discrete gradient is undertaken and an $H^{-1}$ inner product estimate of this nonlinear numerical error is derived to establish convergence. In addition, an a-priori $W^{1,\infty}$ bound of the numerical solution at the discrete level is needed in the error estimate. Such a bound can be obtained by performing a higher order consistency analysis by using asymptotic expansions for the numerical solution. Following the technique originally proposed by Strang (e.g., 1964), instead of the standard comparison between the exact and numerical solutions, an error estimate between the numerical solution and the constructed approximate solution yields an $O( s^3 + h^4)$ convergence in $\ell^\infty (0, T; \ell^2)$ norm, which leads to the necessary bound under a standard constraint $s \le C h$. Here, we also prove convergence of the scheme in the maximum norm under the same constraint.
• ### Partial Procedural Geometric Model Fitting for Point Clouds(1610.04936)

Oct. 17, 2016 cs.CV, cs.GR
Geometric model fitting is a fundamental task in computer graphics and computer vision. However, most geometric model fitting methods are unable to fit an arbitrary geometric model (e.g. a surface with holes) to incomplete data, due to that the similarity metrics used in these methods are unable to measure the rigid partial similarity between arbitrary models. This paper hence proposes a novel rigid geometric similarity metric, which is able to measure both the full similarity and the partial similarity between arbitrary geometric models. The proposed metric enables us to perform partial procedural geometric model fitting (PPGMF). The task of PPGMF is to search a procedural geometric model space for the model rigidly similar to a query of non-complete point set. Models in the procedural model space are generated according to a set of parametric modeling rules. A typical query is a point cloud. PPGMF is very useful as it can be used to fit arbitrary geometric models to non-complete (incomplete, over-complete or hybrid-complete) point cloud data. For example, most laser scanning data is non-complete due to occlusion. Our PPGMF method uses Markov chain Monte Carlo technique to optimize the proposed similarity metric over the model space. To accelerate the optimization process, the method also employs a novel coarse-to-fine model dividing strategy to reject dissimilar models in advance. Our method has been demonstrated on a variety of geometric models and non-complete data. Experimental results show that the PPGMF method based on the proposed metric is able to fit non-complete data, while the method based on other metrics is unable. It is also shown that our method can be accelerated by several times via early rejection.
• ### Distributed Optimal Gas-Power Flow Using Convex Optimization and ADMM(1610.04681)

Oct. 15, 2016 math.OC
This paper proposes a convex optimization based distributed algorithm to solve multi-period optimal gas-power flow (OGPF) in coupled energy distribution systems. At the gas distribution system side, the non-convex Weymouth gas flow equations is convexified as quadratic constraints. The optimal gas flow (OGF) subproblem is solved by an iterative second-order cone programming procedure, whose efficiency is two orders of magnitudes higher than traditional nonlinear methods. A convex quadratic program based initiation scheme is suggested, which helps to find a high-quality starting point. At the power distribution system side, convex relaxation is performed on the non-convex branch flow equations, and the optimal power flow (OPF) subproblem gives rise to a second order cone program. Tightness is guaranteed by the radial topology. In the proposed distributed algorithm, OGF and OPF are solved independently, and coordinated by the alternating direction multiplier method (ADMM). Numerical results corroborate significant enhancements on computational robustness and efficiency compared with existing centralized OGPF methods.