In Ref.[1], the Chern-Simons correction to the osculating orbital elements due to the dipole-dipole interaction for the scalar field was derived in Eq. (130) by integrating out the interaction Lagrangian, which was taken to be the kinetic term for the scalar field [Eq. (122) in [1]], using the scalar field solution. However, such kinetic term only leads to the left-hand side of the evolution equation for the scalar field, i.e.where J is the effective source term of the scalar field equation [2]. In order to obtain the correct interaction Lagrangian, one also needs to include the source term (L source ) that reproduces the right-hand side of the above evolution equation. The Lagrangian density for such an additional term is given by −ϑJ. Therefore, Eq. (126) needs to be changed towhere L source,int is the interaction part of the source term Lagrangian density. Notice that, by accounting for such additional term, the sign of the interaction Lagrangian now flips. This means that the sign in Eqs. (127), (129) and (130) needs to be reversed. The results in Sec. IX B are unaffected since we have neglected contribution from the dipole-dipole interaction because the spin of the secondary pulsar in the double binary pulsar is much smaller than that of the primary. We also report errors in App. B of [1]. In this appendix, we calculated dissipative corrections to the energy and angular momentum energy flux emitted from neutron star binaries based on calculations in [2]. However, some terms in the perturbed scalar and metric field equations were not properly accounted for in [2]. Taking these terms into account, we found that corrections to the gravitational energy and angular momentum flux due to metric perturbations enter at third post-Newtonian order relative to general relativity, which is of 1 post-Newtonian order higher than the effect due to the scalar field. Therefore, Apps. B.1.b and B.2.b in [1] are now irrelevant. We discuss this in more detail in [3].On the other hand, the scalar energy flux in App. B.1.a is affected as follows. First, the far-zone solution for the scalar field in Eq. (B3) becomesThen, the scalar energy flux in Eqs. (B5) and (B6) becomesrespectively. Finally, Eq. (B8) is corrected to T. Tanaka, Phys. Rev., D85, 064022 (2012), arXiv:1110.5950 [gr-qc].[3] K. Yagi, L. C. Stein, N. Yunes, and T. Tanaka, Phys. Rev., D93, 029902 (2016D93, 029902 ( ), arXiv:1110. We study isolated and binary neutron stars in dynamical Chern-Simons gravity. This theory modifies the Einstein-Hilbert action through the introduction of a dynamical scalar field coupled to the Pontryagin density. We here treat this theory as an effective model, working to leading order in the Chern-Simons coupling. We first construct isolated neutron star solutions in the slow-rotation expansion to quadratic order in spin. We find that isolated neutron stars acquire a scalar dipole charge that corrects its spin angular momentum to linear order in spin and corrects its mass and quadrupole moment to quadratic order in spin, as measured by ...