def METHOD_NAME(c):
    '''
    Run circuit c on qubits initialized in the 0^n state and return the
    density matrix of the output state as qutip Qobj.
    '''
    state = qsim.sparsedm.SparseDM(c.get_qubit_names())
    c.apply_to(state)
    return qutip.Qobj(state.full_dm.to_array(), dims=[[2, 2], [2, 2]])