In this paper, we study the leading-twist distribution amplitude (DA) of the heavy pseudoscalars (HPs), such as η c , η b , and B c , within the QCD theory in the background fields. New sum rules up to dimension-six condensates for both the HP decay constants and their leadingtwist DA moments are presented. From the sum rules for the HP decay constants, we obtain f η c = 453 ± 4 MeV, f B c = 498 ± 14 MeV, and f η b = 811 ± 34 MeV. Based on the sum rules for the HPs' leading-twist DA moments, we construct a new model for the η c , η b , and B c leadingtwist DAs. Our present HP DA model is also adaptable for the light pseudoscalar DAs, such as the pion and kaon DAs. Thus, it shall be applicable for a wide range of QCD exclusive processes. As an application, we apply the η c leading-twist DA to calculate the B c → η c transition form factor f B c →η c + (q 2 ). At the maximum recoil region, we obtain