A theory for thermomechanical behavior of homogeneous DNA at thermal equilibrium predicts critical temperatures for denaturation under torque and stretch, phase diagrams for stable B-DNA, supercoiling, optimally stable torque, and the overstretching transition as force-induced DNA melting. Agreement with available single molecule manipulation experiments is excellent.PACS numbers: 87.14.gk, 87.15.ad, 87.15.Zg,64.60.De DNA is a highly refined nanomechanical object. The interplay between strong covalent bonds of the backbone and weak hydrogen interactions between bases [1], the thermal bath in which it is immersed, and the proximity of physiological conditions to the denaturation temperature, make DNA highly and non-linearly responsive to mechanical and thermal changes, and render any solely mechanical approach inviable. This is critical for nanotechnology, where DNA is the basis for novel materials [2], but understanding double helix thermomechanics would also illuminate biology, where enzymes involved in replication and repair are viewed as molecular motors.In the past twenty years, direct single molecule manipulation [3,4] has revolutionized our understanding of key aspects of DNA, revealing new couplings and transitions between different structures, whose nature and forms, however, are still speculative. When a few micrometer long strand of DNA is stretched to a tension of the order of pico-Newtons (pNs) to avoid formation of plectonemes, sharp transitions are activated at positive and negative torques, while an overstretching transition is observed for DNA under tension of 60 pN at zero torque [5]. Tentative tension-torque phase diagrams for the stability of B-DNA, and various phenomenological theories, have been proposed to explain these effects. The robustness of the Peyrard-Bishop-Dauxois approach (PBD) [6][7][8] has been corroborated by Cocco and Barbi [9, 10]: they incorporated torque and successfully reproduced denaturation by unwinding. However, these recent studies do not include tension, do not explain denaturation at overwinding, and do not provide phase diagrams in the tension-torque plane. Also, although much simpler than the molecular structure they describe, their complexity still calls for numerical treatments, and cannot offer analytic equations to more easily guide experiments.We address these issues by modeling the steric dependence of the base bond and the effect of tension and torque in a way suitable for elimination of the angular degrees of freedom via integration of the partition function, and obtain an effective energy for a PBD model which incorporates temperature and external loads. We compute the phase diagram for B-DNA and the dependence of supercoiling on torque, tension and temperature at criticality. Finally, we propose simple algebraic formulae for the observables.In our model i labels nucleotides separated by a distance a along the DNA backbone (Fig. 1), x i is the length of the i th base's bond, ω i = (θ i+1 − θ i−1 )/2 − Ω is the angular shift between nucleotides along th...